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Abstract 


1  A  model  Is  constructed  of  a  nonlinear,  one-dimensional  inhomogeneous  multi  - 
component  plasma  slab  which  has  the  characteristics  of  the  plasma  sheath  surround¬ 
ing  a  typical  hypersonic  re-entry  vehicle  at  200,  000  ft  while  traveling  at  about  18, 000 
ft /second.  A  plane  wave  incident  upon  the  slab  at  arbitrary  angle  can  induce  changes  in 

electron  temperature  without  affecting  the  neutral  gas  temperature.  The  variations 
of  electron  temperature  produce  changes  in  collision  frequency  and  the  rate  coeffi¬ 
cients  describing  the  various  electron  production  and  loss  mechanisms.  Such  changes 
cause  alterations  in  the  effective  dielectric  constant  of  the  medium.  On  the  basis  of 
a  model  of  a  plane-layered  medium  composed  of  a  stack  of  linear  homogeneous  sheets, 
the  reflection  and  transmission  coefficients  of  the  nonlinear  slab  may  be  computed 
by  a  step-by-step  numerical  integration  of  Maxwell's  equations  expressed  in  the  form 
of  difference  equations . 
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Nonfinaar  Transmission  Characteristics 
of  th«  Plasma  Shaath 

1.  DtTRODUCnON 

A  hypertonic  re-entry  vehicle  it  enveloped  by  a  shock-induced  sheath  of  ionized 
gas.  This  re-entry  plasma  sheath  can  greatly  influence  electromagnetic  communica¬ 
tions  frequencies  in  the  radio  and  microwave  range-  Many  calculations  have  been 
made  on  the  reflection  and  transmission  properties  of  the  pladma  sheath  surround¬ 
ing  a  hypersonic  re-entry  vehicle.  Perhaps  the  most  straightforward  analysis, 
which  can  yield  useful  information  on  signal  degradation,  is  based  upon  the  model 
consisting  of  a  plane  wave  Incident  on  a  dielectric  slab.  The  electromagnetic  prop¬ 
erties  of  the  re-entry  plasma  sheath  are  characterized  try  an  equivalent  dielectric 
slab  whose  complex  dielectric  constant  and  thickness  simulate  the  actual  plasma 
sheath.  However  this  type  of  analysis  is  not  capable  of  predicting  the  dominant 
characteristics  of  the  radiation  patterns  of  plasma-coated  slot  antennas  located  on 
a  re-entry  vehicle  surface. 

If  communications  between  the  vehicle  and  ground  are  effected  by  means  of  slots 
or  horns  flush-mounted  on  the  missile  surface,  a  suitable  model  describing  the  an¬ 
tenna  system  consists  of  a  slot-excited  plasma  slab  covering  a  ground  plane.  The 
prominent  features  of  the  radiation  patterns  of  such  a  structure  may  be  obtained  by 
evaluating  the  integral  representations  for  the  field  components  by  the  method  of 
steepest  descent.  Tamir  and  Oliner 1  have  demonstrated  that  the  pear  field  of  a, 
slot -excited  plasma  slab  is  represented  mainly  by  contributions  from  complex  waves. 

(Received  for  publication,  28  February  1963) 
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For  frequencies  below  the  plesma  frequency,  the  complex  waves  are  of  the  spectral 
type  which  do  not  transport  energy  but  represent  stored  energy.  At  frequencies 
above  the  plasma  frequency,  the  complex  waves  are  of  the  leaky-wave  type  and  may 
be  strongly  excited  for  certain  slab  thicknesses.  The  shape  of  the  far-fleld  radiation 
pattern  may  be  obtained  by  integrating  the  near  field  over  a  Kirchoff-Huygens  sur¬ 
face.  The  near-field  components  are  found  by  evaluating  the  integral  representa¬ 
tions  of  the  field  by  the  calculus  of  residues.  For  frequencies  greater  than  the 
plasma  frequency,  the  radiation  pattern  of  a  slot-excited  grounded  plasma  slab  exhibits 
peaks  at  the  critical  angle.  Shore  and  Meltz  have  investigated  the  far-field  compo¬ 
nents  using  Fourier  analysis  and  the  stationary  phase  procedure  for  the  case  of  a 
slot  in  a  ground  plane  covered  by  an  anisotropic  lossy  plasma  layer. 

These  methods  of  analysis  usually  assume  a  uniform  plasma  layer,  although 

inhomogenietles  in  one  dimension  may  be  taken  into  account  by  using  a  plane-layered 

3 

model.  Felsen  and  Marcuvitz  have  pointed  out  that  an  electromagnetic  field  prob¬ 
lem  with  prescribed  sources  may  be  reduced  by  suitable  modal  procedures  to  a  net¬ 
work  problem.  The  equivalent  network  problem  may  then  be  solved  for  a  typical 
mode  and  the  field  solution  is  then  obtained  by  modal  synthesis.  The  problem  of  the 
modal  analysis  of  a  plane-layered  medium  may  be  transformed  into  an  equivalent 
network  problem  of  a  transmission  line  consisting  of  uniform  sections,  where  the 
voltage  and  current  at  the  discontinuities  in  the  line  have  the  same  boundary  condi¬ 
tions  as  the  E  and  H  fields  have  at  the  discontinuities  in  the  plane-layered  medium. 

All  of  these  approaches  to  the  problem  of  electromagnetic  transmission  through 
the  re-entry  plasma  sheath  do  not  consider  the  effect  of  perturbations  induced  in  the 
plasma  medium  by  the  presence  of  an  electromagnetic  wave.  If  the  electromagnetic 
propagation  characteristics  of  the  plasma  medium  are  described  by  a  complex  die¬ 
lectric  constant  K,  then  this  parameter  may  be  related  to  the  electron  density  Ng 
and  electron-neutral  collision  frequency  r: 

K  .  Kr  +  JI^ 

(w_/u>)2  (v/w)(w  /u>)2 

■  1-  E  — T  -i - —K -  .  (1) 

i  +  <v/wr  i  +  (v/wr 

For  relatively  high  field  strengths,  it  is  expected  that  the  degree  of  ionization  and 
the  electron-neutral  particle  collision  frequency  will  change,  thereby  altering  the 
electromagnetic  transmission  characteristics  of  the  plasma  medium.  However, 
any  model  for  describing  such  a  nonlinear  Interaction  between  high-power  electro¬ 
magnetic  radiation  and  the  plasma  medium  will  have  to  include  a  detailed  picture  of 
the  interaction  of  the  local  fields  Inside  the  medium  with  the  plasma  constituents. 
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A  detailed  picture  of  the  various  microscopic  processes,  such  as  the  local  field- 
particle  interactions  and  the  particle-particle  interactions,  give  rise  to  a  medium  that 
may  be  characterized  on  a  macroscopic  level  by  a  complex  dielectric  constant  that 
is  dependent  upon  the  local  field  distribution  in  the  plasma  medium.  This  type  of 
medium  does  not  lend  itself  to  modal  analysis,  since  the  representation  of  an  arbi¬ 
trary  field  in  the  medium  by  a  Fourier  integral  over  the  wave  numbers  depends  upon 
the  validity  of  the  principle  of  superposition.  Bach  mode  is  assumed  to  propagate 
with  a  phase  velocity  that  depends  upon  the  dielectric  constant  of  the  medium,  and 
the  dielectric  constant  is  assumed  to  be  independent  of  the  presence  of  other  modes. 
This  condition  is  no  longer  satisfied  in  a  medium  whose  dielectric  constant  depends 
upon  the  local  field  distribution  in  the  medium.  For  this  reason  a  nonlinear  model 
was  constructed  consisting  of  a  plane  wave  of  arbitrary  angle  of  polarization  inci¬ 
dent  at  an  arbitrary  angle  upon  a  plane-layered  plasma  slab.  The  slab  consists  of 
all  the  constituents  normally  found  behind  a  high-temperature  shock  (5000*K).  The 
thermodynamic  and  gas  kinetic  parameters  were  selected  to  coincide  with  conditions 
found  behind  a  shock  of  a  representative  blunt-nosed  vehicle  traveling  at  18,  000  ft/ 
sec  at  an  altitude  of  200,000  feet. 

The  total  reflection  and  transmission  coefficients  of  a  plane-layered  plasma 
slab  having  the  same  characteristics  as  the  sheath  surrounding  a  typical  re-entry 
vehicle  may  be  computed  for  an  incident  plane  wave  consisting  of  a  single  frequency 
component.  The  nonlinear  response  of  the  plasma  medium  to  the  perturbing  electro¬ 
magnetic  wave  is  found  for  the  steady-state  case,  corresponding  to  an  Incident  plane 
wave  whose  frequency  is  much  greater  than  the  frequency  for  energy  transfer  be¬ 
tween  the  electron  gas  and  the  neutral-particle  gas.  The  condition  for  a  steady- 
state  nonlinear  response  of  a  plasma  to  a  perturbing  electromagnetic  wave  may  be 
written 

w  »Gt'  , 

where'  u  *  frequency  of  impressed  field 

G  *  relative  fractional  energy  loss  of  an  electron  in  colliding 
with  a  neutral  particle  (G  «  2m/M  for  electron  elastic 
collisions  with  neutrals  of  mass  M) 

v  *  electron-neutral  particle  collision  frequency. 

Since  tgjj  *  1/Gv  is  the  energy  relaxation  time  for  the  electron  gas,  this  condition 
may  be  written: 


where  T  *  2t/u  is  the  period  of  the  impressed  field.  Because  the  energy  relaxation 
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-  rap? 


time  for  the  electron  gas  is  much  greater  than  the  period  of  the  impressed  electro¬ 
magnetic  field,  the  electron  temperature  will  settle  down  to  some  average  value 
dependent  upon  the  mean  square  value  of  the  electromagnetic  field.  In  the  steady 
state  the  electromagnetic  propagation  characteristics  of  a  plasma  perturbed  by  a 
relatively  high- power  electromagnetic  field  may  be  determined  by  specifying  the 
complex  dielectric  constant,  which  is  a  function  of  the  following  parameters:  "* 

(1)  An  effective  electron  collision  frequency 


kEFF 


'  ^  (w^)  '  I vM  y4  exp{$r^)dv 


kEFF  '  *'max(Te)  +  3  Te  8T 


8k  (T  ) 
max1  e' 


(2) 


where 


max 


°V  ■ 


»(v)  f_„v  V  dv 
max 


■  it 


{f&r;}12  §vM  y2  exp(wr  )dv  ■ 


v(v)  is  the  electron  collision  frequency  as  a  function  of  electron  velocity  and 
Te  is  the  electron  temperature. 

(2)  The  Interelectron  collision  frequency 


v 


ee 


Ne  In  A _ 

11.4  A1/2Te3^2 


(3) 


where  (f  nA)  is  given  as  a  function  of  Ne  and  by  Spitser.5  A  *  1/1823  for 
electrons. 

(3)  The  electron  density  N#. 

(4)  The  Isotropic  part  of  the  electron  velocity  distribution  function  fQ  (which  is 
very  nearly  Maxwellian  even  in  the  presence  of  a  strong  electromagnetic  field  if 
the  interelectron  collision  frequency  is  much  greater  than  the  collision  frequency 
for  energy  transfer  between  the  electron  gas  and  the  neutral-particle  gas  v#e  »Gx). 

It  turns  out  that  the  condition  K#e  »  Ok  Is  satisfied  for  the  partially-ionised 
(0.  01  to  0. 1%)  plasma  behind  strong  shocks  in  air  over  a  considerable  range  of 


i 


s 


g 

pressures  And  temperature!.  Peekof f  hai  tabulated  the  electron  density,  the  inter- 
electron  collision  frequency,  and  the  electron-neutral  collision  frequency  for  the 
stagnation  region  behind  a  normal  shock  with  a  free-stream  velocity  of  18,  000  ft/sec 
at  altitudes  of  40,  60,  80,  and  100  km.  Peakoffs  tabulations  indicate  that  the  inter- 
electron  collision  frequency  (i>ee)  is  of  the  order  of  the  electron-neutral  collision 
frequency  (*  g);  1.  e.,  vgg  »  »>eN.  Since  the  fractional  energy  lose  of  an  electron 

per  collision  (G)  with  the  atoms  and  molecules  found  behind  a  high-temperature 
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shock  in  air  ranges  between  10  to  10  for  electron  energies  of  0.  5  to  2  ev  (cor¬ 
responding  to  electron  temperatures  of  5000*  to  20,  000*K),  the  condition  |;ee»G>'eN 
is  always  satisfied  behind  a  shock  with  gas  temperatures  at  5000TC  and  electron 
temperatures  of  5000*  to  20,  000*K  (see  Massey  and  Burhop,^  p.  279).  It  will  be 
demonstrated  in  the  next  section  (see  Ginzburg  and  Gurevich,  p.  132)  that  when  the 
condition  v  »  Gv  is  satisfied,  the  Isotropic  part  of  the  electron  velocity  distrlbu- 
tion  function  in  a  partially-ionized  plasma  is  Maxwellian,  even  in  the  presence  of  a 
strong  electromagnetic  field.  This  fact  permits  the  establishment  of  the  electron 
temperature  as  a  proper  thermodynamic  variable.  If  the  additional  assumption  is 
made  that  the  gas  temperature  is  not  changed  by  the  perturbing  electromagnetic 
field  (a  good  approximation),  then  the  steady-state  electron  density  may  be  found 
as  a  function  of  the  various  parameters  that  determine  the  rate  at  which  electrons 
appear  or  recombine.  The  rate  coefficients  for  electron  production  or  loss,  such 
as  the  rate  of  ionization  for  electron  Impact  on  neutrals,  dissociative  recombination, 
and  three-body  recombination  are  known  as  an  explicit  function  of  electron  tempera¬ 
ture  for  the  various  constituents  of  high-temperature  air. 

Since  the  predominant  electron  production  and  electron  loss  mechanisms  are 
believed  to  be  known  for  a  high-temperature  shock  in  air,  even  when  the  electron 
temperature  is  greater  than  the  neutral  and  heavy  ion  temperature,  the  computation 
of  the  change  in  the  electron  density  of  the  plasma  sheath  produced  by  changes  in 
the  electron  temperature  will  require  specific  knowledge  of  the  electron  temperature 
distribution  in  the  slab  of  plasma.  In  the  steady  state  the  electron  temperature  at 
each  point  in  the  plasma  slab  depends  upon  the  rate  at  which  the  electromagnetic 
field  heats  the  electron  gas  and  the  rate  at  which  the  electron  gas  16ses  energy.  The 
various  electron  energy  loss  mechanisms  include: 

(1)  Elastic  electron-neutral  and  electron-ion  collisions. 

(2)  Ionizing  collisions  by  electron  impact. 

(3)  Collisions  between  electrons  and  molecules  that  involve  the  excitation  of 
vibrational  or  rotational  states. 

(4)  Recombination  collisions  (dissociative  and  three  body). 

(5)  Heat  transport  due  to  gradients  of  the  difference  between  the  electron 
temperature  and  the  neutral  gas  temperature.  There  are  twa  heat  transport 
coefficients  (see  Anderson  and  Goldstein,8  p.  77).  One  transport  coefficient 
depends  upon  heat  conduction  due  to  electron-neutral  collisions  and  the  other 
depends  upon  heat  conduction  due  to  electron-electron  collisions  (Spitser  and 
Harm*"). 
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(6)  Heat  loss  through  electron  diffusion  which  involves  a  term  in  the  energy 
balance  equation  of  the  form  64k2  T. 

VN  •  VT. 
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On  the  basis  of  equilibrium-flow  calculations  performed  by  Rotman  and  Meltz1 1 
for  the  electron  density  and  temperature  profile  behind  a  high -temperature  shock  in 
air,  it  may  be  shown  that  the  various  heat  loss  mechanisms  of  the  electron  gas  due 
to  conduction  or  diffusion  are  negligible  in  comparison  with  losseB  due  to  elastic  and 
inelastic  electron-neutral  and  electron-ion  collisions.  The  model  assumed  consists 
of  a  plane  wave  incident  on  a  plasma  slab,  and  owing  to  the  infinite  extent  of  the  plane 
wave  front,  heat  losses  due  to  electrons  flowing  out  of  the  region  heated  by  the 
electromagnetic  field  will  not  be  considered.  These  preliminary  facts  permit  compu¬ 
tation  cf  the  total  reflection  and  transmission  coefficient  of  the  plasma  sheath  for 

o 

relatively  high-power  densities  (about  100  watts /cm  )  as  a  function  of  the  amplitude 
of  the  incident  plane  wave,  angle  of  incidence,  and  angle  of  polarization. 

One  of  the  first  investigations  of  the  nonlinear  interaction  of  electromagnetic 
radiation  with  the  plasma  sheath  surrounding  hypersonic  re-entry  vehicles  was  con¬ 
ducted  by  Sisco  and  Fiskin. 15  In  this  study,  the  assumption  was  made  that  only  the 
collision  frequency  of  the  hypersonically-produced  plasma  would  change  under 
perturbation  by  an  electromagnetic  wave  as  long  as  the  field  strengths  involved  were 
below  the  breakdown  threshold.  However  the  condition  for  breakdown  may  not  be 
sharply  defined  in  the  case  where  ionization  exists  prior  to  the  application  of  high- 
power  electromagnetic  radiation.  It  is  possible  that  the  electron  density  may  change 
by  a  factor  of  two  or  three  in  a  plasma  layer  surrounding  an  aperture  antenna  because 
of  changes  in  the  power  level.  Such  alterations  in  the  electron  density  will  have  as 
pronounced  an  effect  on  the  microwave  conductivity  of  the  plasma  as  changes  induced 

in  the  electron  collision  frequency.  At  power  levels  at  X-band  (10  kMc)  of  the  order 
2 

of  100  watts /cm  ,  it  is  estimated  that  the  electron  temperature  will  change  by  about 
a  factor  of  four  (from  5000*  to  20,  000*K).  This  will  produce  about  factor  of  two 
changes  in  electron  density  and  collision  frequency.  If  the  frequency  of  a  plane  wave 
incident  upon  a  plasma  slab  is  about  equal  to  the  average  plasma  frequency  in  the 
slab  (u  »  Up),  then  there  will  be  little  reflection  from  the  first  air-plasma  interface 
and  a  maximum  of  energy  coupling  between  the  electromagnetic  wave  and  the  lossy 
plasma.  Alterations  in  the  electron  density  and  collision  frequency  of  the  order  of 
a  factor  of  two  are  capable  of  inducing  50  to  100  percent  changes  in  reflection  coef¬ 
ficient  and  absorption  coefficient. 

The  method  of  computing  the  reflection  and  transmission  coefficients  of  the  non¬ 
linear  slab  of  plasma  proceeds  on  the  basis  of  a  step-by-step  numerical  integration 
of  Maxwell's  equations.  The  approximation  consists  in  replacing  the  inhomogeneous 
nonlinear  plasma  slab  with  a  set  of  homogeneous  linear  slabs.  If  Maxwell's  equations 
are  written  in  the  form  of  difference  equations,  and  the  step-by-step  numerical 
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Integration  is  accomplished  by  using  backward  differences  starting  on  the  face  of  the 
slab  from  which  the  transmitted  wave  emerges,  then  the  solution  for  the  field  dis¬ 
tribution  in  the  nonlinear,  nonhomogeneous  slab  may  be  determined  by  assuming  a 
value  for  the  amplitude  of  the  transmitted  wave.  In  this  scheme  the  fields  are 
completely  determined  at  the  point  Z-h,  once  the  fields  are  known  at  the  point  Z 
and  the  dielectric  constant  is  known  at  the  point  Z.  ' h'  is  the  increment  chosen  in 
obtaining  solutions  of  the  difference  equations  and  corresponds  to  a  layer  over  which 
the  dielectric  constant  is  assumed  not  to  vary.  The  dielectric  constant  at  the  point 
Z  in  the  plasma  slab  is  found  by  calculating  the  effective  collision  frequency  and 
the  electron  density  at  this  point  in  the  plasma.  These  parameters  are  determined 
once  the  electron  temperature  is  known.  The  electron  temperature  is  determined 
from  an  energy  balance  equation  for  the  layer  lying  between  the  points  Z  and  Z-h. 
The  energy  balance  equation  is  simply  a  statement  of  the  fact  that,  in  the  steady 
state,  the  rate  at  which  the  electromagnetic  field  supplies  energy  to  the  electron  gas 
is  equal  to  the  rate  at  which  the  electron  gas  loses  energy  to  the  neutral  and  ion  gas 
through  elastic  and  Inelastic  collisions. 

The  step-by-step  numerical  integration  of  the  electromagnetic  field  equations 
12  13  14 

has  been  discussed  by  Penico,  Stickler,  and  Richmond  for  the  case  of  a  plane 
wave  incident  upon  a  dielectric  slab  with  the  dielectric  constant  a  function  of  posi¬ 
tion  in  the  direction  perpendicular  to  the  interfaces.  Most  of  these  analyses  start 
with  the  wave  equation  for  the  field  components.  However  an  advantage  is  gained 

by  starting  from  Maxwell's  equations  directly,  particularly  in  the  nonlinear  problem. 

14 

Specifically,  Richmond  utilizes  the  wave  equation  for  the  field  components  E  and 
H.  In  his  analysis  the  solution  of  the  difference  equation  is  started  by  assuming  a 
value  for  the  amplitude  of  the  transmitted  field.  Since  the  wave  equation  is  of  the 
second  order,  the  field  must  be  known  at  a  second  point  to  uniquely  specify  the 
solution.  The  field  at  the  position  of  the  first  backward  difference  is  usually  ex¬ 
panded  in  a  Taylor  series  with  the  size  of  the  increment  (h)  as  parameter.  Deter¬ 
mination  of  the  field  at  the  position  of  the  first  backward  difference  by  a  Taylor  ex¬ 
pansion  is  not  necessary  if  Maxwell's  equations  (first-order  equations)  are  used 
directly.  It  should  be  noted  that  the  step-by-step  numerical  integration  of  the  elec¬ 
tromagnetic  field  equations  by  backward  differences  is  a  method  that  may  be  applied 
to  the  nonlinear  plasma  slab  problem  only  when  thermal  and  particle  diffusion  from 
layer  to  layer  can  be  neglected.  The  effects  of  the  boundaries  are  neglected  insofar 
as  it  is  assumed  that  the  temperature  and  electron  density  gradients  in  the  boundary 
layer  and  at  the  edge  of  the  shock  arc  so  steep  compared  to  a  wavelength  that  the 
boundaries  may  be  represented  as  plane  interfaces. 
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2.  THE  CONDUCTIVITY  07  A  PARTIALLY  IONIZED  PLASMA  IN  THE 

PRESENCE  07  A  STRONG  ELECTROMAGNETIC  FIELD 

If  an  electromagnetic  field  la  impreaeed  upon  a  partially-ionized  gas,  the  elec¬ 
trons  will  gain  energy  from  the  electromagnetic  field  by  suffering  collisions  with  the 
neutral  species  and  heavy  ions  in  such  a  manner  that  their  ordered  oscillatory  motion 
is  .changed  to  random  thermal  motion.  For  a  partially-ionized  plasma,  such  as 
high  temperature  air  (500010,  the  electrons  lose  only  a  small  relative  fraction  of 

their  energy  (G)  per  collision  with  neutrals  and  ions  even  for  electron  energies  up 
-3  -2 

to  2  ev(G*10  at  electron  energies  of  0.  5  ev  and  G*10  at  electron  energies  of 

2  ev).  Thus  it  is  a  good  approximation  to  assume  that  the  velocity  distribution 

function  for  the  neutral  atoms  and  molecules  and  the  ions  remains  Maxwellian  and 

that  their  temperature  remains  constant  even  in  the  presence  of  relatively  strong 

electromagnetic  fields. 

Initially,  in  the  absence  of  an  externally  impressed  electromagnetic  field,  the 
electrons  are  in  thermal  equilibrium  with  the  neutrals  and  positive  ions  behind  the 
high-temperature  shock  of  a  re-entry  vehicle  (this  is  true  within  about  100  mean 
free  paths  behind  the  shock  front).  Hence  the  electron  velocity  distribution  function 
f  is  Maxwellian  in  the  absence  of  an  electromagnetic  field  and  the  electron  tempera¬ 
ture  is  equal  to  the  gas  temperature  T  «  T.  Consider  a  c.  w.  electromagnetic 
plane  wave  Incident  upon  the  plasma  sheath  with  a  frequency  (u)  much  greater  than 
the  electron-neutral  collision  frequency  for  energy  transfer  (Gv»10  at  T#  *  5000TC, 
Gvarl0+7  at  T  ■  20,  000 *K  for  T*5000*K,  plpQ  ■  10*®  corresponding  to  a  shock  with 
free  stream  velocity  ■  18, 000  ft/sec  at  200,  000  feet).  In  this  case,  the  relaxation 
time  for  energy  transfer  between  the  electron  gas  and  neutral  gas  (ten>  is  much 
greater  than  the  period  of  the  impressed  electromagnetic  wave  rEN»  Tem,  sp  that 
the  average  energy  gained  by  the  electron  from  the  electromagnetic  field  will  depend 
upon  the  mean  square  value  of  the  field.  The  increase  in  the  average  energy  of  the 
electron,  resulting  from  randomization  of  the  electron’s  oscillatory  motion  through 
elastic  collisions  with  neutrals,  will  Induce  changes  in  the  various  electron  energy 
loss  processes  such  as: 

(1)  Collisions  which  can  excite  electronic  levels  of  atom#  and/or  molecules 

(2)  Collisions  which  can  excite  rotational  and  vibrational  levels  of  molecules 

(3)  Ionising  collisions  of  atoms  or  molecules 

(4)  Dissociative  recombination  of  ionic  molecules  and  three-body  recombination 
collisions  with  ionic  atoms. 

In  adAtion  to  the  elastic  and  various  types  of  Inelastic  collisions  of  the  electrons 
with  heavy  particles,  the  effect  of  electron-electron  collisions  must  also  be  considered. 
Electron-electron  collisions  result  in  only  a  small  relative  fractional  energy  loss 
per  collision 
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ee 
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Thia  ia  a  consequence  of  the  fact  .that  the  total  collision  cross  section  for  electron- 
electron  collisions  is  obtained  by  integrating  the  Rutherford  differential  scattering 
cross  over  all  angles  down  to  a  minimum  scattering  angle,  corresponding  to  a  maxi- 
mum  impact  parameter  equal  to  the  Debye  length  (see  Ginzburg,  p.  129). 

r  ™,  i1,! . 
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The  fractional  momentum  loss  for  electron- electron  collisions  is  of  the  same  order 
as  the  fractional  energy  loss  for  electron-electron  collisions.  Spitzer5  gives  the 
following  expression  (Eq.  (3)]  for  the  interelectron  collision  frequency 
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where  A  ■  1/1823  for  electrons  and  I  n  A  ia  tabulated  as  a  function  of  Ne  and  Te- 

For  a  fairly  wide  range  of  re-entry  conditions,  especially  in  and  near  the  stag¬ 
nation  region  of  vehicles  traveling  at  18,  000  ft /sec  or  greater  at  altitudes  between 
300,  000  and  100, 000  ft,  the  interelectron  collision  frequency  (i»ee)  is  several  orders 
of  magnitude  greater  than  the  collision  frequency  for  energy  transfer  between  the 
electron  gas  and  the  gas  of  heavy  particles  (</ee>>  Gv).  When  this  condition  prevails, 
any  perturbation  such  as  an  electromagnetic  wave  applied  to  the  partially-ionized 
plasma  will  result  in  the  redistribution  of  energy  among  the  electrons  in  a  time 
r _ •  l/v _ which  is  much  faster  than  the  time  for  energy  transfer  between  the 

M  06 

electron  gas  and  the  heavy  particle  gas  (l/Gv). 

Such  an  effect  will  have  a  pronounced  influence  on  the  response  of  a  plasma  to 
a  dynamic  perturbation.  The  dynamic  (time  dependant)  response  of  a  re-entry  plasma 
sheath  under  the  Influence  of  high-power  elect  rdmagnetic radiation  has  been  care- 
hilly  analysed  by  King. 18  King's  investigation  proceeds  from  the  WKB  solution  of  the 
wave  equation  for  an  electromagnetic  wave  incident  upon  a  semi -infinite  plasma. 

The' validity  of  the  WKB  solution  is  restricted  to  cases  where  the  changes  in  electron 
density  (or  dielectric  constant  K)  are  small  compared  to  a  wavelength  dK/dZ<<K/\. 
The  WKB  method  was  developed  in  connection  with  obtaining  solutions  to  the 
Schroedlngsr  equation.  In  this  class  of  problems  the  DeBroglle  wavelength  is  almost 
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Always  very  small  compared  to  distances  over  which  the  potential  function  varies. 
Also,  with  reference  to  ionospheric  propagation,  the  criterion  for  the  validity  of  the 
WKB  approximation  is  met  for  frequencies  in  the  megacycle  range  and  higher.  How¬ 
ever  the  criterion  may  not  always  be  met  for  electromagnetic  propagation  through 
the  re-entry  sheath,  where  the  electron  density  may  vary  by  ten  orders  of  magnitude 
over  a  distance  of  about  1  cm  with  wavelengths  in  the  centimeter  range  (Rotman  and 
Meltz11).  The  WKB  solution  contains  no  Information  on  the  amplitude  of  the  reflect¬ 
ed  wave,  because  reflection  is  neglected  in  the  approximation  of  geometric  optics. 

17 

Ginztmrg  has  indicated  a  modification  of  the  WKB  method  to  Include  the  effects  of 
weak  reflection. 

16 

The  starting  point  in  King's  analysis  is  the  WKB  solution  for  the  electro¬ 
magnetic  field  distribution  in  the  re-entry  plasma  sheath.  The  WKB  solution  is 
written  in  the  form 


e  -Eo*/yk 


exp  ± 


where 


k2  *ko(f/‘o  '  ia/u,€o) 


kQ  *  u/c  *  2*/Xq 


(8) 


m(w2+y2) 

a  is  the  real  part  of  the  conductivity  and  the  propagation  direction  in  the  direction 
of  the  electron  density  gradient  (x-direction).  This  form  for  the  field  distribution 
is  then  substituted  into  the  expression  for  energy  deposition  per  unit  time  in  the 
plasma  (a  E2). 

The  fraction  of  the  total  energy  deposited  in  the  plasma  which  is  responsible 

o 

for  producing  increased  ionization  is  given  by  fcE  where  f  *  ojo  and  Oj  is  the 
ionization  cross  section  while  a  is  the  total  cross  section.  If  this  quantity  is  then 
divided  by  the  energy  required  to  form  an  ion  pair  in  the  plasma,  «,  the  result  yields 
the  local  rate  of  growth  of  electrons  due  to  the  electromagnetic  field: 
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“St®  *  •£ ffE  •  (7) 

The  electron  density  appears  explicitly  on  the  left-hand  aide  and  implicitly  on  the 

16 

right-hand  side  of  this  equation.  King  has  obtained  analytic  expressions  for  the 
time  dependence  of  density  by  making  a  suitable  transformation  of  variables.  The 
effects  of  diffusion  and  electron  attachment  may  also  be  taken  into  account,  although 
the  solutions  may  no  longer  be  expressed  in  terms  of  known  functions.  The  Implicit 
assumption  of  such  an  approach  rests  upon  the  relation. 
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where  K  is  the  effective  dielectric  constant 
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For  the  case  of  a  medium  with  a  time-dependent  dielectric  constant.  Maxwell's 
equations  are: 


curl  E 


curl  ft 
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This  leads,  in  general,  to  a  wave  equation  of  the  form: 
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The  WKB  approximation  is  a  valid  solution  of  the  wave  equation  only  if 


w  » 


91  nN. 
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which  permit*  the  term* 


to  be  neglected  end  e  solution  to  be  obtained  in  the  form  E  «  Aft)  e^w*  .  The 
condition  u  »  SinN./8t  implies  that  the  characteristic  time  for  electron  density 
buildup  is  much  longer  than  the  period  of  the  Impressed  field.  This  condition  Is  al¬ 
most  always  satisfied  for  most  plasma-microwave  Interactions.  A  final  point  should 
be  made  that  no  provision  for  changes  in  the  collision  frequency  have  been  Included 
in  the  analysis,  although  such  changes  may  generally  be  overshadowed  by  changes 
in  the  electron  density. 

Only  the  steady-state  response  of  a  re-entry  plasma  sheath  to  high-power  electro¬ 
magnetic  radiation  will  be  considered  in  this  report.  The  question  of  the  relative 
importance  of  the  relaxation  time  for  equilibration  of  energy  among  the  electrons 
(l/i'ee)  and  the  relaxation  time  for  equilibration  of  energy  between  the  electrons  and 
neutrals  and  ion  (1/Gv)  is  Important  from  the  point  of  view  of  the  dynamic  response 
of  the  plasma  to  the  perturbation  of  an  electromagnetic  wave.  However,  even  in  the 
steady  state,  electron-electron  collisions  can  p’»y  a  significant  role  in  determining 
the  form  of  the  Isotropic  part  of  the  electron  velocity  distribution  function.  The 
general  form  for  the  isotropic  part  of  the  electron  distribution  function  in  a  partially- 
ionised  plasma  in  the  presence  of  a  strong  electromagnetic  field  is  not  always  given 
by  the  Margenau18  expression 


fo  ‘  c  e*P 


mvdv 
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2e2E2 


3mG(K2+w2)J 


(9) 


Where  C  ■  a  normalizing  factor 
T  ■  gas  temperature 

v  ■  electron-neutral  collision  frequency  (a  function  of  electron  velocity). 

The  Margenau  expression  for  the  electron  distribution  function  was  derived 

under  the  assumption  that  electron-electron  collisions  are  negligible  (v«  Gv). 

10  8  " 

However,  as  pointed  out  by  Cahn  and  by  Ginzburg  and  Gurevich,  when  the  inter¬ 
electron  collision  frequency  is  of  the  order  of  the  collision  frequency  for  energy 
transfer  between  electrons  and  neutrals  (veew  Gv),  the  form  of  the  Isotropic  part  of 
the  electron  velocity  distribution  function  is  given  by  a  rather  complex  expression 
which  reduces  to  the  Margenau  expression  when  v#e  «  Gv  and  to  the  Maxwellian 
form  when  v##  »  Gv. 
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i.  1  Hie  Boltmann  Equation  for  a  Partially  Ionised  Gas 


The  propagation  characteristics  of  electromagnetic  waves  in  a  finite  layer  of 
plasma  may  be  described  by  a  complex  dielectric  constant  (or  conductivity)  which  is 
a  scalar  quantity  in  the  absence  of  a  magnetic  field.  The  conductivity  (o)  is  com¬ 
puted  by  taking  the  first  velocity  moment  of  the  electron  velocity  distribution  function: 


o  « 


l^f^v 


(10) 
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Ginzburg  and  Molmud  have  shown  that  Eq.  (10)  for  the  complex  conductivity 
reduces  to  the  standard  expression 


0  *  ar  *  1  ai 
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only  if  the  collision  frequency  in  Eq.  (11)  is  replaced  by  an  'effective'  collision 
frequency  (•'gyp)  given  by  Eq.  (21  The  effective  collision  frequency,  as  given  by 
Eq.  (2),  is  valid  only  if  u2  »  v2  .  For  the  low-frequency  limit  (w2«  v2),  the 
effective  collision  frequency  is  defined  differently.  At  intermediate  frequencies,  the 
right-hand  side  of  Eq.  (11)  must  be  multiplied  by  a  parameter  that  is  a  function  of 
<*>/« EFF'  To  utilize  the  exact  Expression  (10)  for  the  conductivity,  the  form  of  the 
electron  velocity  distribution  function  must  be  found  by  solving  the  Boltzmann  equa¬ 
tion. 

The  hill  description  of  a  multicomponent  system  such  as  a  partially-ionized 
plasma  Involves  a  series  of  coupled  Boltzmann  equations  for  the  distribution  function 
of  each  component. 


(12) 
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where  ((,„_)  ie  the  distribution  (unction  for  the  I th  ion  constituent,  f  .  is  the  distri- 
ion/  th  nj 

button  function  for  the  j  neutral  constituent,  the  subscript  e  refers  to  electrons,  and 
the  subscript  n  refers  to  neutrals.  Bmn  represents  the  contribution  to  the  collision 
integral  for  the  distribution  function  of  the  m4*1  type  of  particles  in  collision  with  an 
nt4>  type.  The  Boltzmann  equations  for  the  electron,  ion,  and  neutral  components  are 
coupled  through  the  collision  integral  terms.  The  usual  assumption  made  in  investi¬ 
gations  of  the  interaction  of  radio  frequency  or  microwave  energy  with  a  partially- 
ionized  gas  is  that  only  the  electron  distribution  function  changes,  and  the  neutral 
and  ion  distribution  functions  are  assumed  to  remain  Maxwellian  at  a  constant  tem¬ 
perature.  This  assumption  is  valid  when  the  collision  frequency  is  relatively  small 
compared  to  the  plasma  frequency,  the  plasma  is  less  than  1%  ionized  and  when 

the  fractional  energy  loss  of  an  electron  per  collision  with  a  neutral  is  small 
-2 

(G  <  10  for  electron  energies  below  2.  5  ev  for  most  gases). 

The  macroscopic  fields  that  appear  in  Eq.  (12)  are  self-consistent,  so  that  the 
electrons  do  not  interact  directly  with  one  another  on  a  microscopic  level  (except  for 
distances  less  than  the  Debye  length,  f  D).  Instead,  the  ensemble  of  electrons  give 
rise  to  a  macroscopic  electromagnetic  field  which  then  can  interact  with  the  individual 
electrons.  The  macroscopic  fields  that  appear  on  the  left-hand  side  of  the  Boltzmann 
equation  remain  constant  during  a  collision.  This  implies  that,  when  considering 

electron-electron  collisions,  the  sphere  of  interaction  must  be  less  than  a  Debye 
21 

length  (see  Drummond,  p.  12). 

The  great  majority  of  papers  on  the  kinetic  theory  of  electromagnetic  wave 
interaction  with  a  partially-ionized  gas  neglect  ion-electron  and  electron-electron 

collisions  (Holstein,22  Margenau,1®'  22  Reder  and  Brown,24  Allis  and  Brown2®). 
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Cahn  has  discussed  the  effects  of  electron-electron  and  electron-neutral  collisions 

g 

on  the  distribution  function.  Ginzburg  and  Gurevich  have  clearly  indicated  the 
conditions  under  which  it  is  permissible  to  neglect  electron-electron  and  electron- 
ion  collisions  in  a  partially-ionized  gas.  Elastic  collisions  between  electrons  and 


neutrals  may  be  treated  in  the  most  straight-forward  manner.  The  expression  for 

the  collision  integral  for  electron-neutral  collisions  is  obtained  in  terms  of  the 

neutral  particle  density,  temperature,  elastic  collision  cross  section,  and  the  zero- 

order  electron  energy  distribution  function.  The  derivation  of  these  expressions 

usually  involves  the  assumption  that  the  electron  suffers  only  a  small  chahge  in 

energy  on  collision  with  a  neutral  (which  is  not  necessarily  true  for  excitation  or 
26 

ionizing  collisions).  Allis  has  treated  the  case  where  account  is  taken  of  recoil 
of  the  heavy  molecule  under  elastic  collisions  with  electrons.  Inelastic  collisions 
between  electrons  and  neutrals  are  accompanied  by  the  excitation  of  rotational,  vi¬ 
brational,  or  electronic  levels  and  also  by  ionization  and  recombination.  In  addition. 


charge  exchange  between  ions  and  second -order  Impacts  are  possible,  in  which  the 
energy  of  an  excited  state  of  the  molecule  is  transferred  to  incoming  electrons. 


IS 


However,  an  exact  calculation  taking  account  of  all  these  Inelastic  processes  would 
be  exceedingly  complex.  The  cross  sections  for  many  of  these  processes  are  known 
exactly  only  in  a  few  cases.  In  this  report,  the  effects  of  charge  exchange  and 
second-order  impacts  (inelastic  collisions  of  the  second  kind)  are  neglected. 

The  Boltzmann  equation  for  the  electron  velocity  distribution  function  for  a 
partially-ionized  plasma  (the  re-entry  sheath)  in  the  absence  of  a  magnetic  field 
may  be  written: 


f  V  +  S*-  V"Bee+?Bei+f  1^  *  (13> 

where  E  is  the  amplitude  of  the  externally  impressed  electromagnetic  wave  and  the 
forces  exerted  on  the  electrons  by  the  magnetic  component  of  the  wave  are  neglected. 
Since  diffusion  is  neglected  in  comparison  with  other  electron  Iobb  mechanisms,  the 
Boltzmann  equation  may  be  solved  for  the  case  of  a  homogeneous  plasma  within  each 
section  of  the  plane-layered  medium.  Hence,  the  Boltzmann  equation  may  be  written: 

f  +  3lSV’Bee+£Bei+fBenj  ‘  (14) 

Ginzburg  and  Gurevich®  have  shown  that  when  the  electron  energy  loss  parameter  G 
is  much  less  than  one,  a  spherical  harmonic  expansion  of  the  electron  distribution 
function  converges  rapidly,  so  thai  the  expansion  may  be  terminated  after  the  first 
two  terms 


f(r ,  v,t) 


? 


f|  (r ,  v,t)  Pj  (cos  o)  *  fQ  +  fj  cos  a  , 


(15) 


where  the  electric  field  vector  is  taken  to  lie  along  the  Z-axis  and  a  is  the  angle 
between  2  and  v . 

Using  the  spherical  harmonic  expansion  [  Eq.  (15)J  ,  together  with  the  expres¬ 
sion  for  partial  derivatives 


(16) 


e2  « 

m  *  v 


♦  5 


sin3  u 


8P. 


1_ 


8(cos  a) 


one  obtains 


IS 


l 


si  r_L_  J-V2f  +  (* a 

m  ^y2  dv  v  X1  ^  dv 


H 


_1  „s, 

dv  2 


cob  a  + 


(17) 


If  Eq.  (17)  and  Eq.  (IS)  are  substituted  into  Eq.  (14),  and  if  then  the  resulting 
equation  is  multiplied  successively  by  PQ  (cos  a),  Pj  (cos  a),  P2  (cos  a)  etc. ,  and 
integrated  over  2*  sin  a  d  a ,  the  resulting  chain  of  coupled  equations  is  obtained 


where 


is1  +  [i* +  ^5  *  (v%)]  *  Bi 

"if  +  [f  v£  (-^)+  *  (y4f a)]  “ b2 

Bl  ■?irJ,[Bee+?B«i+5BCnJ-pi  <c°"“>dQ 


and 


(18a) 

(18b) 

(18c) 


dQ  ■  2s  sin  a  da 

-2 

Since  the  electron  energies  will  be  considered  only  up  to  2  ev,  for  which  G  »  10 
for  each  of  the  constituents  of  high -temperature  air,  only  the  equations  for  fQ  and 
f j  need  be  considered* 


2.  1.  1  COLLISION  INTEGRAL  FOR  ELASTIC  COLLISIONS  BETWEEN 
ELECTRONS  AND  NEUTRALS  OR  IONS 

The  collision  Integral  may  be  evaluated  by  considering  the  rate  of  change  of 
the  number  of  particles  in  a  velocity  volume  element  d3v.  This  rate  of  change  may 
be  found  by  taking  the  difference  between  the  number  of  particles  scattering  into 
that  volume  element  and  the  number  scattered  out: 

*  Actually,  G  •  10* 1  at  2  ev  for  NO,  but  the  NO  concentration  in  a  shock  at  5000*K 
is  less  than  or  equal  to  1/10  the  concentration  of  N2  and  02> 


t 


I 
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BdSv  «  (V  (^')f(v')c  o  (6,  c)  d2Ud3V’d3v 

-  rF(^)f(?>ca(e.c)d2Qd3Vd3v  , 


F(V)  *  distribution  function  for  scattering  particles 
f(v)  *  distribution  function  for  scattered  particles 
c  *  relative  speed  of  particles 

a  >  differential  scattering  cross  section 
d2Q  •  2s  sin  6  de 

0  *  scattering  angle  ■  angle  between  v  -  V  and  v 1  -  ty. 

The  primes  refer  to  the  quantities  before  a  collision  and  the  unprlmed  quantities 
refer  to  after  a  collision.  By  Liouville*  s  theorem 

d3v'd3V'  «  d3vd3V  . 

*  . 

so  that  the  collision  integral  may  be  written: 

B  »  (* [F(V‘ )  Uv  ')  -  F(Q)  f(v»  c  ed20d3V.  (20) 

%) 
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Allis  has  indicated  how  the  spherical  harmonic  components  of  the  collision  integral 
may  be  evaluated.  If  the  spherical  harmonic  expansion  for  f(v)  [  Eq.  (15)]  is  sub¬ 
stituted  into  Eq.  (20),  and  use  is  made  of  the  addition  theorem  for  spherical  har¬ 
monics, 

t(v')  *  £  tt  (V)  jp4  (a)PJ(e)  +  2  £  {{  ■  P|m(0)cos  m(o-*j  ,  (21) 

where  *  ■  azimuthal  angle,  then  the  Boltzmann  collision  integral  may  be  written 
B  »  T.  B.  P,  (cos  o) 


Bj  -  toY  F(V')  tf  <v'>  P<  (0)  •  cffdaQd3V'  (v)ced20d3V. 


where 
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The  lero-order  term  In  the  spherical  harmonic  expanaion  of  the  electron-heavy 
particle  collision  integral  representing  elastic  collisions  may  be  written 


B 


elas. 

o 


2m  1 
*  M+m  v 


2v 


8 

W 


(23) 


where 


m  *  electron  mass 

M  *  heavy  particle  mass 

T  ■  temperature  of  heavy  particle  gas 

G  *  2m /M+m  *  relative  fraction  of  the  energy  loss  of  an 
electron  per  collision  with  heavy  particle 

v  ■  N  v  Cff(v,  9)(l-coa  0)2rsinOdO. 
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Allis  has  derived  Eq.  (23)  on  the  basis  of  the  assumption  that  the  electron  loses 
only  a  small  fraction  of  its  energy  in  colliding  with  a  heavy  neutral  particle  (G«l). 
Also,  it  is  assumed  that  the  molecules  recoil  under  electronic  Impact  and  they  pos¬ 
sess  a  Maxwellian  distribution  in  velocity. 

The  first  spherical  harmonic  component  of  the  collision  integral  representing 
elastic  electron-neutral  collisions  may  be  immediately  derived  by  using  the  fact  that 
the  magnitude  of  the  electron's  velocity  changes  only  slightly  during  an  elastic  col¬ 
lision  with  a  heavy  particle  (v*  ■  v  and  V  *  V),  so  that 


B*“-  -  ^a(e.c)  c{fj(v')  P(V')  Pj(coa  9)  -  fj(v)  F(V)}  d20d3V 


-fl(v)J' 


d3Vd2OevF(V)(l-cos  9) 


*  -v(v)  fj{v) 


(24) 


where 


v  a  Nv^»(v,  ©Hl-cos  ©)  d80 

where  N  ■  neutral  particle  density  and  s  is  the  collision  frequency  for  momentum 
transfer. 
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2.  1.  2  COLLISION  INTEGRAL  FOR  INELASTIC  COLLISION  BETWEEN 
ELECTRONS  AND  NEUTRALS 

Inelastic  collisions  between  electrons  and  neutral  particles  result  in  the  exci¬ 
tation  of  rotational,  vibrational,  and/or  optical  levels  of  molecules.  For  a  partially- 
ionized  plasma  consisting  of  diatomic  gaseB  such  as  oxygen,  nitrogen,  and  NO,  the 
dominant  electron  energy  loss  mechanisms  for  electrons  of  about  1  ev  or  lower  en¬ 
ergy  are  excitation  of  rotational  and  vibrational  levels.  The  excitation  of  rotational 

-2  -4 

levels  results  in  an  energy  loss  of  10  to  10  ev;  whereas  the  excitation  of  vibra¬ 
tional  levels  results  in  electron  energy  losses  of  about  0.  1  to  0.  5  ev.  If  the  electron 
energy  distribution  function  is  Maxwellian,  then  ionizing  collisions  can  become  im¬ 
portant  when  the  average  electron  energy  is  1/10  or  greater  then  the  ionization  po¬ 
tential  of  any  atomic  or  molecular  constituent  (NO  has  the  lowest  ionization  potential, 

9.  25  ev,  of  any  constituent  in  high-temperature  air).  This  is  due  to  the  long  tail  at 
high  energies  of  the  Maxwellian  distribution  function. 

In  the  present  report,  electromagnetic  fields  incident  on  the  re-entry  sheath 
will  be  considered  with  intensities  sufficient  to  raise  the  electron  temperature  from 
5,  000*K  (0.  5  ev)  up  to  20,  000*K  (2  ev).  Measurements  of  the  mean  energy  of  electrons 
in  swarm  tube  experiments  together  with  Luxembourg  (cross  modulation)  type  ex¬ 
periments  make  it  possible  to  deduce  the  energy  loss  for  electrons  in  this  energy 

07  28 

range.  Such  experiments  have  been  performed  by  Healey  and  Heed,  Harries, 

28  30  7 

Haas,  and  Huxley.  Table  1  (from  Massey  and  Burhop  )  gives  the  fractional 

energy  loss  of  electrons  \  for  the  various  constituents  of  air  for  electron  energies 

between  0  and  8.  0  ev. 

Since,  for  an  average  electron,  only  a  small  part  of  the  energy  is  lost  below 
and  up  to  excitation  of  rotational,  vibrational,  and  optical  levels,  the  integral  for 
such  inelastic  collisions  in  a  molecular  plasma  may  be  represented  in  the  form 

g 

(Ginzburg  and  Gurevich  ); 

■>r~  •  j.  &  {w1-1  •*  («  %  *  4 

where  GEFF(V)  describes  the  total  relative  fractional  energy  loss  of  an  electron  per 
collision  due  to  excitation  of  rotational,  vibrational  and  optical  levels. 

vinel  ■  N  v  J  <rlnel  (v,  9)  (1-cos  0)  d20 

and  *  cross  section  for  such  inelastic  collisions.  It  may  be  noted  that  the 
relative  fractional  energy  loss  parameter  (G)  is  known  from  experiment  only  as  a 
function  of  electron  temperature  fT#),  and  not  electron  velocity. 
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Here,  X  »  G(1  -  4eo/3e) 

«_  *  mean  energy  of  gas  molecules 
eQ  =  1.  3  X  10‘3  ev 

e  »  mean  energy  of  Impacting  electron 
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However  this  amount  of  information  ia  sufficient  for  the  present  investigation,  since 
the  Boltsmann  equation  [  Eq.  (18))  will  be  integrated  over  the  zero-order  distribution 
function,  which  is  Maxwellian.  If  cases  are  to  be  considered  where  the  interelectron 
collision  frequency  (vee)  is  not  sufficiently  high  compared  with  the  collision  frequency 
for  energy  transfer  (Gv)  to  maintain  a  Maxwellian  form  for  the  isotropic  part  of  the 
distribution  function,  then  G  must  be  known  as  a  function  of  electron  velocity  v  ■ 
This  functional  dependence  may  be  obtained  by  solving  the  integral  equation  for  G(v): 


Cg<v>  v(v)v4exp  [_■$£] 


dv  » 


GfT e>';EFF(Te)  (26) 


where  the  right-hand  side  of  Eq.  (26)  is  known. 

The  first  spherical  harmonic  component  of  the  inelastic  collision  Integral  repre¬ 
senting  electron  impact  collisions  which  excite  rotational,  vibrational,  and  optical 
levels  in  the  molecular  plasma  may  be  written  in  the  form: 

B  gneiss  ,  _  ylnelaSf  ^  (27) 

where 


ylnelas  .  N  v  (*  (v#  g)  (i  .  COb  0)  d2Q  , 

and  <rinel(v,  0)  is  the  differential- scattering  cross  section  for  inelastic  scattering. 


2.  1.  3  COLLISION  INTEGRAL  FOR  ELASTIC  SCATTERING  BETWEEN 
ELECTRONS  AND  POSITIVE  IONS 

The  primary  assumptions  made  in  the  derivation  of  Eqs.  (23)  and  (24)  for  elastic 
collisions  between  electrons  and  heavy  neutral  particles  were:  (a)  the  fractional 
energy  loss  of  the  electron  per  collision  is  small,  (b)  the  electron  mass  is  very  small 
compared  with  the  neutral  mass  m  «  M,  and  (c)  the  temperature  of  the  gas  con¬ 
sisting  of  the  heavy  particles  remains  constant  regardless  of  the  average  electron 
energy.  Hence,  Eqs.  (23)  and  (24)  for  the  zero-  and  first-order  spherical  harmonic 
components  of  the  collision  integral  may  also  be  used  to  describe  elastic  collisions 
between  ions  and  electrons  when  the  appropriate  expression  for  the  ton-electron 
collision  frequency  is  used.  The  equation 

|  ^on-Nionvf»<v.8)(l-cos0)d2Q 

j 

j 

I 

tr 

\ 
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governs  the  functional  dependence  of  the  electron-ion  colliclon  frequency  on  electron 
velocity,  where  c(v,  6)  ie  given  hy  the  Rutherford  formula: 

..Mr)’  — J— 

v.2mv  /  sin  9/2 


where  6  »  scattering  angle  and 

Nj  «  particle  density  of  ions  •  Nfi. 

Thus, 


where  9^  la  the  minimum  angle  of  acattering  corresponding  to  the  maximum  im¬ 
pact  parameter.  Since,  as  shown  fay  Drummond,21  collisions  between  charged  par¬ 
ticles  which  occur  at  distances  greater  than  the  Debye  length  are  described  by  macro¬ 
scopic  fields  which  appear  on  the  left-hand  side  of  the  Boltzmann  equation,  the  maxi¬ 
mum  Impact  parameter  will  be  taken  equal  to  the  Debye  length  (<D>  Eq.  (5)]  : 


<Dtan9min/2 


(29) 


or 


W2Un 


■1  e2  . 

a  "■  ■  ■  i 

mv  <  D  mvzl  D 


Substituting  Expression  (29)  for  9mln  into  Eq.  (28)  yields: 


vion(v)  *  2*Ne 


In  1  + 


I  2  mV 


(30) 
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Bachynski  et  al.  have  derived  the  expression  for  the  contribution  to  the  con¬ 
ductivity  of  a  plasma  by  electron-ion  collisions.  Their  computation  proceeds  from 
the  formula 


4*e2Ne 
ffion  *  ”  Sin 


JfjV3dv _ 

Jfgzd20v2dv 


(31) 


where  the  vector  g  represents  the  flow  velocity  in  velocity  space  (the  collision 
integral  is  related  to  g  through  the  relation  B  *  -  ■  g).  The  vector  g  is  evaluated 

by  considering  the  diffusion  in  velocity  space  and  computing  the  average  of  the  veloc¬ 
ity  change  of  a  teat  particle  after  an  encounter.  The  computation  depends  upon  the 

31 

form  of  the  electron  distribution  function.  Bachynski  et  al.  have  presented  a 
treatment  somewhat  similar  to  that  of  Spltzer  and  Harm, 10  who  have  included  the 
effects  of  electron-electron  collisions  on  the  conductivity.  Both  treatments  are  based 
upon  the  calculation  of  a  transport  coefficient  (conductivity)  in  the  presence  of  a  DC 
electric  field.  The  AC  conductivity  may  be  related  to  the  DC  conductivity  by  the 
approximate  formula: 


ffAC  " 


Ng  e 


m 


Ne  e 


m  a. 


DC 


(32) 


It  may  be  noted  that  Eq.  (32)  is  valid  only  if  i>ion  "  (•'lon>EpF.  where  (vion)EFF 
may  be  found  by  substituting  [  Eq.  (30)]  into  Eq.  (2): 


^ion’EFF 


2  [&ir 

3  </"m 


e4N 


<KTe> 


372 


(33) 


2  2 

which  is  valid  when  w  »  >'EFF' 

Bachynski  et  al.31  have  expressed  the  total  conductivity  of  a  parti  ally -ionized 
gas  (ctotft^)  as  the  geometric  mean  of  two  conductivities  -  one  due  to  electron-ion 
collisions  (<7jon>  and  another  due  to  electron-neutral  collisions  (trn): 


_ L_  .  _L 

‘'total  ‘'ion 


(34) 


This  relation  is  valid,  of  course,  only  for  the  DC  case,  and  is  based  upon  the  assump¬ 
tion  that  the  total  collision  frequency  for  electrons  is  given  by 
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''total  *  ''ion  +  ’'neutrals  ’ 


(35) 


where  »'neutrals  ig  the  electron-neutral  collision  frequency. 

As  noted  by  Shkarofsky  et  al.,32  In  a  later  paper,  a  more  satisfactory  procedure 
than  the  relationship  expressed  by  Eq.  (34)  exists  for  the  computation  of  the  total 
conductivity  of  a  partially-ionized  plasma  taking  into  account  electron-ion  and  electron- 
neutral  collisions.  If  the  expression  for  in  Eq.  (35)  is  used  in  the  formula 

for  the  first  spherical  harmonic  component  of  the  collision  integral  Bj  *  -  ,  f,, 

then  this  expression  for  Bj  may  be  substituted  into  Eq.  (18b)  for  the  first  spherical 
harmonic  component  of  the  distribution  function  fj.  Equation  (18b)  is  then  solved 
for  fj  in  terms  of  v,  and  This  solution  is  then  substituted  into  Eq.  (10), 

so  that  the  conductivity  is  expressed  as  a  ‘•motion  of  tQ,  the  frequency  of  the  impressed 
field  to,  and  otot  The  electron-electron  collisions  produce  two  effects  on  the 
conductivity.  First,  if  vee  »  Gv,  than  fQ  is  Maxwellian  even  in  the  presence  of 
strong  electromagnetic  fields.  Second,  there  is  a  first  spherical  harmonic  compon¬ 
ent  of  the  electron-electron  collision  integral  which  introduces  a  correction  to  the 

total  collision  frequency  appearing  in  Eq.  ( 10)  for  the  conductivity.  This  correction 

17 

appears  in  the  form  (see  Ginzburg,  p.  83): 


■'total-  ‘'total  11  + 


(36) 


where  v  is  given  by  Eq.  (3). 
ee  32 

Shkarofsky  et  al.,  have  computed  the  electron  collision  frequency  of  high- 
temperature  air  as  a  function  of  gas  pressure  and  electron  velocity  by  adding  the 
electron-neutral  collision  frequencies  for  each  constituent  and  then  adding  the  ion- 
electron  collision  frequency  [  Eq.  (30)}  to  obtain  the  total  collision  frequency: 


"total*  "ion Vj^v. 


(37) 


If  w  <  v _ ,  it  should  be  noted  that  the  correction  term  for  interelectron  collisions 

~  ee 

[Eq.  (36)]  is  appreciable,  and  should  be  added  before  substitution  into  the  expres¬ 
sion  for  the  conductivity,  Eq.  (10). 


2.  1.  4  COLLISION  INTEGRAL  FOR  ELECTRON-ELECTRON  COLLISIONS 

Electron-electron  collisions  are  characterized  by  the  long-range  Coulomb  forces 
which  result  in  weak  scattering.  One  of  the  primary  differences  between  electron- 
neutral  and  electron-electron  collisions  consists  in  the  i  atio  between  fractional 
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energy  to  fractional  momentum  loet  per  collielon.  For  electron -neutral  colUslona, 
this  ratio  is  G;  while  for  electron-electron  collisions,  this  ratio  is  one.  For 
electron -elect ran  collisions,  both  the  energy  and  momentum  are  changed  only  slightly 
per  collision,  where  the  fractional  energy  loss  of  an  electron  per  collision  (G##)  is 
given  by  Eq.  (4). 

An  important  consequence  of  the  small  change  in  energy  and  momentum  of  an 
electron  in  electron-electron  encounters  is  the  possibility  of  representing  the  col¬ 
lision  Integral  as  the  divergence  of  a  particle  flux  in  velocity  space  (the  Fokker- 
Planck  expression): 


B  «  -  1 


where  the  particle  flux  density  is  given  by: 

g  ■  \  ^c4?«(c,W  {f(7')  fO?')  -  f<7)  f(V)>  d20d3V  (38) 

c  *  |v  -  V  | 

^ 

Av  *  v'  •  v  . 

Since  the  velocity  of  an  electron  changes  only  slightly  during  an  electron-electron 
encounter,  the  difference  in  the  distribution  functions  before  and  after  the  collision 
may  be  written: 

f(7')  f(V’)  -  f(v)  f(V)  -  [AV  Vy?(V)]  f(v')  -[Av  •  V;  f(v)J  f(V)  ,  (39) 

where  Av  *  v’  -  v  and  AV  ■  V'  -  7.  If  Eq.  (39)  is  substituted  into  Eq.  (38)  and  use 
is  made  of  the  fact  that  e< c,  9)  has  a  maximum  at  9  •  0* ,  then  the  particle  flux  density 
may  be  written: 

g  ■  2 r d3y  •',«(«){«  if<v)  o  ^i<?)  -  fc^>c  • 

+  c2  [  f(v )  Vy}  (7)  -  f(V)  V  1(7)1)  (40) 

♦  j 

where  c  *v 


and 


a« 


vttic) 


2*N, 


e4 


in 


g 

Gins  burg  and  Gurevich  have  shown  how  Eq.  (40)  for  the  particle  flux  density 
may  be  uaed  to  derive  the  zero-order  spherical  harmonic  component  of  the  electron- 
electron  collision  integral: 


C  ■  jt 

*o  +  ¥o)-^]} 


(41) 


where  terms  of  order  f^  have  been  neglected  compared  with  terms  f^.  Here, 
4*  v  (v) 


4*  V  (V)  (*  , 

AjCv)  -  — jjSfi—  \  V2f0(V)  dV 


and 


a2(v)  *  ~T  ^  Vf0(V)dV  +  v3C Vfo(V)  dV]  . 

The  collision  integral  identically  equals  zero  when  tQ  is  Maxwellian.  . 

The  first -order  spherical  harmonic  component  of  the  electron-electron  collision 
Integral  may  be  obtained  directly  from  Eq.  (34): 


where  9  is  the  angle  between  the  relative  velocities  before  and  after  the  collision: 

c'  ■  V'  -  v 1 

?.?-v 
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d 20  *  2w  Bin  9d6 
d2k  =  2»  sin  /9d£  . 

0  is  the  angle  between  v  and  V  (in  the  laboratory  frame).  In  the  center-of-mass 
reference  frame,  the  relative  velocities  of  the  two  particles  (before  and  after  the 
collision)  are  equal  in  magnitude  and  opposite  in  sense:  each  particle  in  the  center 
of  mass  having  a  velocity  equal  to  y/c  /.  As  a  result  of  conservation  of  momentum 
and  energy  /c '  /  *  / c  /  and  \?'  +  v 1  =  V  +  v .  The  velocity  of  the  center  of  mass  is 
equal  to  +  v).  In  the  laboratory  reference  frame,  the  kinetics  of  the  electron- 
electron  collision  may  be  represented  by  the  three  vectors:  c'  *  V'  -  v',  which  is 
the  relative  velocity  of  the  two  electrons  before  the  encounter,  c  »  V  -  v ,  which  is 
their  relative  velocity  after  the  encounter  (equal  in  magnitude  to  their  relative 
velocity  before  the  encounter,  but  rotated  by  the  angle  0),  and  the  vector 

3  .  v*  +  v'  *  V  +  v, 

which  is  the  sum  of  their  velocities,  and  remains  fixed  in  magnitude  and  direction 
before  and  after  the  encounter.  <r(c,  9)  is  given  by  the  Rutherford  formula.  Ginz- 
burg  has  evaluated  the  vector  integral  represented  by  Eq.  (42),  and  demonstrated 
that  the  result  may  be  written: 

5**  v  fj(v)  (43) 

where  v  *  v  +  v.  This  result  is  noteworthy,  for  it  leads  to  a  correction  term 
en  ion 

for  the  total  collision  frequency:  the  sum  of  electron-neutral  and  electron-ion  col¬ 
lision  frequencies  must  be  multiplied  by  the  factor  (1  +  >'ee/‘»)  when  Eq.  (10)  is  used 
in  the  calculation  of  the  conductivity. 


2.  1.  5  COLLISION  INTEGRALS  FOR  IONIZING  AND  RECOMBINATION 
COLLISIONS 

For  most  gas  discharge  phenomena,  such  as  arc  and  glow  discharges,  the 

primary  electron  production  mechanism  is  ionization  of  neutrals  by  electron -impact. 

34 

Bond  has  found  that  electron-impact  of  neutrals  is  the  chief  ionization  process 

behind  shock  waves  in  noble  gases.  On  the  basis  of  ionization  cross  sections  of 
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N2,  Oj,  and  NO  deduced  by  Massey  and  Burhop  from  data  taken  by  Tate  and  Smith, 
Lin  and  Tears4  have  computed  the  specific  ionization  rate  due  to  electron  impact 
on  the  constituents  of  high-temperature  air  as  a  function  of  distance  behind  a  normal 
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one-dimensional  shock.  By  examining  the  relative  importance  of  various  electron 

4 

production  mechanisms  behind  strong  shocks  in  air,  Lin  and  Teare  have  demonstrated 
that  electron-impact  is  a  relatively  inefficient  process  compared  with  atom-atom  and 
atom-molecule  collisions.  However,  if  in  the  presence  of  an  electromagnetic  field 
the  electron-distribution  function  remains  Maxwellian  and  the  electron  temperature 
becomes  equal  to  or  greater  than  10, 000*K  (1  ev),  then  ionization  by  electron-impact 
must  be  considered.  The  dominant  electron  loss  processes  are  dissociative  recom¬ 
bination  (with  positive  molecular  ions)  and  three-body  recombination  (with  positive 
atomic  ions). 

23 

Margenau  has  considered  the  collision  Integral  representation  for  ionization 

and  recombination  collisions.  To  obtain  an  explicit  expression  for  the  ionization 

collision  integral,  it  is  necessary  to  make  an  assumption  regarding  the  division  of 

23 

energy  between  the  two  electrons  after  the  ionizing  impact.  Margenau  assumes 

a 

that  Hie  energy  is  divided  equally  (note  that  Peskoff  takes  into  account  the  case 
where  one  electron  after  the  ionizing  collision  comes  off  with  zero  energy): 

(v1)2  *  2v2  +  Uj 


v'dv'  =  2v  dv 

where  v'  equals  the  velocity  of  the  incident  electron,  ut  *  Ej/m  and  Ej  is  the  ioni¬ 
zation  energy.  Margenau  expresses  the  collision  integral  for  ionizing  collisions  by 
electron  impact  on  a  neutral  constituent  with  a  particle  density  (N^)  in  the  following 
form: 


CVBioniz 


4iv  dv  *  4»N 


nr 

-  *  v 


foaionjv3dv  +  r<)fo<Tionj  v3dv 


(44) 


Here,  (7jon^  is  the  ionization  cross  section  for  the  j**1  neutral  constituent.  The  first 
integral  in  Eq.  (44)  represents  the  production  of  two  electrons  with  velocities  less 
than  v  due  to  an  impacting  electron  with  velocity  greater  than  v.  The  second 
integral  represents  the  appearance  of  one  electron  due  to  ionization  with  an  impact¬ 
ing  electron  having  a  velocity  less  than  v.  When  Eq.  (44)  is  integrated  over  all 
velocities,  the  totel  rate  of  electron  production  due  to  electron  impact  on  the  j**1 
neutral  constituent  is  given  by: 


gioniz  4lrv2dv  .  4)rN  ^ 
0  J  *o 


fo  aionj  y3dv 


(45) 
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Lin  and  Tear*4  have  written  the  expression  tor  the  rate  of  electron  production 

12 

in  term*  of  the  electron  energy  c  «  ^mv  : 

(**) ,  ■  "•  (™) w  “i  ,xpfe)  V’*4* 

where  is  the  threshold  energy  for  the  ionisation  cross  section  of  the  j**1  con¬ 
stituent,  and  fQ  is  assumed  to  be  Maxwellian: 

f  \  3/2  /  2\ 

fo  *  Ne  *xp(na£~J  (47) 

Since  NO  has  the  lowest  ionisation  potential  of  any  constituent  (equal  to  8.  25  ev)  and 
the  maximum  electron  temperature  to  be  considered  will  be  20, 000*K  (2  ev),  kT@  is 
always  much  less  than  the  threshold  energy  for  ionisation  (W^)  so  that  an  asymptotic 
expansion  of  Eq.  (46)  may  be  made: 


WiJ  'ion,  <V 


(48) 


where  cjon^  (Wy)  is  the  slope  of  the  ionisation  cross-section  curve  vs  energy  at  the 
threshold  energy.  By  examining  the  curves  of  ionisation  cross  section  vs  electron 
energy  given  by  Massey  and  Burliop7  for  N^and  Oj,  Lin  and  Teare4  conclude  that,  for 
these  molecules. 


ej  (W.,)W..  «#a*  ■  0.  87  X  10* 16  cm2 
loiij  lj  13  o 

where  aQ  is  the  Bohr  radius. 

Since  the  ionisation  cross  sections  vs  electron  energy  curves  for.  NO,  Nj,  Oj. 
and  the  noble  gases  all  exhibit  the  same  shape  (except  for  a  displacement  in  the 
threshold  energy),  Lin  and  Teare  take 

<onj<Wij,Wi)‘°-87Xl0‘16c“a 
for  the  constituents  N  and  O  also. 

The  total  rate  at  which  electrons  disappear  due  to  recombination  with  positive 
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ion*  of  particle  denaity  Nj  ia  given  by  an  expression  aimllar  to  Eq.  (45): 


where  Nj  *  particle  concentration  of  1th  conatltuent  of  poaitive  ion* 

and  ■  recombination  croaa  aection  for  1th  ion  conatltuent. 

Ri 

The  laotropic  part  of  the  electron  dietrlbutlon  function  ia  aaaumed  to  be  Maxwellian, 
ao  that  the  rate  of  disappearance  of  electron*  due  to  recombination  with  the  1th  con¬ 
atltuent  of  poaitive  Iona  la  expreaaed  aa  a  function  of  electron  denaity  (N^),  poaitive 

ion  denaity  (N.)  and  electron  temperature  CT ).  Now,  the  reaction  rate  for  a  two- 
*  *  + 
body  interaction  auch  aa  dlaaoclative  recombination  (e  +  XY  —  X  +  Y)  la  defined  aa: 


By  examining  Eq.  (49)  and  comparing  it  to  Eq.  (50),  it  may  be  noted  that  the  reaction 
rate  for  dlaaoclative  recombination  la  a  function  only  of  the  electron  temperature. 
TMa  important  reault  la  a  consequence  of  the  fact  that  the  electron  velocity  ia  much 
greater  than  the  velocity  of  the  heavy  ion.  The  moat  general  expreaaion  for  a  reac¬ 
tion  rate,  which  include*  the  caae  where  the  electron  gas  is  not  necessarily  in  ther¬ 
mal  equilibrium  with  the  ion  gas,  may  be  written  in  the  form: 

^  *  k  N#Ni  •  yCf(V,)f(V1)a(Vr)Vrd3Ved3Vi  ,  (51> 

where  V#  ■  electron  velocity 
Vj  « ion  velocity 

Vr  *  ion-electron  relative  velocity. 

Because  the  electron  maaa  is  alwaya  much  less  than  the  ion  mass,  Vr  •  V  .  Hence, 
k  •  Cf(V,>  a(V#)  V#d3V#  . 


(52) 
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4 

Lin  and  Teare  claim  that  the  predominant  electron  production  mechanism  be¬ 
hind  strong  shocks  in  air  is  atom-atom  impact.  Several  of  the  most  important  reac¬ 
tions  include: 


(a)  N  +  O  +  2.  8  ev  -  NO+  +  e 

(b)  N  +  N  +  5.  8  ev  -  Nj  +  e 

(c)  O  +  O  +  6.  9  ev  -  Og  +  e 


(53) 


It  may  seem  surprising  that  reactions  (a),  (b),  and  (c)  are  more  important  for  elec¬ 
tron  production  in  air  shocks  than  ionization  by  electron  impact,  since  electron  im¬ 
pact  was  found  to  be  the  predominant  mechanism  (several  mean  free  paths  behind  the 
shock  front  where  an  appreciable  electron  concentration  has  built  up)  in  Argon  shocks 
(Bond  ).  The  reactions  [  Eq.  (S3))  are  efficient  electron  producers,  even  though 
the  average  kinetic  energy  of  the  atoms  is  considerably  less  than  the  energy  of  the 
endothermic  reactions,  because  the  presence  of  crossing  of  the  potential  energy 
curves  of  the  colliding  atomic  system  lead  to  a  large  cross  section.  However,  direct 
electron  impact  on  neutral  constituents  will  start  to  compete  with  and  even  predomi¬ 
nant  over  neutral-neutral  impact,  as  an  electron  production  mechanism  when  a  high- 
intensity  electromagnetic  wave  is  impressed  upon  the  ionized  flow  field.  From  an 
analysis  of  potential  energy  curve  crossing,  in  conjunction  with  a  somewhat  heuristic 

4 

curve-fitting  procedure,  Lin  and  Teare  arrive  at  an  estimate  for  the  forward  reac¬ 
tion  rate,  (a)  of  Eq.  (53).  The  equilibrium  constant  is  then  determined  from  the 
partition  function.  The  assumption  that  the  detailed  paths  of  the  atomic  states  are 
the  same  in  the  forward  and  backward  directions  (probably  not  a  justified  assumption, 
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see  Biondi  )  permits  Lin  and  Teare  to  write  the  dissociative  recombination  rate 
corresponding  to  the  reverse  reaction  (53a)  as: 


where  kj  is  the  forward  rate  constant 

k  is  the  reverse  rate  constant 
r 

is  the  equilibrium  constant. 

For  Te  <  10*  *K,  Lin  and  Teare  arrive  at  a  dissociative  reaction  rate 
k^  »3X  10'3  Te'3/Z  cm3 /sec  . 


(54) 


(55) 
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(a)  -  7  3 

The  rather  fast  rate  for  dissociative  recombination  predicts  that  k  *  6X10  cm  /sec 

at  T  •  300*,  which  appears  to  agree  with  the  rate  coefficient  of  2  X10  deduced  by 
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Doering  and  Mahan  from  photolysis  experiments  of  NO.  This  piece  of  experi- 

44 

mental  data,  together  with  Sugden's  observation  at  flame  temperatures,  forms 
the  basis  of  Lin  and  Teare's  acceptance  of  Eq.  (55)  as  the  correct  expression  for  the 
dissociative  recombination  coefficient  for  NO+.  Unfortunately  this  particular  choice 
of  temperature  dependence  is  unrealistic  on  physical  grounds.  This  may  be  seen  hy 
writing  Eq.  (52)  in  the  form: 

k  -  4t  (2%) 3/2  £  °(v)  y3  exp(^f] dv 

where  the  electron  velocity  distribution  function  has  been  assumed  Maxwellian.  It 
may  be  noted  that  Eq.  (52)  and  Eq.  (55)  are  compatible  only  if  the  velocity-dependent 
cross  section  assumes  the  form: 


which  is  not  plausible  from  a  physical  point  of  view,  especially  since  the  cross 
section  must  assume  a  finite  value  as  the  electron  velocity  approaches  zero. 

In  the  present  analysis  an  accurate  knowledge  of  the  dissociative  recombination 
rates  corresponding  to  the  reverse  reactions,  (a),  (b),  and  (c)  of  Eq.  (53),  is  neces¬ 
sary  not  only  for  the  computation  of  the  steady -state  electron  density  achieved  in  the 
re-entry  plasma  sheath  under  the  influence  of  high-power  electromagnetic  waves  in 
the  microwave  range,  but  also  the  electron  temperature  must  be  determined  from 
an  energy  balance  equation  by  a  procedure  to  be  outlined  in  Section  2.  3.  2.  The 
terms  in  the  energy  balance  equation  which  represent  energy  loss  of  the  electron 

gas  due  to  dissociative  recombination  Involve  derivatives  of  the  recombination  rate 

-3/2 

coefficients  with  respect  to  electron  temperature.  The  Te  temperature  depend¬ 
ence  for  the  dissociative  recombination  coefficients  of  the  reverse  reactions  (53a, 
b,  c),  as  quoted  by  Lin  and  Teare,4  predict  a  zero  energy  loss  by  the  electron  gas 
on  dissociative  recombination;  a  result  which  is  not  realistic  from  a  physical  point 
of  view.  Since  there  is  a  paucity  of  experimental  data  on  dissociative  recombination 
coefficients  over  extended  temperature  ranges,  the  procedure  to  be  adopted  in  this 
report  will  consist  of  assuming  a  functional  form  for  the  associative  ionization  rate 
[Eq.  (53a)],  and  then  making  a  least-squares  fit  to  the  experimental  data  of  Lin's 
determination  of  the  associative  ionization  cross  section  for  02  -  Ng  mixtures.45 
The  dissociative  recombination  rate  will  then  be  determined  from  Eq.  (54)  (Section 
2.  3.  2). 
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The  recombination  of  electrons  with  the  positive  atomic  ions  0+  and  N+  will 
proceed  primarily  by  the  process  of  three-body  recombinations.  Radiative  recom¬ 
bination  will  be  neglected.  The  concentrations  of  0+  and  N+  become  appreciable  at 
distances  of  ten  or  more  mean  free  paths  behind  the  shock  front  (see  Lin  and  Tears4). 
In  the  absence  of  electromagnetic  field,  the  chief  mechanisms  for  production  of  N+ 
and  0+  ions  are  atom-atom  and  atom-molecule  collisions: 

0+M-M  +  0++e  (56a) 

N+M-M  +  N++e  .  (56b) 

Unfortunately  there  is  at  present  no  reliable  experimental  information  on  the  rates 
of  reactions  of  Eqs.  (56a)  and  (56b).  Even  though  Lin  and  Tears4  Indicate  that  esti¬ 
mates  of  the  three-body  electron-ion  recombination  coefficients  based  on  the  clas- 
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steal  Thomson  theory  yield  rates  that  seem  to  be  several  orders  of  magnitude  too 
IsLrge,  calculations  based  on  the  Thomson  theory  will  be  utilized  in  this  report  tor 
want  of  more  accurate  Information. 

7 

Massey  and  Burhop  have  derived  the  three-body  electron-ion  recombination 
coefficient  by  computing  the  total  probability  that  an  electron  will  suffer  a  collision 
with  a  neutral,  lose  most  of  its  kinetic  energy,  and  simultaneously  remain  within  a 
distance  'ro'  of  a  positive  ion.  ,ro'  is  the  maximum  distance  of  electron-ion  separa¬ 
tion  for  which  the  two  particles  can  describe  a  closed  orbit: 


where  v  >  electron  velocity,  which  is  much  greater  than  the  ion  velocity.  The  Thomson 
cross  section  for  three -body  electron-ion  recombination  is: 


eR(v) 


4»  3  %  1 

'  T  ro  G(v>  Mvj 


(57) 


where  Mv)  •  mean  free  path  for  electron-neutral  collisions.  Using  the  fact  that  the 
electron-neutral  collision  frequency  v(v)  is  given  by  »>(v)  *  v/k(v),  the  rate  of  dis¬ 
appearance  of  electrons  due  to  three-body  electron-ion  recombination  is  given  fay: 


(3*). 


4rv^dv 
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■  Ni  £ $  OjM  ■  4'v2dv  (5B) 

■  NiH.  (337)“  *6  f  1". 

where  G^(v)  is  the  fractional  energy  loss  of  an  electron  in  collisions  with  the  j1*1 
neutral  conatituent 

and  «<j(v)  is  the  electron  collision  frequency  for  collisions  with  the  jttl  neutral  con¬ 
stituent. 
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Eschenroeder  has  calculated  the  three-body  electron-ion  recombination  co¬ 
efficient  according  to  the  Thomson  theory,  but  has  misinterpreted  G  as  the  frac¬ 
tional  energy  loss  of  an  electron  for  electron-ion  collisions,  instead  of  electron- 
neutral  collisions.  For  a  given  neutral  species,  G(v)  may  be  computed  by  solving 
the  integral  Eq.  (26).  The  collision  frequency  v(v)  may  be  computed  for  each  species 

by  using  the  total  cross  section  vs  electron  energy  curves  as  tabulated  by  Shkarofsky 
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et  al.  The  relative  concentration  of  each  neutral  specie  may  be  found  from  Figure 
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4.  1  of  the  report  by  Bachynski  et  al.  It  should  be  noted  that  Eq.  (58)  is  a  generali¬ 
zation  of  the  Thomson  three-body  electron-ion  recombination  formula  for  the  case 
where  the  electron  temperature  is  not  necessarily  equal  to  the  gas  temperature  and 
the  role  of  the  third  body  is  played  by  a  mixture  of  neutral  constituents. 


2. 2  Solutions  of  the  Zero  Order  and  First  Spherical 
Harmonic  Components  of  the  Boltzmann  Equation 


The  zero-order  and  first  spherical  harmonic  components  of  the  collision  Integral 
which  appear  on  the  right-hand  side  of  Eq.  (18a)  and  (18b)  represent  the  sum  of  all 
types  of  electron  collisions: 


B  ■  B®1  +  Blnel  +  Bi0n  +  B*e 

OOO  00 

Bj  »  B®1  +  Bl"el  +  B*°n  +  B®*  . 


Substituting  Eq.  (59a)  into  (18a)  and  (59b)  into  (18b)  yields: 

^*-3  ^  w«rt)  ■  "2[f  £**3} 

*££{v2[A><,M°*A*w*'l} 


(59a) 

(59b) 


(60a) 


+  gioniz  +  Qrecomb. 
o  o 
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8?,  -  3f  v 

™  =-<»»+  >  7,  — f®-  (v  +  »<,  )  r 

at  m  ov  ion  1  uj  ion  1 


(60b) 


where  v  =  S  j"v  cr^  (0,  v)(l-cos0)  d2R 


NT  is  the  particle  density  of  the  jtb  neutral  constituent  and  <r  is  the  cross  section 
for  elastic  and  inelastic  scattering. 

The  term 


v2  to  V2  [_A|(*>  *0  +  A2(v) 


is  of  the  order  (v  f  ).  Hence  the  form  of  the  function  f  depends  upon  the  relative 
ee  o  o 

magnitude  of  the  parameters  vee  and  Gv. 

The  condition  v  »  Gr  is  sufficient  to  guarantee  that  f  be  Maxwellian,  even 
when  ionizing  and  recombination  collisions  are  present  for  the  conditions  considered 
in  this  analysis.  TheBe  conditions  include  the  facts  that  (1)  the  plasma  is  less  than 

0.  1  per  cent  ionized,  and  (2)  the  fractional  energy  loss  of  an  electron  per  collision 

-2 

G  is  always  less  than  or  equal  to  10 

g 

Ginzburg  and  Gurevich  define  a  strongly  ionized  plasma  as  a  plasma  for  which 

the  condition  i»eg  »  Gr  pertains;  and  a  weakly  ionized  plasma  as  one  in  which  the 

condition  v  «  Gv  is  satisfied.  Under  the  assumption  that  B*oniz  and  Biecomb 
ee  r  o  o 

are  zero,  the  Boltzmann  equation  for  the  isotropic  part  of  the  distribution  function 
(Eq.  (60a)]  may  be  solved  for  the  general  case  (v  K  Gv)  by  the  method  of  sue- 

0C  g 

cessive  approximations.  As  noted  by  Ginzburg  and  Gurevich,  this  method  converges 
rapidly  because  variations  in  f  cause  only  small  changes  in  the  parameters  A ^  and 
A2  [defined  by  Eq.  (41)] .  In  the  first  approximation,  a  Maxwellian  distribution 
function 


oo 


N 


3/2 


ej 


(61) 


is  substituted  into  Eq.  (60a).  By  'a  method  t.o  be  described  in  detail  in  Section  2.  3, 
the  electron  particle  density  is  found  by  taking  the  zero-order  velocity  moment  of 
Eq.  (60a)  and  the  electron  temperature  is  obtained  by  taking  the  second  velocity 
moment.  Then,  the  second  approximation  to  Eq.  (60a)  may  be  written: 


f01  »  C  exp 


■  vdv  [z 


ZGjVj  + 


'o  jkT/m  + 


2m 

_M.. 


’ion..- 


2m  „ 

M  io  ty 


»t>1 

_ 

2e2E2r 

3m2(w2+v2) 


+  2 


AS 


(62) 
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where  v  *  < v +  vlon)(l  +  v^Ju)  and  Clu  normalizing  constant.  Here, 

A2  *  IT  A1  1  m  vee<v)  [*W  -  fry  exp<-y2)] 
where 


(error  function) 


In  the  derivation  of  Eq.  (62),  It  hae  been  assumed  that  Bf^/Bt  *  0  (steady-state 
case,  u  »  GjVj).  An  expression  for  f .  has  been  obtained  from  Eq.  (60b),  where 

f.  has  been  assumed  proportional  to  e*  .  This  result  was  first  obtained  in  a  some- 

*  19  8 

what  different  manner  by  Cahn  and  later  by  Ginzburg  and  Gurevich. 

Equation  (62)  is  the  most  general  expression  for  the  electron  velocity  distribu¬ 
tion  function  in  a  partially-ionized  multicomponent  plaBma  (in  which  ionizing  and 
recombination  collisions  are  neglected).  For  positions  on  a  blunt-nosed  re-entry 
vehicle's  surface  at  distances  of  two  or  three  nose  radii  back  from  the  tip,  particu¬ 
larly  for  vehicles  traveling  slower  than  Mach  16  at  200,  000  ft,  the  temperatures 
and  pressures  behind  the  shock  are  such  that  i»ee  “  Gv.  For  this  case,  Eq.  (62) 
should  be  used  in  calculating  the  electrical  transport  properties  (conductivity)  of 
the  plasma.  However,  in  this  report,  only  regions  in  and  near  the  stagnation  region 
of  re-entry  vehicles  at  200,  000  ft  and  traveling  at  about  Mach  16  or  greater  will  be 
considered,  for  which  »  Gv.  In  this  case  it  may  be  shown  that  Eq.  (62)  reduces 

to  Maxwellian  form  (f  ,  *  f _ ).  For  the  other  extreme  [v  _  «  Gv)  it  may  be  shown 

o l  oo  ee 

that  the  general  expression  for  the  isotropic  part  of  the  distribution  function  reduces 
to  the  Margenau  form  [Eq.  (9)J  .  Only  for  the  case  vee  »  Gv  is  the  computation 
of  the  changes  in  electron  temperature,  electron  density,  and  collision  frequency 
due  to  the  presence  of  an  electromagnetic  field  relatively  straightforward,  because 
the  ionization  rates,  recombination  rates  and  fractional  energy  loss  parameters  are 
known  directly  in  terms  of  electron  temperature. 

The  solution  of  the  equation  for  the  first  spherical  harmonic  component  of  the 
electron  velocity  distribution  function  [  Eq.  (60b)]  may  be  written: 


m  |Jw  +  (v  +  »'lon)(l  +  vM/«)J  9v 
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Utilizing  Eqs.  (63)  and  (10),  it  ia  possible  to  find  an  expression  for  the  AC  con¬ 
ductivity  of  a  partially-ionized  plasma  in  terms  of  the  isotropic  part  of  the  electron 
distribution  function  (assumed  to  be  Maxwellian  in  this  case),  the  frequency  of  the 
applied  field  (o>),  the  total  electron-neutral,  electron-ion,  and  interelectron  collision 
frequencies,  and  the  electron  density: 


a 


vfd3v 


‘f  J  v^-  7)  2*  sin  0d0vdv  *  i.  (x)  Pi 


v3dv 


(64) 


where  f }  lies  along  the  E  field,  taken  along  the  Z-axis.  Substituting  Eq.  (63)  into 
Eq.  (64)  yields: 


8e2N„  I"p«  v(u)(l  +  u4e‘u2du  p*°  u4e'uZdu 

\  o  o  o  " 22 


0  =  3^m 


-i-0  a,2  +  i/2  (1  +  vee/u>)2 


o  w2  +  v2  (1  +  v  /w)^J 
ee 


(65) 


where 


v(v)  *  2  N .  ^  v  +  v,on 


/  m 

u  s  V  2kTe  v 

4  'll  m2v4N 
v_(v)  =  2x  N_  -—T  In,  - 


ee 


e 

m  v 


The  conductivity  is  related  to  the  dielectric  constant  (MKS)  through  the  relation: 


K  =  1  + 


JW£o  ' 


(66) 


Typical  parametern  for  regions  in  and  near  the  stagnation  region  of  blunt-nosed 
re-entry  vehicles  at  200,000  ft  at  Mach  16  are: 

Ne  m  5  X  1012  elec/sec 

u  w  io 1 1  rad/sec 

P 


v  *  5  X  10s  coll/sec 


v#e  »  10H  coll/sec 


G  »  5  X  10 


»  50Gy 


For  the  caae  of  a  monochromatic  plane  wave  incident  on  a  slab  of  plasma  (having 

the  characteristics  of  the  re-entry  sheath),  most  of  the  electromagnetic  energy  will 
be  coupled  into  the  plasma  only  if  the  frequency  of  the  plane  wave  (w)  lies  within  a 
factor  of  two  or  three  times  greater  than  or  less  than  the  plasma  frequency 
u_/3  <  w  <  3u  .  Otherwise,  most  (BO  per  cent  or  more)  of  the  energy  will  be  either- 
reflected  from  the  first  interface  or  transmitted  through  the  slab.  However,  since 
the  nonlinear  properties  of  the  plasma  sheath  (changes  in  electron  density  and  col¬ 
lision  frequency)  will  be  manifested  only  when  there  is  appreciable  coupling  between 
the  electromagnetic  wave  and  the  plasma,  the  following  conditions  will  pertain  to  the 
range  of  parameters  where  maximum  nonlinearity  occurs: 


2^2 
O)  »  v 


and  v _ /w  «  1  . 

ee 


This  Implies  that  the  dielectric  constant  may  be  written  in  the  form: 


where  vEFF  is  defined  by  Eq.  (2): 


kEFF  WT 


and 
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2  2 

Hence,  for  the  high-frequency  case  (u>  »  v  )  the  total  effective  collision  frequency 
is  the  sum  of  the  effective  collision  frequency  between  electrons  and  neutrals  and  the 
effective  collision  frequency  between  electrons  and  ions: 

"EFF  =  ‘•'en^pp  +  (,WEpp 


where  (v.  )  is  given  by  Eq. 

*on  EFF 


(33): 


<pion) 


EFF 


-  2  /8» 
”  3  V  m 


N?e2 

(kT  )& 
e 


In 


/kvD) 

V  e2  J 


Thus  Eq.  (67)  for  the  dielectric  constant  will  form  the  basis  for  the  computations 

of  the  nonlinear  transmission  characteristics  of  the  re-entry  plasma  sheath.  The 

parameter  (v  )_,__  may  be  calculated  from  a  knowledge  of  the  electron  temperature 
en  Lr 

Te<  electron  density  Nfi,  and  the  concentrations  of  the  neutral  constituents  Ny  to¬ 
gether  with  their  respective  cross  sections  The  concentrations  of  the  neutral 
constituents  are  calculated  from  the  graphs  of  the  total  relative  density  (p/po)  pro¬ 
file  and  gas  temperature  (T)  profile  behind  the  shock  front  found  in  the  AFCRL 
report  by  Rotman  and  Meltz 1 1  and  the  curves  of  the  relative  concentrations  of  the 
neutral  constituents  O,,  N,,  O,  N,  and  NO  plotted  as  a  function  of  p  Ip  and  T 

located  in  the  report  by  Bachynski  et  al.  (Figure  4. 1).  The  total  cross  sections 
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for  electron  scattering  <^(v)  may  be  found  in  the  report  by  Shkarofsky  et  al.  for 

the  various  constituents  O,.  N_,  N,  O,  and  NO  plotted  as  a  function  of  electron 
12  "  * 

energy  (^  mv  ).  The  procedure  for  the  determination  of  Te  and  Ng  at  each  point 
in  the  plasma  slab  will  be  discussed  in  Section  2.  3. 


2. 3  Particle  Conservation  and  Energy  Conservation: 

The  Zero  and  Second  Velocity  Moments  of  the  Boltzmann  Equation 


2.  3.  1  THE  STEADY -STATE  ELECTRON  CONCENTRATION 

When  v _ »  Gp,  f  is  Maxwellian  and  an  expression  for  the  electron  density 

wv  o  2 

may  be  obtained  by  multiplying  both  sides  of  Eq.  (60a)  by  4wv  dv  and  integrating. 

The  following  expressions  are  the  only  nonvanishing  terms  which  result  from  the 

Integration: 

8  C  2  8Ne 
■W  }  fo  4,v  dv  '  ~W 


(68a) 
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(68b) 


■  ®v  £  [?  °)<v)  -i<v07 

Equation*  (68b),  (88c),  and  (68d)  have  been  derived  in  Section  2.  1.  5.  Since,  in  the 
steady  state,  the  electron  density  is  equal  to  the  sum  of  the  ion  densities  (N0  *  LN^), 
the  zero-order  moment  of  the  Boltzmann  equation  for  the  isotropic  part  of  the  dis¬ 
tribution  function  may  be  written: 


It  should  be  noted  that  Eq.  (88)  does  not  include  the  effects  of  electron  production 
resulting  from  'external'  agencies,  such  as  photoionization,  and  neutral-neutral 
impact.  >, 

On  the  basis  of  Lin  and  Teare's  Investigation,  the  examination  of  the  many 
possible  electron  production  and  recombination  mechanisms  (such  as  atom-atom 
and  atom-molecule  impact,  electron-neutral  and  electron-ion  impact,  photoioniza¬ 
tion,  electron-neutral  attachment,  dissociative  recombination,  and  three-body  recom- 
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bination)  indicates  that  the  following  reactions  are  the  predominate  processes  in  the 
presence  of  an  electromagnetic  field  which  raises  the  electron  temperature  from  O.S 
to  2  ev: 


N  +  0  -  NO+  +  e 

(1) 

N  +  N  -  N  2  +  e 

(2) 

O  +  0  -  Oj  +  e 

(3) 

X  +  N-N++X+e 

(4) 

X  +  0-0++X+e 

(5) 

NO  +  e  -  NO+  +  2e 

(6) 

02  +  e  -  02  +  2e 

(7) 

N2  +  e  -  N2  +  2e 

(8) 

N  +  e  -  N+  +  2e 

(9) 

O  +  e  -  0+  +  2.  e 

(10) 

Ne  »  N+  +  0+  +  N2  +  02  +  NO+  . 

(11) 

Then,  the  rate  of  production  for  each  ionic  specie  is  given  by  the  expressions: 


-  k,  (N)(0)  +  k,  (NO)(N  )  -  k  (NO+)(N  ) 
dt  fj  f6  e  r,  e 


ml)  2  + 

“df  *  kf2  <">  +  \  *  kr2  <N2)<Ne) 


d<02) 


kf<0)3  +  kf,<02)(Ne>-kr  .<°>e> 


i^-k,  (X)(0)  +  k.  (0)(N  )  -  k  (0+)(X)(N  ) 
m  ^5  '10  *  r5  * 

.  kj  (N)(X)  +  kf  <N)(Ne>  -  kp  (N+)00(Ne)  . 
4  9  4 


(71) 


Using  particle  conservation,  charge  neutrality,  and  the  steady-atate  condition 
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(^■5sni'°)  ' 

an  equation  for  the  steady-state  electron  density  may  be  obtained  from  expressions 
(71): 


where  is  the  particle  concentration  of  the  jtl1  neutral  constituent.  The  relative 

concentrations  of  N.  O,  N,,  O,.  and  NO  as  a  function  of  plpn  and  gas  temperature 
«  •  38  ° 

T  are  taken  from  Figure  4.  1  of  Bachynski  et  al.  The  ionization  potentials  of  these 
constituents  are  given  in  Table  2: 

TABLE  2 


The  relative  density  (p/pQ)  and  temperature  profiles  along  the  line  located  at 
45*  with  respect  to  the  axis  of  the  hemisphere  cylinder  are  given  in  Rotman  and 
Mettz. 1 1  These  profiles  were  determined  for  equilibrium  flow  conditions  in  which 
the  boundary  layer  is  neglected.  The  nonlinear  electromagnetic  transmission 
characteristics  of  the  re-entry  sheath  computed  in  this  report  will  be  based  upon 
the  calculations  presented  in  the  report  by  Rotman  and  Mettz11  for  the  line  located 
at  45*  with  respect  to  the  axis  of  the  hemisphere  cylinder. 
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The  forward  rate  for  reaction  (1)  will  be  obtained  by  taking  a  theoretically  de- 

rived  expresaion  for  the  associative  ionization  cross  section  and  making  a  least- 

45 

squares  fit  of  this  functional  form  to  the  experimental  data  of  Lin.  This  procedure 
will  be  explained  in  detail  in  Section  2.  3.  2  in  connection  with  the  energy  loss  of  the 
electron  gas  in  dissociative  recombination.  This  will  give  the  expression  for  kf  . 
The  expression  for  the  dissociative  recombination  coefficient  kr^  will  be  obtained 
through  the  equilibrium  constant 


where  K^*3  «  (1.  4  X  10‘8T  +  1.  2  X  10‘12T2  +  1.  4  X  10‘16T3)  exp  (-32.  500/T). 

For  want  of  better  experimental  data,  it  will  be  assumed  that  the  dissociative  recom¬ 
bination  coefficients  of  the  reverse  reactions  1,2,  and  3  will  all  be  equal: 


k  »  k  *  k 
rl  r2  r3 


Then  the  forward-rate  coefficients  for  the  associative  ionization  reactions  (2)  and  (3) 
will  be  determined  from  their  respective  equilibrium  constants: 


and 


k  «<Kf3>k 
f2  2  rl 


k,  *  (K?3)  k 
f3  3  rl 


where  the  equilibrium  constants  for  reactions  (2)  and  (3)  have  been  obtained  from 
the  partition  functions  (Lin  and  Teare4): 

K^3  -  (3  X  10"®T  +  4  X  10’  12T2  +  10'  l5T3  -  3  X  10*2aT4>  •  exp(-67,  300/T) 
K®3  •  (1.6X  10‘®T  +  1.  2X  10'1ZT2+  3.  5X  10*18T3)  •  exp(-80,  100/T)  . 


The  forward  rates  for  reactions  (4)  and  (5)  were  estimated  by  the  simple  col- 
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lision  theory  presented  by  Bortner  and  tabulated  in  the  report  by  Atallah: 
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<n!  k") 
J  rj 


(1) 

(N)(0)  kf 
*1 

(2) 

(N)(N)  K?Qkr 

2  rl 

(3) 

(OXO)  K^k 

3  rl 

(4) 

(X)(N)  4.  86  X  10‘12T^exp<-333,  000/RT) 

(5) 

(XXO)  4.  62  X  10‘12T^exp(-313,  000/RT) . 

Lin  and  Teare 

*  have  pointed  out  that  there  is  no  direct  experimental  data  available 

in  the  literature  for  the  forward  rate*  (4)  and  (5).  The  backward  rates  may  be 

computed  on  the  baaia  of  the  Tbomaon  theory.  It  ahould  be  noted  that  the  terma 

<NJk£  ,  which  appear  in  expreaaions  (67)  for  the  electron  denalty  (N  ),  all  play  the 
J  1  ® 

role  of  a  constant.  Theae  terms  are  constant  for  they  are  functions  only  of  the 

neutral  concentrations  and  the  gas  temperature  (and  not  electron  temperature). 

Both  the  gas  temperature  and  the  concentrations  of  the  neutral  constituents  are  as¬ 
sumed  to  be  unchanged  in  the  presence  of  the  electromagnetic  field.  This  is  a  good 
assumption  because  the  degree  of  ionisation  of  the  plasma  sheath  ia  between  0.  1  and 
0.  01  per. cent  or  less  (p/pQ  *  10"3  to  10"2,  N#  »  1012  or  1013  elec/cm3). 

The  terms  in  Eq.  (72)  refer  to  the  backward  rates  of  reactions  1  through  5. 

The  backward  rates  1  through  3  refer  to  dissociative  recombination.  Since  the 

-3/2 

temperature  dependence  (T#  )  of  the  dissociative  recombination  rates  correspond¬ 

ing  to  the  backward  rates  1  through  3  as  given  by  Lin  and  Teare  Is  inadequate  for 
the  purposes  of  the  present  analysis  because  it  predicts  zero-energy  loss  of  the 
electron  gas  on  dissociative  recombination,  the  dissociative  recombination  rate  for 
the  reverse  of  reaction  1  will  be  determined  by  making  a  least-squares  computation 
to  be  explained  in  Section  2.  3.  2.  The  dissociative  recombination  coefficients  cor¬ 
responding  to  the  reverse  reactions  1  through  3  will  all  be  assumed  equal.  The 
backward  rates  of  reactions  4  and  5  refer  to  the  Thomson  three-body  recombination 
coefficient,  generalised  to  the  case  of  a  multicomponent  plasma  where  the  electron 
temperature  is  not  equal  to  the  gas  temperature  [  Eq.  (58)  J : 


Gj(v)  is  the  fractional  energy  loss  of  an  electron  per  collision  with  the  neutral 
constituent.  The  parameter  G  is  usually  determined  as  a  function  of  electron 
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temperature  (T^),  and  la  liated  In  Table  1  for  varloua  gaaea.  Notice  that  the  three- 
body  lon-electron  recombination  ratea  [  Eq.  (58)]  do  QOJ,  have  to  be  calculated  fay 
flrat  finding  G(v).  Thia  would  involve  aolving  the  integral  Eq.  (26)  for  G(v)  for 
each  constituent,  and  then  aubstitutlng  the  reault  into  Eq.  (58).  A  much  more  ele¬ 
gant  technique  consists  in  taking  the  right-hand  side  of  Eq.  (26)  fimown  for  each 
constituent)  and  integrating  four  tlmea  with  reapect  to  electron  temperature.  Thia 
will  yield  the  correct  form  for  the  generalised  Thomson  three-body  recombination 
rate: 


C"  -^5  A  C  W^2G,(w)r.  (wWwdvdude. 

*  o  u*  *  o  v*  *  o  J  JEFF 


The  effect  of  charge  exchange  has  been  neglected  in  the  derivation  of  Eq.  (72)  for  the 
electron  density.  This  mechanism  appears  to  have  a  relatively  small  effect  on  the 
total  electron  concentration.  Under  the  action  of  the  charge  exchange  mechanism, 
the  relative  concentrationa  of  the  various  ionic  species  may  change.  This  may 
produce  a  change  in  the  total  electron  concentration,  because  the  electron-ion  re¬ 
combination  rates  are  specie  dependent. 


2.  3.  2  THE  ENERGY  BALANCE  EQUATION 

When  the  relaxation  time  for  energy  transfer  between  the  electron  gas  and  the 
neutral  gas  is  much  greater  than  the  period  of  a  monochromatic  electromagnetic 
wave  impreased  upon  the  plasma  (t^  *  ^  »  T  ■  2t /a>),  then  the  electron  gas  will 
no  longer  be  in  thermodynamic  equilibrium  with  the  neutral  gas  but  will  achieve  a 
temperature  that  depends  upon  the  mean  square  value  of  the  impressed  field.  In 
the  steady  state  the  temperature  of  the  electron  gas  ia  determined  from  an  energy 
balance  equation.  This  energy  balance  equation,  obtained  by  taking  the  second 
velocity  moment  of  the  Boltzmann  equation  for  the  Isotropic  part  of  the  distribution 
function  ( Eq.  (80a)] ,  is  an  expression  of  the  fact  that  in  the  steady  state  the  energy 
acquired  by  the  electron  gas  due  to  ohmic  heating  by  the  electromagnetic  field  is 
equal  to  the  various  electron  energy  loss  processes.  The  electron  energy  loss 
processes  for  a  multicomponent  plasma  Include: 

<1)  Elastic  electron-neutral  collisions 

(2)  Elastic  electron-ion  collisions 

(3)  Electron-molecule  collisions  which  excite  rotational,  vibrational,  and 
electronic  levels 

(4)  Ionizing  collisions  by  electron  Impact  on  neutrals 

(5)  Ion-electron  dissociative  and  three-body  recombination  collisions. 
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All  of  these  processes  will  be  taken  into  account  in  the  computation  of  the  steady- 
state  electron  temperature.  Other  electron  energy  loss  processes  for  a  particular 
point  in  the  plane-layered  plasma  slab  that  may  be  neglected  in  the  energy  balance 
equation  Include: 

(6)  Heat  flow  due  to  electron-electron  collisions 

I 

(7)  Heat  flow  due  to  electron-neutral  collisions 
f  KenfTe)  V\  -  KenfT)  V2T1 

(8)  Heat  flow  due  to  particle  diffusion 

VN  •  VT  -  SA& - VN  VT  . 

J^EFF  e  e  fl'm,,EFF  e  J 

These  last  three  electron  energy  loss  processes  will  be  shown  to  be  negligible  in 
Section  2.  4. 

There  are  several  important  facts  that  bear  careful  consideration.  First,  the 
term  (2m/M)vion  in  Eq.  (60a)  refers  to  elastic  collisions  between  electrons  and  a 
single  ion  specie  of  'average'  mass  M.  The  general  term  for  a  plasma  containing 
many  ionic  species  would  be  written: 


where  Mt  is  the  mass  and  Nj  the  particle  concentration  of  the  ith  ionic  specie.  The 
approximation  consists  of  replacing  by  an  average  mass  M  and  XNj  *  Ng  (steady 
state).  This  results  in  a  great  simplification  of  the  equations  while  introducing  only 
a  small  error. 

Another  important  point  involves  the  precise  balancing  of  all  energy  gain  and 
loss  processes  of  tne  electron  gas.  It  was  previously  mentioned  that  electron 
production  mechanisms  due  to  external  agencies,  such  as  neutral-neutral  impact, 
are  not  included  in  the  collision  integral  terms  appearing  on  the  right-hand  side  of 
the  Boltzmann  equation  [Eq.  (80a)].  Processes  such  as  associative  ionization  cor¬ 
responding  to  the  forward  reactions  (1)  through  (3)  and  ionization  by  neutral-atom 
impact  corresponding  to  the  forward  reactions  (4)  and  (5)  have  been  taken  into  ac¬ 
count  in  the  expression  for  the  steady-state  electron  concentration,  Eq.  (72). 
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Moreover  the  forward  reaction*  (1)  through  (5)  constitute  a  net  energy  gain  fay  the 
electron  gas  which  is  not  changed  by  the  presence  of  an  electromagnetic  field  since 
the  forward  reactions  do  not  depend  upon  the  electron  temperature.  However  the 
dissociative  and  three-body  recombination  processes  corresponding  to  the  backward 
reactions  (1)  through  (S)  constitute  an  energy  loss  to  the  electron  gas  which  does 
depend  upon  the  electron  temperature. 

The  consideration  of  the  energy  gain  and  loss  terms  due  to  ionizing  collisions 
by  electron  impact  corresponding  to  the  forward  reactions  (6)  through  (10)  is  more 
subtle.  The  creation  of  a  new  electron  by  electron  impact  on  a  neutral  constituent 
represents  a  gain  of  energy  by  the  electron  gas,  which  may  be  computed  by  taking 
the  second  velocity  moment  of  the  collision  integral  for  ionizing  collisions.  The 
impacting  electron,  which  is  responsible  for  the  ionization,  suffers  a  loss  of  energy, 
which  is  included  in  the  parameter  G.  This  is  because  G  represents  the  total 
fractional  energy  loss  of  an  electron  due  to  all  electron-neutral  collision  processes- 
elastic  electron-neutral  collisions,  collisions  which  excite  vibrational,  rotational, 
and  electronic  levels  of  molecules,  and  ionizing  collisions.  In  the  absence  of  an 
electromagnetic  field,  when  the  electron  gas  is  in  thermodynamic  equilibrium  with 
the  neutral  gas,  each  of  these  numerous  processes  occurs  at  a  rate  such  that  there 
is  no  net  gain  or  loss  of  energy  by  the  electron  gas.  The  application  of  an  electro¬ 
magnetic  field  to  the  ionized  flow  field  results  in  an  elevation  of  the  electron  tem¬ 
perature  with  a  consequent  change  in  some  of  the  rate  processes.  The  energy  lost 
t»y  the  electromagnetic  field  to  the  plasma  is  dissipated  only  by  those  energy  loss 
processes  over  and  above  what  is  suffered  by  the  electron  gas  when  in  equilibrium 
with  the  neutral  gas.  This  implies  that  each  energy  gain  or  loss  process  of  the 
electron  gas  must  be  represented  in  the  energy  balance  equation  as  the  difference 
between  two  terms.  One  term  refers  to  the  energy  gained  or  lost  by  the  electron 
gas  in  the  presence  of  the  electromagnetic  field  at  the  elevated  electron  temperature 
(T  )  and  electron  density  N  (T  ),  whereas  the  second  term  refers  to  the  energy 
gained  or  lost  by  the  electron  gas  when  in  equilibrium  with  the  neutral  gas  at  tem¬ 
perature  <T)  and  electron  density  NeC T).  Thus,  it  may  be  seen  that  each  term  in 
the  energy  balance  equation  refers  to  an  energy  gain  or  loss  by  the  electron  gas 
corresponding  to  the  difference  between  two  states  of  the  electron  gas:  one  state  of 
the  electron  gas  is  at  the  temperature  T#  and  electron  density  Ne(Te>  in  the  presence 
of  the  electromagnetic  field,  the  other  corresponds  to  the  equilibrium  state  at  tem¬ 
perature  CT)  and  electron  density  NefT).  The  appropriate  terms  in  the  energy 
balance  equation  are: 

The  time  rate  of  change  of  the  total  kinetic  energy  in  the  electron  gas: 


(74a) 
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The  ohmic  heating  of  the  electron  gas  by  the  electromagnetic  field: 

ip  2  •  f  v3ridv  .  .  E2  . 

3  *0  1  r  Bltf 


(74b) 


where 

v  »  (v  )  + 

EFF  <lonEFF  yJT 


.CSy4 «.(*$*■ 


'Electron  energy  lose  due  to  elastic  electron-ion  collisions  and  Inelastic 
electron-neutral  collisions. 

2#m  C  If  >'ion  +  *°J<V>,'J<V3  y3  ^  8^"  +  *'fo)  dv  " 


(74c) 


!  W  N(Te>  WVT>(,Weff  +  I  SGj<V  [l'EFFfre)]jkNefre)fre-T) 


where 


C  ,  .2/1?  W?  J^\ 

*on  EFF  ^  ** m  (Jr  )®^  v  e2  J 

Energy  gain  of  the  electron  gas  due  to  creation  of  new  electrons  by  ionizing 
electron-neutral  collisions, 

£  (2*m)  £v4  |s*oni*  (Te)  -  (T)]  dv 

.  2wm  l  M,  f  [W  -  fGf T)\  aionj  v5dv 

Z  (2»m)  (*” v4  jBp0ni*(Te)  -  Bp0nl*rr)]  dv  *  §  k  £  rjv,o°nU(Te)'T  Bo0”1*0’] 

•  4w2dv  +  kT2  Z  C Bp0ni*(Te)  4»v2dv-kT2  jjt  £  r*Bo°"lla’)4w2dv 
•  *  0  1  0 

* B°°niIfr«>  ^  +  ot2  ©it  2  £ B°0nUfr) 


(74d) 


2 

•4fv  dv 
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/8kT  \  ^  f 

*  ZNj  «p(- Kfj)  <«J 

/  g^p  N  ^  /"^j\ 

+  W^FStJ  ZNJWij  exp^r-j  («J) 

•«™.<n(S)WZ»,«p(T^  <-4 

*Nefr)(^)1/2  ZNjWij  exp^)  («J)  . 

Energy  loss  of  the  electron  gas  due  to  dissociative  recombination 


2irm  l  r  fel8,0C-  reCOmb  ff  )  -  B*bboc-  recomb  (T>|  v4dv 
*  o 

*  2«n  ^  Nt  HW  -  f  (T)]  cDR  V5dv 

4  O  1 

■  2  kTe  W  l  W  W  +  ^^e'  l  Wlr^  <74«> 

,  9  dk.CT) 

-  |  kTNe(T)  2]NjfT)  k.(T)  -  kT‘  NefT)  ^NjCr)  —*£-  . 
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As  mentioned  previously,  Lin  and  Teare's  choice  of  fTe  )  for  the  tempera¬ 
ture  dependence  of  the  dissociative  recombination  rates  (corresponding  to  the  reverse 
reactions  1,  2,  and  3)  is  unrealistic  from  a  physical  point  of  view  because  it  requires 

the  velocity-dependent  cross  section  for  dissociative  recombination  to  assume  the 

3  -3/2 

form  c(v)  *  8(v)/v  ,  where  fl(v)  is  the  Dirac  delta  function.  A  (Te  )  temperature 

dependence  of  the  rate  coefficient  predicts  zero  energy  loss  by  the  electron  gas  for 
dissociative  recombination.  This  Implausible  result  is  a  consequence  of  the  curve - 
fitting  procedure  Lin  and  Teare  have  used  for  the  determination  of  the  associative 
ionization  rate  corresponding  to  the  forward  reaction  1.  Lin  and  Teare  have  ex¬ 
pressed  the  general  functioned  form  of  the  associative  ionization  cross  section  on 
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the  basis  of  the  Bates  and  Massey  potential  energy  curve-crossing  model: 


Qe>ex  '  «’<Vo)2 


so 


where  Ex  equal*  activation  energy  and  g  is  the  ratio  between  the  statistical  weight 
associated  with  the  initial  potential  energy  curve  which  leads  to  the  crossing  point 
and  the  sum  of  the  statistical  weights  associated  with  all  possible  initial  potential 
energy  curves  of  the  colliding  atomic  system,  p  is  the  reduced  mass  of  the  col¬ 
liding  system  and  R  a  is  the  distance  between  the  two  centers  of  mass  of  the  col- 
*  ®  12 
liding  system  at  the  point  of  crossing.  E  is  the  impact  energy  E  *  ^  pv  (in  electron 

volts)  and  is  the  autoionization  lifetime  of  the  molecular  complex.  While  there 

is  no  definite  Information  on  the  value  of  the  autoionization  lifetime  of  complex 

- 14 

atomic  systems,  the  assumption  that  t^j  «  10  leads  to  a  velocity  dependent  cross 
section  expressible  in  the  form: 

QE  >  Ex  *  AE  ^ 

.where  A  is  a  constant.  This  leads  to  a  velocity -averaged  cross  section  of  the  form: 
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Lin  has  measured  the  rate  for  the  associative  ionization  reaction  1  in  the  tem¬ 
perature  range  4000*  to  5000TC  by  monitoring  both  the  DC  conductivity  and  the  micro- 
wave  attenuation  in  shock  heated  0.  25%  C>2  -  99.  75%  N2  mixtures.  The  velocity- 
averaged  associative  ionization  cross  sections  for  this  data  are  plotted  as  a  function 

4 

of  temperature  in  Figure  5  of  the  report  by  Lin  and  Teare.  In  their  report  the 
velocity -averaged  associative  ionization  cross  section  for  the  extended  temperature 
range  300  <T  <  30,  000*K  is  determined  by  fitting  the  functional  form  Q  *  AT*1* 
exp(-Ex/kT)  to  the  center  of  the  experimental  points  in  Figure  5  and  choosing  the 
activation  energy  equal  to  its  lowest  possible  value,  the  heat  of  the  reaction  E%  * 

2.  8  ev.  When  the  dissociative  recombination  coefficient  corresponding  to  the 

reverse  reaction  1  is  obtained  by  dividing  the  forward  rate  by  the  equilibrium  con- 

-3/2 

start,  this  procedure  leads  to  a  recombination  rate  with  a  T0  temperature  de¬ 
pendence.  Since  this  approach  leads  to  the  physically  unrealistic  situation  where 
there  is  no  electron  energy  Iosb  on  dissociative  recombination,  a  different  proce¬ 
dure  will  be  adopted  in  the  present  analysis.  A  least-squares  fit  will  be  made  of 
the  functional  form  for  the  velocity-averaged  associative  ionization  cross-section 
Q  *  AT"  ^  exp(-B/T)  corresponding  to  tAj  »  10" 14  sec  to  the  experimental  points 
given  in  Figure  5  of  Lin  and  Teare.  Since  the  activation  energy  must  be  equal  to 
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or  greater  than  the  heat  of  the  reaction,  the  restriction  B  >  32,  500  will  be  added. 
A  least-squares  fit  may  be  accomplished  by  minimizing  the  function 


W  *  ^  {  Q(A,  B.T. )  -QJ  2 
k«l  K  K 

where  Q(A.  B,  Tfc)  *  ATk’^  exp(-B/Tk>  and  Qk  (a  constantlcorreaponds  to  the  kth 
experimental  point  at  the  temperature  Tk  (Figure  5).  The  minimization  is  achieved 
from  the  relations  8W/8A  *  0  and  8W/8B  *  0,  which  lead  to  the  expressions: 

Jj  [Q(A,  B.T^)  -  Qk)  Q(A,B,Tk)=0 


and 


J  1<5<A.  B.T.  )  -  Q.  1  Q(A,  B.T.)  *  0. 

k»l  K  K  * 


It  is  relatively  easy  to  eliminate  the  parameter  A  from  these  two  equations  so  that 
a  single  equation  remains  to  determine  the  parameter  B.  Of  course  the  result  of 
the  calculation  must  yield  a  value  for  B  such  that  B  >  32,  500.  Once  the  parameters 
A  and  B  are  determined  from  the  least -squares  fit,  the  associative  ionization  rate 
is  determined  from  the  formula 


kfj  Q(A,B,T)  AT_1/2exp(-|)  . 

The  dissociative  recombination  rate  is  then  obtained  from  the  equilibrium  constant 
J».  SSj  A  exp  (-  ® ) 


krirre> 


_ii_ 

k?Q- 


where 


K®6*-  -  (1.4X  10'^  +  1.2  X  10"12T2  +  1,4  X  10'IftT3)  exp(-32.  500/T). 


It  should  be  mentioned  that  this  procedure,  while  leading  to  a  more  realistic  func¬ 
tional  form  for  the  dissociative  recombination  rate,  suffers  from  the  same  weakness 


52 
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as  the  development  of  Lin  and  Teare.  Bates  has  specifically  pointed  out  that 
attempts  to  use  the  principle  of  detailed  balancing  to  determine  the  cross  section 
for  associative  ionisation  from  measurements  o*  the  dissociative  recombination 
cross  section  are  not  justified.  This  is  because  the  systems  in  their  initial  states 
are  usually  in  their  ground  state,  whereas  the  final  system  is  expected  to  be  in  an 
excited  state.  Nevertheless  the  principle  of  detailed  balancing,  though  not  strictly 
applicable,  will  be  used  in  this  report  to  obtain  the  dissociative  recombination  rate. 
Finally  it  will  be  assumed  that  the  dissociative  recombination  coefficients  of  the 
reverse  reactions  1,  2,  and  3  will  all  be  'equal.  The  values  thus  obtained  for  these 
rate  coefficients  will  be  utilised  until  more  reliable  experimental  information  is 
available. 

Energy  loss  of  the  electron  gas  due  to  three -body  recombination: 

i. 

2*m  £  "  [gl*,ree"hody  recomb.  .  gthree-body  recomb^  y4dv 
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Since  only  steady-state  conditions  are  to  be  considered,  the  term  d/dt  (3/2  kTeN#) 
will  equal  sero  and  the  energy  balance  equation  becomes: 
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A  summary  of  symbols  used  In  the  energy  balance  Eq.  (75)  follows:  N  ( T  )  Is  the 

®  •  2 

electron  density  as  a  function  of  electron  temperature  given  by  Eq,  (72).  E  is  the 
mean  Square  value  of  the  total  field  at  a  particular  point  in  the  plasma. 
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N.  refers  to  the  concentration  of  the  j  neutral  constituent.  W, .  refers  to  the  ioni- 

J  V 

zation  potential  of  the  j  neutral  constituent  (Table  2) 

»a2  *  0.876  X  10*l6cm2 
o 


M  is  the  averaged  ionic  mass  =  24  amu 

GjCiy  is  the  fractional  energy  loss  of  an  electron  in  collision  with  the  jth  neutral 
constituent  (Table  1) 


t  ‘,EFF(Te)1j  °  3<f  (kT  )  ^  Nj<Tj(v>V 
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N.  refers  to  the  concentration  of  the  i4*1  molecular  ionic  constituent  (NO+,  O^.  N*) 


as  a  function  of  electron  temperature. 
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refer*  to  the  concentration  of  the  k**1  atomic  ionic  constituent  (0+,  N+)  a*  a 
function  of  electron  temperature. 
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For  a  given  value  of  the  mean-square  electromagnetic  field  (E  ),  the  electron 
temperature  at  a  particular  point  in  the  plasma  may  be  determined  by  solving  the 
energy  balance  equation  [  Eq.  (75))  for  (Te).  All  of  the  parameters  which  appear  in 
the  energy  balance  equation  are  capable  of  being  expressed  as  functions  of  electron 
temperature  in  analytic  form,  except  the  terms  (*'en)EFF-  GjfTg).  and  the  energy 
loss  of  the  electron  gas  associated  with  three-body  recombination.  However  these 
terms  are  known  or  may  be  evaluated  in  graphical  form  as  functions  of  electron 
temperature.  Then,  it  is  always  possible  to  construct  a  relatively  simple  analytic 
expression  to  fit  the  known  graphical  form  of  these  functions.  This  permits  the 
energy  balance  equation  to  be  solved  for  the  electron  temperature  once  the  intensity 
of  the  electromagnetic  field  is  fixed.  When  the  electron  temperature  has  been 
determined,  the  electron  density  may  be  computed  from  Eq.  (72).  As  soon  as  the 
electron  temperature  and  density  have  been  determined,  the  dielectric  constant  (K) 
for  this  particular  point  in  the  plasma  may  be  found  by  using  Eq.  (67).  Finally,  on 
the  basis  of  the  model  of  a  plane-layered  plasma  medium,  the  electromagnetic  field 
distribution  in  the  nonlinear  plasma  slab  may  be  computed  step-by-step  by  utilising 
the  energy  balance  Eq.  (75),  the  expression  (72)  for  the  electron  density,  and  ex¬ 
pression  (67)  for  the  dielectric  constant  in  conjunction  with  Maxwell's  equations 
expressed  in  the  form  of  difference  equations. 


2. 4  Heat  Transport  Through  the  Electron  Gas  Due  to 
Conduction  and  Particle  Diffusion 


The  various  energy  loss  processes  of  the  electron  gas  due  to  heat  conduction 
and  particle  diffusion  have  been  carefully  investigated  hy  Anderson  and  Goldstein,®'  ** 
These  investigation*  were  devoted  mainly  to  the  study  of  cross  modulation 
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(Luxembourg)  effects  in  the  afterglow  of  pulsed  gaseous  discharge  plasmas.  Elec¬ 
tron-neutral  collision  frequency  changes  are  produced  by  heating  the  plasma  in  the 
afterglow  through  the  application  of  a  heating  microwave  pulse.  These  electron  col¬ 
lision  frequency  changes  are  observed  by  the  variations  in  attenuation  they  produce 

Q 

on  a  second  or  'wanted'  probing  electromagnetic  wave.  Anderson  and  Goldstein 

properly  note  that  an  energy  balance  equation  for  a  particular  region  of  the  electron 

gas  must  include  terms  which  account  for  energy  loss  due  to  heat  transport  through 

(1)  electron-electron  collisions  (Spitzer  and  Harm10  coefficient),  (2)  electron-neutral 

collisions,  and  (3)  particle  diffusion.  Hence,  in  addition  to  the  enumerated  electron 

gas  energy  gain  and  loss  terms  which  account  for  energy  exchange  between  the 

electromagnetic  field,  the  electron  gas,  and  the  neutral  and  ion  gas  [Eq.  (75)]  , 

energy  transport  from  one  part  of  the  electron  gas  to  another  due  to  gradients  in 

electron  temperature  and  electron  density  must  also  be  examined. 

□ 

Anderson  and  Goldstein  have  given  an  explicit  expression  for  the  heat  transfer 
coefficient  due  to  electron-neutral  and  electron-ion  collisions  which  was  first  derived 
by  Gould:4"* 
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where  ■'gpp  is  given  by  Eq.  (2)  and  v(v)  =  Z>eN  +  ^  .  The  heat  conductivity  co¬ 

efficient  due  to  electron-electron  collisions  (Kee>  is  larger  than  because  the 
interelectron  collision  frequency  is  high  and  the  Coulomb  forces  have  a  long  range. 
Spitzer  and  Harm10  have  derived  an  expression  for  the  heat  conduction  coefficient 
for  the  case  of  relatively  small  temperature  gradients  in  a  plasma: 
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where  c  5  (3k  T^m)1^2  ,  fn(q  c2)  "6  when  Ng  =  1012/cm3 

«T  *  0.  225  Te  =  3000*K 

Finally,  heat  transport  may  occur  due  to  the  existence  of  electron  density  gradients 
in  the  plasma.  This  term  assumes  the  form: 
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It  should  be  remembered  that  there  will  be  considerable  heat  transport  at  equi¬ 
librium  across  the  sharp  density  and  temperature  gradients  which  exist  at  the  shock 
front  and  at  the  boundary  layer.  This  heat  flow  will  proceed,  of  course,  even  in 
the  absence  of  an  electromagnetic  wave.  This  heat  flow  is  of  no  consequence  in  so  far 
as  the  determination  of  the  steady  state  electron  temperature  in  the  presence  of  the 
electromagnetic  field  is  concerned,  because  it  is  only  those  density  and  temperature 
gradients  established  by  the  electromagnetic  field  over  and  above  the  gradients  which 
had  originally  existed  at  equilibrium  which  will  be  responsible  for  loss  at  a  particular 
point  in  the  plasma  of  the  heat  deposited  by  the  electromagnetic-field.  These  con¬ 
siderations  lead  to  heat  transport  terms  which  may  be  written: 

(1)  Flow  of  heat  deposited  by  electromagnetic  wave  due  to  electron-neutral  and 
electron-ion  collisions 
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(2)  Flow  of  heat  deposited  by  electromagnetic  wave  due  to  electron-electron 
collisions 
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(3)  Flow  of  heat  deposited  by  electromagnetic  wave  due  to  density  gradients 
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The  nonlinear  reflection  and  transmission  coefficients  for  a  plane  wave  incident 
at  an  arbitrary  angle  upon  a  plane-layered  inhomogeneous  plasma  slab  having  the 
characteristics  of  the  re-entry  sheath  may  be  computed  by  determining  the  electro¬ 
magnetic  field  distribution  in  the  slab  by  integrating  Maxwell's  equations  by  a  step- 
by-step  procedure.  The  step-by-step  numerical  integration  of  the  field  equations 
entails  the  solution  of  an  energy  balance  equation  at  each  step.  The  most  general 
form  for  the  energy  balance  equation  at  an  arbitrary  point  in  the  plasma  slab  in¬ 
cludes  the  following  terms: 


[Energy  deposited  by  I 
e.  m.  field  in  plasma  > 
[jay  ohmic  heating  ! 


+ 


(“Energy  gain  of  electron”! 
J  gas  due  to  new  particles  l 
1  created  by  electron-  f 
Lneutral  impact  J 
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("Energy  loss  of  electron  gasl 
s due  to  dissociative  re-  ’ 
[combination  J 


Energy  loss  due  to  elastic 
electron-neutral  and  electron- 
ion  collisions  and  inelastic 
electron-neutral  collisions 
I  such  as  collisions  which  excite 
rotational,  vibrational,  and 
optical  levels  and  ionizing 
collisions  _ 
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J  Energy  loss  of  electron 
sgas  due  to  three-body 
[recombination  J 

[Energy  transport  due 
jto  electron-neutral  s 
I  and  electron-ion 
!  collisions  j 


[Energy  transport  dueH 
sto  electron-electron  t> 
[collisions  J 


(Energy  transport  duel 
]to  particle  diffusion  J 


If  the  last  three  energy  transport  terms  in  the  above  equation  are  negligible 
compared  to  the  terms  that  represent  heat  flow  between  the  electron  gas  and  neutral 
gas,  then  the  fact  that  the  neutral  gas  constitutes  an  infinite  heat  reservoir  implies 
that  there  is  no  coupling  of  energy  between  the  different  layers  of  the  plane-layered 
medium.  In  spite  of  the  fact  that  the  heat  transport  coefficients  may  be  large  (espe¬ 
cially  the  heat  conductivity  coefficient  due  to  electron-electron  collisions),  if  the 
electron  temperature  gradients  established  by  the  electromagnetic  field  are  not 
drastically  different  from  the  gradients  which  existed  in  the  absence  of  the  field, 
then  the  steady-state  electron  temperature  within  an  arbitrary  layer  of  the  plane¬ 
layered  medium  will  not  be  a  function  of  the  temperature  at  adjacent  points.  Rather, 
the  electron  temperature,  determined  from  the  energy  balance  Eq.  (75),  will  depend 
solely  upon  the  rate  of  heat  flow  into  the  electron  gas  due  to  ohmic  heating  and  the 
heat  flow  out  of  the  electron  gas  to  the  neutral  and  ion  gas  through  the  various  elas¬ 
tic  and  Inelastic  collision  processes. 

It  should  be  emphasized  that  the  effects  of  heat  flow  from  the  electron  gas 
through  the  boundary  layer  to  the  vehicle's  surface  have  been  entirely  neglected  in 
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the  present  analysis.  Since  most  of  the  streamtube  approximations  (Lin  and  Teare, 
Rotman  and  Meltz11)  neglect  diffusion  from  streamtube  to  streamtube  (constant 
enthalpy  along  a  streamtube)  in  the  process  of  calculating  electron  density  for  equi¬ 
librium  and  nonequilibrium  flow  regimes,  the  standard  techniques  for  computing 
electron  density  profiles  about  hypersonic  re-entry  vehicles  contain  no  inherent 
mechanism  which  can  account  for  heat  transport  across  the  ionized  layer.  Bond*4  has 
investigated  the  approach  to  equilibrium  ionization  in  argon  shocks.  His  analysis 
leads  to  the  prediction  of  charge  separation  across  the  shock  front  due  to  the  dif¬ 
ference  in  mobilities  of  positive  ions  and  electrons.  The  electron  density  gradient 
was  responsible  for  the  electron  diffusion.  This  effect  will  also  be  neglected  in 
this  report.  In  addition,  for  the  well-defined  shock  fronts  characteristic  of  blunt- 
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nosed  re-entry  vehicles  at  200,  000  (t  traveling  at  18,  000  ft /sec,  it  will  be  assumed 
that  the  gradients  of  electron  temperature  and  electron  density  in  the  boundary  layer 
and  at  the  shock  front  will  be  so  steep  compared  to  a  wavelength  (3  cm  for  X-band) 
that  they  may  be  represented  as  plane  interfaces.  By  examining  the  heat  transport 
coefficients  represented  by  Eqs.  (76),  (77),  (78),  together  with  the  electron  tempera¬ 
ture  and  electron  density  gradients  usually  encountered  behind  the  Bhock  of  a  blunt- 
nosed  re-entry  vehicle  at  200,  000  ft  traveling  at  18,  000  ft /sec  (see  Figures  18,  19, 
and  20  in  the  report  by  Rotman  and  Mettz11),  it  is  easily  seen  that  the  largest  heat 
transport  term  is  given  by  Kee(Te)  V2^  -  KeeCT)V2T,  representing  heat  conduc¬ 
tivity  due  to  electron-electron  collisions.  Assuming  that  the  electron  temperature 
and  electron  density  gradients  established  by  the  electromagnetic  field  are  not 
greater  than  ten  times  the  corresponding  gradients  that  prevailed  in  the  absence  of 
the  field,  it  may  be  demonstrated  that  the  largest  heat  transport  term  corresponding 
to  the  greatest  change  in  temperature  gradient  is  at  least  two  orders  of  magnitude 
less  than  the  term  representing  heat  loss  of  the  electron  gas  through  elastic  and 
inelastic  collisions. 

The  orders  of  magnitude  of  these  quantities  are  listed  below: 
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The  effects  of  heat  conduction  may  be  taken  into  account  as  a  second  order 
correction  to  the  electron  density  and  temperature  profile  across  the  layer  of  ionized 
gas.  The  procedure  involves  an  Iterative  technique  to  be  described  in  Section  5. 
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2.5  Diffusion  Effects 

It  should  be  noted  that  no  mass  motion  of  the  plasma  is  assumed  in  the  analysis. 
This  is  essentially  the  result  of  the  fact  that  the  plane  wave,  which  is  incident  upon 
the  plasma  slab,  has  an  infinite  extent.  Hence  there  is  no  flow  of  plasma  out  of  a 
finite  region  over  which  the  ohmic  heating  occurs.  The  electrons  which'  are  created 
on  the  spot  by  the  induced  changes  in  the  reaction  rates  are  assumed  not  to  diffuse 
across  the  layer  of  ionized  gas.  This  is  a  very  good  assumption  since,  within  the 
layer  of  ionized  gas,  strong  electrostatic  restoring  forceB  will  severely  inhibit  the 
flow  of  electrons  out  of  the  region  in  which  they  were  created.  Diffusion  is  charac¬ 
terized  by  strong  coupling  between  electron  and  ion  motions  (ambipolar),  since  the 
Debye  length  I Q  is  much  less  than  the  shock  layer  thickness.  The  positive  ion  motion 
is  completely  determined  by  the  structure  of  the  flow  field  of  the'  neutral  gas.  The 

Debye  length  for  electron  densities  of  10  elec/cm  and  electron  temperatures  of 
-3  -4 

10,  000*K  is  about  10  cm,  or  less  than  10  times  the  shock  thickness.  The  condi¬ 
tions  for  ambipolar  diffusion  are  satisfied,  which  implies  that  there  will  be  a  space 
charge  sheath  of  electrons  at  the  boundaries  of  the  plasma  slab  (shock  front  and 
boundary  layer).  It  will  be  assumed  that  this  electron  layer  buildup  at  the  boundaries 
has  a  negligible  effect  on  the  microwave  transmission  characteristics  of  the  re-entry 
plasma  sheath. 

3.  KINETKS  OF  HUB  TEMPERATURE  AIR  IN  THE 
PRESENCE  OF  AN  ELECTROMAGNETIC  FIELD 

There  have  been  some  efforts  devoted  to  the  calculation  of  enhanced  ionization 

of  the  re-entry  plasma  sheath  under  the  Influence  of  intense  microwave  radiation. 

8 1  82  83 

King  and  Gray,  King,  and  Sodha  have  computed  the  enhanced  degree  of  ioniza¬ 
tion  of  the  plasma  sheath  under  the  influence  of  high-power  electromagnetic  radia¬ 
tion  on  the  basis  of  the  Saha  equation.  The  Saha  equation  relates  the  electron  density 
in  a  plasma  at  the  (j  +  l)th  stage  of  ionization  to  the  temperature,  partition  function, 

and  ionization  potential  of  the  jth  stage  of  ionization.  The  usual  form  of  Saha's 
38 

equation  (Bachynski  et  al.  )  applies  only  to  the  case  of  thermodynamic  equilibrium, 

and  is  derived  under  the  assumption  that  the  distribution  function  for  particle  energies 

84 

is  Maxwellian  and  that  the  electron  temperature  equals  the  gas  temperature.  Dew  an 
has  generalized  the  Saha  equation  to  Include  any  steady-state  equilibrium  situation 
in  a  plasma.  Dewan  demonstrates  that  the  steady-state  equilibrium  condition  implies 
that  the  rate  of  collision  and  radiation  ionisation  processes  for  a  given  stage  of  ioni¬ 
zation  must  equal  the  rate  of  collision  and  radiative  recombination  processes  for 
that  stage.  This  leads  to  an  expression  relating  the  ion  densities  to  the  radiation 
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and  particle  energy  distribution  functions.  Dewan®*  has  specifically  pointed  out  that 
attempts  to  apply  the  usual  form  of  the  Saha  equation  to  physical  situations  where 
thermodynamic  equilibrium  does  not  pertain  can  lead  to  totally  erroneous  results. 
Dewan's  approach,  which  ignores  collisions  between  heavy  particles,  may  be  applied 
to  a  multicomponent  plasma  when  charge  exchange  is  negligible.  For  the  situation 
considered  in  the  present  analysis,  where  the  plasma  is  transparent  and  ionization 
due  to  recombination  radiation  is  negligible,  Dewan's  approach  reduces  to  the  steady- 
state  particle  conservation  equation  obtained  by  taking  the  zero-order  moment  of  the 
Boltzmann  equation.  An  accurate  picture  of  microwave-enhanced  ionization  of  the 
plasma  sheath  may  be  obtained  only  by  studying  the  numerous  particle-particle  and 
field  particle  interactions.  Such  a  detailed  study  will  then  reveal  the  relative  im¬ 
portance  of  the  competing  ionization  processes. 

The  complete  description  of  the  kinetics  of  the  appropriate  collisional  processes 
and  the  kinetics  of  chemical  reactions  occurring  in  a  multicomponent  plasma,  such 
as  shock-heated  air,  is  an  exceedingly  difficult  problem.  The  complexity  is  due  in 
part  to  the  simultaneous  occurrence  of  translational,  rotational,  and  vibrational 
relaxation  behind  the  shock  front,  vibrational  and  radiative  excitations,  dissociation 
and  ionization.  When  an  electromagnetic  field  of  high  intensity  is  impressed  upon 
an  ionized  flow  field,  the  rates  of  reaction  of  many  of  these  processes  will  be  changed. 
Nevertheless  the  problem  is  not  intractable,  primarily  because  many  of  these  proc¬ 
esses  can  be  assumed  to  occur  independently  of  the  others  to  a  high  degree  of  ac¬ 
curacy.  Translational  and  rotational  equilibrium  of  and  C>2  are  achieved  within 
a  distance  of  about  one  upstream  mean  free  path  behind  a  normal  shock  in  air  (1  .  = 

1  4 

0.  05  cm  at  an  altitude  of  about  200,  000  feet).  The  computations  of  Lin  and  Teare 

demonstrate  that  translational/vibrational  equilibrium  is  established  at  a  distance 

of  about  100  upstream  mean  free  paths .  The  processes  of  vibrational  relaxation 

and  dissociation  cannot  be  uncoupled,  because  molecules  in  excited  vibrational  states 

49  4 

dissociate  more  readily  than  unexcited  ones.  Lin  and  Teare's  investigations  of 
the  various  competing  ionization  processes  behind  the  shock  front  show  that  the 
relative  importance  of  these  processes  follows  the  order: 

(1)  Atom-atom  impact 

(2)  Photoionization 

(3)  Electron  impact  on  neutrals 

(4)  Molecule-molecule  and  atom-molecule  impact. 

This  order  is  preserved  for  most  of  the  region  behind  an  air  shock  at  0.  02  mm  Hg 
<200,  000  ft)  at  a  velocity  of  about  20,  000  ft /sec,  although  as  Lin  and  Teare  point 
out,  the  rapid  initial  rates  of  change  in  temperature  and  chemical  composition  be¬ 
hind  the  shock  front  can  cause  a  shift  in  the  relative  order  of  importance  of  the 
competing  ionization  processes.  Since  atom-atom  impact  is  the  dominant  electron 
production  mechanism,  the  rate  of  production  of  electrons  behind  the  shock  depends 
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upon  the  translational/vibrational  equilibrium  and  the  degree  of  dissociation. 

Lin  and  Teare  consider  the  following  processes  to  be  the  chief  chemical  reac¬ 
tions  for  high  temperature  air: 


(1) 

02  +  X  +  5. 

t  ev  -  O  +  O  +  X 

(2) 

N2+X+9.  8ev-N+N+X 

(3) 

NO  +  X  +  6. 

5ev  —  N  +  O  +  X 

(4) 

O  +  N2  +  3. 

3  ev  -  NO  +  N 

(5) 

N  +  o2  -  1. 

4  ev  —  NO  +  O 

(6) 

N2  +  o2  +  1 

.  9  ev  -  NO  +  NO 

Since  the  degree  of  ionization  for  shocks  up  to  30,  000  ft /sec  is  always  small  (less 
than  1  per  cent),  the  electron  production  and  recombination  processes  will  not  af¬ 
fect  any  of  the  neutral  reactions  (translational-vibrational  equilibration  and  dis¬ 
sociation).  Hence  Lin  and  Teare  were  able  to  compute  the  vibrational  and  transla¬ 
tional  temperatures  and  the  particle  densities  as  a  function  of  distance  behind  a  one¬ 
dimensional  shock  by  solving  the  chemical  rate  Eq.  (79)  simultaneously  with  the 
vibrational  energy  relaxation  equations  for  and  N.,  together  with  the  mass,  mo¬ 
mentum,  and  energy  conservation  equations.  A  knowledge  of  the  reaction  rate  co¬ 
efficients  for  the  numerous  electron  production  processes  enable  them  to  compute 
the  electron  and  ion  densities  as  a  function  of  distance  behind  the  shock  front. 

If  a  high-power  electromagnetic  wave  is  impressed  upon  the  ionized  flow  field, 
the  translational-vibrational  equilibration  process  and  the  dissociation  reactions 
(79)  will  not  be  influenced.  This  is  due  essentially  to  what  Lin  and  Teare  term 
'the  one-way  coupling1  between  the  dissociation  process  and  the  ionization  processes. 
Since  the  fractional  energy  loss  of  an  electron  per  collision  with  a  neutral  (G)  is 

_  9 

always  less  than  10  even  up  to  electron  temperatures  of  20,  000TC  and  the  degree 
of  ionization  is  0.  1  per  cent  or  less,  the  electron  gas  does  not  affect  the  energy 
balance  or  the  reaction  rates  of  the  atomic  and  molecular  processes  of  the  neutral 
constituents.  Only  the  ionization  and  electron  recombination  process  will  be  affected 
by  the  presence  of  the  high-power  electromagnetic  field.  For  example  an  intense 
microwave  field  will  alter  the  relative  importance  of  the  competing  ionization  proc¬ 
esses  by  Increasing  the  electron  temperature  so  that  ionization  by  electron  impact 
on  neutrals  will  become  more  important  them  photoionization,  and  even  as  important 
or  more  important  than  atom-atom  collisions.  This  would  be  most  pronounced  at 
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least  ten  upstream  mean  free  paths  behind  the  shock  front  where  the  electron  density 

4 

is  built  up  to  an  appreciable  level  (see  Figure  13,  Lin  and  Teare  ).  One  important 
point  that  should  be  emphasized  is  that,  in  the  present  analysis,  a  basic  considera¬ 
tion  rests  upon  the  assumption  that  the  electron  gas  is  in  thermal  equilibrium  with 
the  neutral  gas  in  the  absence  of  an  electromagnetic  field.  This  assumption  is 
almost  always  made  in  computations  of  the  electron  density  profiles  behind  a  shock 

front  (Ldn  and  Teare,4  Rotman  and  Meltz, 1 1  and  Bortner50).  Petschek  and  Byron51 
52 

and  Bond  have  investigated  the  onset  of  ionization  and  the  approach  to  equilibrium 

in  argon  shocks.  Here  the  dominant  electron  production  mechanism  appears  to  be 

electron-neutral  impact.  It  was  found  that  during  the  approach  to  equilibrium,  the 

electron  temperature  is  less  than  the  gas  temperature.  The  distance  over  which 

this  equilibration  takes  place  behind  a  shock  front  in  air  may  be  estimated  in  the 

following  manner.  The  fraction  of  the  average  energy  transferred  between  an  elec- 

_o 

tron  and  neutral  (G)  per  collision  is  of  the  order  10  Hence,  in  about  100  elec¬ 
tron-neutral  mean  free  paths  the  electrons  should  be  in  equilibrium  with  the  neutral 
particles.  For  a  5000*K  shock  corresponding  to  a  blunt-nosed  re-entry  vehicle  at 
200,  000  ft  traveling  at  20,000  ft  /sec,  this  distance  would  extend  approximately  1 
cm  behind  the  shock  front.  Since  this  distance  is  short  compared  to  a  wavelength 
(3  cm  at  X-band),  the  representation  of  the  gradients  and  discontinuities  at  the 
shock  front  by  plane  interfaces  is  a  good  approximation.  This  assumption  is  better 
justified  at  100, 000  feet. 

The  various  ionization  processes  which  may  occur  behind  a  shock  in  air  include: 

(1)  Ionization  by  atom-atom,  atom -molecule,  and  molecule -molecule  collisions 

(2)  Photoionization 

(3)  Electron  impact  on  neutrals  and  ions. 

Also,  charge  exchange  can  Influence  the  total  electron  density  because  electron-ion 
recombination  is  specie  dependent.  Electron  attachment  to  neutrals  is  another  pos¬ 
sible  loss  mechanism,  because  a  heavy  negative  ion  has  little  influence  on  the  micro- 
wave  conductivity  of  the  plasma. 

3. 1  Neutral  -  Neutral  Impact 

Perhaps  the  most  surprising  result  of  Lin  and  Teare' s  study  of  the  ionization 

processes  behind  an  air  shock  is  the  fact  that  atom-atom  impact  is  the  dominant 

electron  production  mechanism  for  shock  velocities  up  to  about  30,  000  ft /second. 

This  is  an  unexpected  result,  for  studies  of  ionization  processes  in  argon  shocks 
31  34 

(Petschek  and  Byron  and  Bond  )  Indicated  that  electron-atom  impact  was  the 
dominant  mechanism.  The  minimum  activation  energy  for  producing  ionization  py 
atom -atom  impact  corresponds  to  the  reaction: 

N+0+2.8ev  -  NO+  +  e. 
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The  concentration  of  atomic  apecles  will  not  reach  appreciable  levels  until  the  dis¬ 
sociation  of  Og  and  Ng  is  fairly  complete.  For  example,  for  a  shock  corresponding 
to  an  altitude  of  200,  000  ft  at  a  velocity  of  20,  000  ft/sec,  dissociation  of  Og  and  Ng 
have  attained  7S  per  cent  and  SO  per  cent  of  their  respective  final  equilibrium  level 
at  a  distance  of  ten  mean  free  paths  behind  the  shock  front.  At  this  distance  the 
translational  temperature  of  the  gas  particles  has  dropped  to  about  1  ev.  In  spite 
of  activation  energies  which  are  considerably  greater  than  1  ev,  ionization  by  atom- 
atom  impact  is  an  important  process  because  of  the  relatively  large  velocity- 
dependent  cross  sections.  The  cross  sections  are  large  because  of  the  crossing  of 
the  potential  energy  curves  representing  the  interaction  potential  of  the  atomic 
system  before  and  after  the  collision.  This  cross  section  may  be  determined  (see 
Bates  and  Massey  )  by  considering  the  fact  that  there  is  a  relatively  large  proba¬ 
bility  for  a  transition  to  occur  between  the  initial  and  final  states  of  the  atomic 
system  at  the  point  where  their  respective  potential  energy  curves  cross.  This 
probability  is  related  to  the  eigenfunctions  of  the  quasi -molecule  formed  by  the  col¬ 
liding  systems,  the  Hamiltonian  of  the  colliding  system,  and  the  lnternuclear  dls- 
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tance  at  which  the  curve  crossing  occurs  (Landau-Zener  formula  '  ). 

Lin  and  Teare4  have  considered  the  following  ionization  processes  due  to 
neutral-neutral  Impact: 


(1)  N  +  O  +  2.  8  ev  -  NO+  +  e 

(2)  N  +  N  +  5.  8  ev  -  Ng  +  e 

(3)  O  +  O  +  6.  9  ev  -*  Og  +  e 

(4)  X  +  O  +  13.  6  ev  -  X  +  0+  +  e 

(5)  X  +  N  +  14.  6  ev  -  X  +  N+  +  e 

(6)  N  +  Og  +  6.  5  ev  —  NO*  +  e 

(7)  O  +  NO  +  7.  9  ev  -  NOg  +  e 

(8)  N  +  NO  +  7.  9  ev  -  NgO+  +  e 

(9)  X  +  NO  +  9.  3  ev  -  X  +  NO+  +  e 

(10)  O  +  Ng  +  11.  2  ev  -  NgO+  +  e 

(11)  O  +  Og  +  11.  7  ev  -  Og  +  e 

(12)  X  +  Og  +  12.  1  ev  -  X  +  Og  +  e 

(13)  X  +  Ng  +  16.  6  ev  -  X  +  Ng  +  e 

(14)  Ng  +  Og  +  11.  2  ev  -  NO  +  NO+  +  e  . 


(80) 
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By  examining  the  rates  reported  in  the  literature,  taken  in  conjunction  with  theoreti¬ 
cal  considerations  on  the  potential  curve  crossing  of  atomic  systems  and  the  equi¬ 
librium  relationship,  Lin  and  Teare  have  selected  values  for  most  of  the  forward 
and  backward  rates  of  the  neutral-neutral  impact  ionization  reactions  (80).  By  and 
large,  the  values  chosen  for  these  rates  appear  to  be  reasonable,  except  for  reac¬ 
tions  (1),  (2),  and  (3)  of  Eq.  (80).  The  temperature  dependence  of  selected 

for  the  dissociative  recombination  rates  corresponding  to  the  backward  reactions 
(1),  (2),  and  (3)  requires  that  the  velocity  dependent  cross  section  assume  a  physi¬ 
cally  unrealistic  form.  Hence,  in  the  present  analysis,  the  rate  for  the  forward 

reaction  (1)  is  determined  by  making  a  least-squares  fit  of  the  functional  form 

- 1/2  45  4 

AT  exp  (-B/T)  to  the  experimental  data  of  Lin  (see  Figure  5,  Lin  and  Teare  ). 

This  procedure  is  described  in  detail  in  Section  2.  3.  2.  The  backward  rate  for  reac¬ 
tion  (1)  is  obtained  from  the  equilibrium  relationship  The  backward  rates  for  the 
dissociative  reactions  1,  2,  and  3  are  assumed  equal.  Then,  the  forward  rate  con¬ 
stants  for  the  associative  ionization  reactions  2  and  3  were  obtained  through  the 
equilibrium  relationship. 

In  the  present  investigation,  atom -molecule  and  molecule-molecule  impact 
ionization  processes  are  neglected  (except  reactions  4  and  5).  This  is  justified  on 
the  basis  of  Lin  and  Teare1  s  estimations  of  the  specific  ionization  rates  for  each  of 
the  various  ionizing  processes.  The  specific  ionization  rate  is  defined  as  the  abso¬ 
lute  value  of  the  time  rate  of  change  of  the  normalized  electron  density  multiplied 
by  the  quantity  (1 j/Ug),  where 

l ^  *  upstream  mean  free  path  *  0.  05  cm  at  an  altitude  of  200,  000  ft 
Ug  «  shock  velocity. 

4 

Figure  8b  from  the  report  of  Lin  and  Teare  indicates  that  for  a  shock  at 
200,  000  ft  traveling  at  20,  000  ft/sec,  atom-molecule  and  molecule -molecule  impact 
ionization  processes  are  quite  negligible  compared  with  ionization  by  atom-atom 
impact  for  most  of  the  region  behind  the  shock  front.  Consequently,  in  this  report, 
only  the  ionization  processes  due  to  neutral-neutral  impact  corresponding  to  reac¬ 
tions  (1),  (2),  (3),  (4),  and  (5)  of  Eq.  (80)  will  be  considered.  The  reaction  rate 
coefficients  for  the  forward  reactions  (4)  and  (5)  are  taken  from  the  estimates  made 

by  means  of  the  collision  theory  used  by  Bortner.  The  three-body  recombination 

36 

rates  for  the  reverse  reactions  (4)  and  (5)  are  based  upon  the  Thomson  theory 
generalized  to  the  case  of  a  multicomponent  plasma  in  which  the  electron  tempera¬ 
ture  is  not  equal  to  the  gas  temperature.  In  conclusion  it  might  be  mentioned  that 

39  40 

in  addition  to  the  forward  reactions  (4)  and  (5),  both  Bortner  and  Atallah  have 
estimated  that  the  following  reactions  are  capable  of  producing  0+  and  N+  at  quite 
a  rapid  rate: 
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(1) 

N  +  NO  -  0+  +  N2  +  e 

(81) 

(2) 

O  +  NO  -  N+  +  02  +  e 

• 

However,  Lin  and  Teare  have  neglected  these  reactions,  and  since  reliable  rate 
constants  for  these  reactions  are  not  available,  they  will  be  neglected  in  this  report. 


3.2  Photoionization 

The  rate  at  which  electrons  are  produced  per  unit  volume  for  a  particular  fluid 
element  behind  a  shock  due  to  photoionization  depends  upon  the  number  densities  of 
the  atomic  and  molecular  species,  their  photoionization  cross  sections,  and  the 
photon  number  density.  The  far  ultraviolet  appears  to  be  the  spectral  region  which 
is  most  effective  in  producing  photoionization  of  N2,  Og,  and  NO.  It  is  only  through 
a  consideration  of  excitation,  emission,  and  radiative  transport  phenomena  that  the 

spectral  photon  density  for  a  particular  fluid  element  may  be  determined.  This 

4 

difficult  problem  has  not  been  solved  by  Lin  and  Teare.  Instead  a  somewhat  crude 

estimate  has  been  made  of  the  electron  production  rate  due  to  photoionization.  On 
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the  basis  of  measurements  made  by  Camm  et  al.  ,  Hammerling  has  deduced  that 

the  major  contribution  to  the  far  ultraviolet  radiation  emanating  from  an  air  shock 

+  4* 

may  be  attributed  to  the  b'  '  2  —  X  '2_  transition  of  N„. 

U  g  4  £ 

This  fact  has  enabled  Lin  and  Teare  to  otrtain  an  expression  for  the  total  in¬ 
tensity  of  ultraviolet  radiation  per  unit  area  emanating  from  the  shock.  Photoioni¬ 
zation  cross  sections  of  N„,  O,,  and  NO  have  been  taken  from  the  data  of  Weissler 
55  **06  " 

et  al.,  and  Wainfan  et  al.  The  energy  of  the  photons  corresponding  to  the 

b'  '2+  -*  X  '2*  transition  of  N„  (12.  9  ev)  is  such  that  this  ultraviolet  radiation  is 
u  g  2. 

capable  of  ionizing  only  the  molecules  O,  and  NO.  These  facts  lead  to  the  following 

^  4 

approximate  expression  for  electron  production  due  to  nhotoionization  (Lin  and  Teare  ). 


dNe  »a9  [  (02>  +  (NO)]  {Np  L  exp(-E*  /kT* ) 
dt  “  8  (Tr  +  r e> 

where 

Tr  »  radiative  life  time  of  upper  electronic  state  of  N2 
«  4  X  10'9  sec 

Tg  =  collisional  de-excitation  time 


(82) 
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E  =  energy  difference  between  two  electronic  states  of  Ng 
=  12. 9  ev 

L  *  effective  length  of  a  particular  fluid  element 

T  *  =  an  effective  temperature  that  determines  the  population  of  the 
electronic  states  of  Ng. 

4 

By  using  Eq.  (82),  Lin  and  Teare  were  able  to  illustrate  the  effect  of  photo¬ 
ionization  on  the  electron  density  profile  behind  an  air  shock  corresponding  to  a 
velocity  of  about  20,  000  ft/sec  at  an  altitude  of  200,  000  feet.  Figure  12  of  the  re- 

4 

port  of  Lin  and  Teare  demonstrates  that  for  distances  behind  the  shock  front  greater 
than  one  upstream  mean  free  path,  photoionization  produces  a  negligible  contribution 
to  the  total  electron  density. 

3. 3  lonizatton  by  Electron  Impact 

Since  the  Maxwellian  distribution  function  has  a  long  tail  at  high  electron  veloci¬ 
ties,  ionization  produced  by  electron-neutral  impact  will  become  comparable  with 
associative  ionization  by  atom-atom  impact  when  the  electron  energy  becomes  equal 
to  or  greater  than  about  one-sixth  the  ionization  potential  of  a  neutral  constituent 
(NO  has  the  lowest  ionization  potential:  9.  25  ev).  Equation  <68b)  describes  the  rate 
at  which  electrons  are  created  by  electron-neutral  impact.  The  right-hand  side  of 
Eq.  (68b)  has  been  obtained  from  the  integral  expression  for  impact  ionization  by 

expanding  the  ionization  cross  section  in  powers  of  (e  -  W.),  and  neglecting  second- 

12  * 

order  terms  in  (CT/Wj).  Here,  c  =  -j  mv  is  the  electron  energy,  Wj  *  ionization 
potential  of  neutral  constituent,  and  kT  has  been  assumed  to  be  much  less  than  the 
ionization  potential  (kT  «  W^. 


dN  _  /w  , 

St  \#  *0*i0njv  dv 

“  Ne(wf~)  (irmkTe)  £Nj  ,^w  exp(lar)  ffion}  (t>€de 

"ion/V"  "  Jtd£ 
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’ion,  <Wij> 


where  <rion  j  (Wjj)  *  0  and  e|on  ^  (W^)  ■=  *a2.  This  result  was  first  obtained  by 
Lin.57  The  neutral  constituents  which  will  be  considered  in  this  analysis  include: 
Ng.  Og,  N,  O,  and  NO.  The  right-hand  side  of  Eq.  (68b)  should  give  a  fairly  ac¬ 
curate  estimation  of  the  electron  production  rate  due  to  electron-neutral  Impact 
for  electron  energies  up  to  about  2  ev  (20,  000*K). 

It  might  be  mentioned  that  the  mechanism  of  electron  production  by  electron 
impact  on  neutrals  in  excited  metastable  electronic  states  has  been  neglected. 


3. 4  Effects  of  Charge  Exchange 

The  effects  of  charge  exchange  on  the  total  electron  density  behind  an  air  shock 
will  be  manifested  indirectly.  The  presence  of  an  electromagnetic  field  in  the 
ionized  flow  field  does  not  influence  the  rates  of  reaction  for  the  charge  transfer 
processes,  for  the  gas  temperature  is  not  changed.  Charge  exchange  involves  an 
alteration  of  the  relative  concentrations  of  the  various  ionic  species.  The  electron 
density  will  be  affected  if  the  electron-ion  recombination  rates  are  different  for 
different  species.  The  somewhat  meagre  experimental  data  on  dissociative  recom¬ 
bination  rates  seem  to  indicate  that  the  recombination  coefficient  is  the  same  for 
the  species  NO+,  Nj,  and  Oj.  This  implies  that  the  electron  density  would  be  af¬ 
fected  only  by  charge  exchange  between  the  molecular  and  atomic  species: 

X+  +  YZ  -  YZ+  +  X  .  (83) 

This  type  of  reaction  (B 3) is  termed  asymmetrical  nonresonant  charge  transfer,  for 

4*  -20  2 

which  the  cross  section  at  2  ev  for  an  Li  -  Argon  system  is  about  10  cm 

IIQ 

(Fogel  et  al.  ).  Most  nonresonant  charge  transfer  processes  are  quite  Inefficient 
at  thermal  energies.  However  there  is  another  type  of  reaction,  termed  asym¬ 
metrical  resonant  charge  transfer,  which  may  have  a  larger  cross  section  at  ther¬ 
mal  energies.  An  example  of  asymmetric  resonant  charge  transfer  has  been  sug- 
59 

gested  by  Omholt: 

0+(2D)  +  N2  (X  'E*  v  ■  0)  -  0(3P)  +  Nj  <A2»u,  v  -  1)  .  (84) 

60 

According  to  Bates  and  Lynn,  it  is  improbable  that  this  reaction  has  a  large 
cross  section  el  thermal  velocities. 
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But,  as  Bates  has  indicated,  even  though  charge  transfer  is  rather  inefficient 
at  thermal  energies,  ion-atom  interchange  may  be  quite  rapid: 

X+  +  XZ  -  XY+  +  Z  (85) 

It  is  possible  that  the  rate  coefficient  for  ion-atom  interchange  may  be  as  large  as 
-9  3 

10  cm  /sec  at  thermal  energies,  provided  the  activation  energy  is  small. 

Hence,  it  would  seem  that  charge  exchange  could  influence  the  steady-state 
electron  density  if  the  following  typical  chain  of  reactions  occurred. 


(a)  0+  +  N2  -  NO+  +  N 

<b)  NO+  +  e  -  N  +  O 


(86) 


The  ion-atom  interchange  reaction  (86a)  could  be  fairly  rapid,  and  the  dissociative 
recombination  reaction  (86b)  proceeds  much  more  rapidly  than  the  three-body  re¬ 
combination  of  electrons  with  O4 .  Only  the  rate  of  the  ion-atom  interchange  reac¬ 
tion 


0+  +  Oz  -  02  +  O  (87) 
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has  been  reposed  in  the  literature.  Dickinson  and  Sayers  have  measured  a  rate 
- 1 1  3  - 18  2 

of  10  cm  / sec  at  300*K,  which  implies  a  cross  section  of  about  2  X  10  cm  . 

Lin  and  Teare*  have  listed  the  following  charge  exchange  reactions: 

(1)  N2  +  N  -  N2  +  N+  +  1.  0  ev 

(2)  N2  +  O  -  N2  +  0+  +  2.  0  ev 

(3)  N2  +  02  -  N2  +  02  +  3.  5  ev 

(4)  N2  +  NO  -  N2  +  NO+  +  6.  3  ev 

(5)  N+  +  O  -  N  +  O*  +  0.  9  ev 

(6)  N+  +  02  -  N  +02  +  2.  5  ev 

(7)  N+  +  NO  -  N  +  NO+  +  5.  3  ev 

(8)  0+  +  02  -  0+  02  +  1.  6  ev 
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(9)  0+  +  NO  -»  0+  NO+  +  4.  4  ev 

(10)  Og  +  NO  -  02  +  NO+  +  2.  8  ev 

These  reactions  are  all  of  the  asymmetric  nonresonant  type,  except  the  reverse 

reaction  (2),  which  may  be  of  the  asymmetric  resonant  type.  Lin  and  Teare  have 

taken  the  cross  section  for  all  the  exothermic  reactions  to  equal  3  X  10  cm  . 

This  is  a  couple  of  orders  of  magnitude  larger  than  the  asymmetric  nonresonant 

58  "f 

charge  transfer  cross  sections  measured  by  Fogel  et  al.  in  Li  -  Argon  systems 
at  2  ev.  However,  Lin  and  Teare4  have  investigated  the  effects  of  varying  the  exo¬ 
thermic  charge  exchange  cross  section  on  the  electron  density  and  positive  ion 
density  profiles  behind  an  air  shock  at  about  200,  000  ft,  corresponding  to  a  velocity 

A 

of  20,  000  ft/second.  Figure  1 1  of  the  report  by  Lin  and  Teare  demonstrates  that 
changes  in  the  charge  transfer  cross  sections  of  factors  of  100  in  either  direction 
produce  changes  of  about  10  per  cent  or  less  in  the  density  profiles  of  electrons  and 
all  positive  ions  except  .  For  distances  between  1  and  100  upstream  mean  free 
paths  behind  the  shock  front,  the  ionic  specie  02  exhibits  about  a  400  per  cent  change 
in  density  when  the  charge  exchange  cross  sections  are  varied  by  a  factor  of  100. 

Lin  and  Teare4  have  listed  the  following  ion-atom  interchange  reactions  and 
charge  rearrangement  reactions: 

(1)  N j  +  O  -•  NO  +  N+  -  2.  2  ev 

(2)  N*  +  O  -  NO+  +  N  +  3.  0  ev 

(3)  N2  +  0+  -  NO  +  N+  -  4.  2  ev 

(4)  Nz  +  0+  -  NO+  +  N  +  1.  1  ev 

(5)  No  +  O,  -»  NO  +  NO+  +  4.  5  ev 

*  (89) 

(8)  N2  +  O*  -  NO  +  NO+  +  0.  9  ev 

(7)  02  +  N  -  NO  +  0+  -  0.  2  ev 

(8)  Oj  +  N  -  NO+  +  O  +  4.  2  ev 

(9)  02  +  N+  -  NO  +  0+  +  2.  3  ev 

(10)  Oz  +  N+  -  NO+  +  O  +  6.  7  ev 

Although  no  experimental  data  on  the  rates  of  reactions  (89)  has  been  reported 
in  the  literature,  the  ion-atom  interchange  reaction  (87)  has  a  measured  cross 


71 


- 18  2 

section  of  about  10  cm  at  0.  03  ev.  This  result  may  be  interpreted  to  imply  that 
the  ion-atom  Interchange  processes  [  Eq.  (89)]  have  a  cross  section  which  is  one  or 
perhaps  two  orders  of  magnitude  greater  than  the  cross  section  for  asymmetric  non- 

4 

resonant  charge  transfer  [  Eq.  (88)] .  Lin  and  Teare  have  taken  the  cross  section 

o 

for  the  exothermic  charge  rearrangement  reactions  (89)  to  be  equal  to  about  *aQ  . 
The  cross  sections  for  the  reverse  reactions  (88)  and  (89)  were  computed  from  the 
equilibrium  constant  using  the  principle  of  detailed  balancing. 

Because  of;  (1)  the  apparent  insensitivity  of  the  electron  density  and  most  of 
the  ion  specie  densities  to  changes  in  the  charge  exchange  rates,  (2)  the  extreme 
paucity  of  experimental  data  and  the  uncertainty  in  the  charge  exchange  rates,  and 
(3)  the  enormous  increase  in  complexity  of  the  expressions  for  electron  and  ion 
densities,  the  effect  of  charge  exchange  will  not  be  included  in  the  present  analysis. 
However,  if  more  sufficient  experimental  data  on  charge  exchange  reactions  become 
available,  then  Eq.  (71)  for  the  rate  of  production  of  each  ion  species  should  be  modi¬ 
fied  to  include  these  effects; 

-T1  *  kf70.  1  <N><°>  +  kf70.  6  -  kr70.  1«°+>W.> 

+  kf88.  4  (N2HNO)  *  kr88.  4  ^2)(NO+)  +  kf88.  7(N+)<NO) 

*  kr88.  7  <N)(NO+)  +  kf88.  9  (°+)(NO)  *  kr88.  9 (0)(N0+> 

+  kf88.  10  *  kr88.  10  (°2^NO  *  +  kf89.  2*N2^°^ 

(90) 

*  kr89.  2  W»6+)  +  *f89.  4  <N2HO+>  "  kr89.  4(N)(NO+> 

+  kf89.  5  (N2)<°2)  '  kr89.  5  <NO)<NO+)  +  k^  6<0+)<N2) 

-  kr89  6  (NO) (NO+)  +  kfg9  8  (O+XN)  -  kr89  8(0)(NO+) 

+  kf89.  10  "  kr89.  10  ) 

The  expressions  for  the  rate  of  production  of  the  ionic  species  n1,  0+,  and  N+  in- 
eluding  the  effects  of  charge  exchange  are  as  complex  as  the  equation  for  NO  [  Eq. 
(90)]  .  They  may  be  obtained  from  inspection  of  Eqs.  (70),  (88),  and  (89). 

3, 5  Effect*  of  Electron  Attachment 

There  are  primarily  four  types  of  electron  attachment  processes; 
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(1)  Radiative  attachment 

e  +  x  -*  x"  +  hv 

(2)  Three-body  attachment 

e  +  x  +  y-x’  +  y 

(3)  Dielectronic  attachment 

e  +  x  -  (x‘)* 

(4)  Dissociative  attachment 

e+xy-*x’  +  y 

Of  the  various  constituents  of  high-temperature  air,  O,  0„,  and  NO  have  the  largest 
66  * 

electron  affinity  (Massey  ).  The  electron  affinity  of  O  is  1.  465  ev;  of  O^,  0.  43  ev; 
of  NO.  small  but  positive;  of  N2>  negative;  of  N,  negative.  In  air  at  thermal  ener¬ 
gies  of  about  0.  5  ev  corresponding  to  translational  temperatures  of  about  5000*K, 
none  of  these  negative  ions  will  be  stable.  Hence  the  effects  of  electron  attachment 
will  be  neglected. 

4.  THE  FIELD  DBnUBUTON  AND  THE  REFLECTION  AND 
TRANSlflSSBN  COEFFICIENTS  OF  AN 
INHOMOGENEOUS,  NONLINEAR  PLASMA  SLAB 

The  previous  sections  of  this  report  have  been  devoted  to  a  detailed  description 
of  the  method  for  computing  changes  in  the  electron  density,  electron  collision  fre¬ 
quency,  and  electron  temperature  which  have  been  brought  about  by  the  presence  of 
a  high-power  electromagnetic  wave  in  the  ionized  flow  field  surrounding  a  hypersonic 
re-entry  vehicle.  The  analysis  of  the  transmission  characteristics  of  the  re-entry 
sheath  under  the  Influence  of  high-intensity  microwave  radiation  is  based  upon  the 
model  of  a  plane  monochromatic  wave  incident  at  an  arbitrary  angle  upon  a  non- 
uniform  multicomponent  plasma  slab.  Part  of  the  complexity  of  the  problem  arises 
from  the  fact  that  the  reflection  and  transmission  coefficients  of  the  plasma  slab 
depend  upon  the  conductivity  of  the  ionized  gas,  which  is  a  function  of  the  electron 
density  and  collision  frequency.  These  parameters  depend  upon  the  electron  tem¬ 
perature  distribution  in  the  plasma  slab.  The  electron  temperature  is  determined, 
in  turn,  by  solving  the  energy  balance  equation  which  contains  the  local  field  as  a 
parameter.  The  technique  developed  for  the  computation  of  the  transmission  and 
reflection  coefficients  of  the  plasma  slab  consists  in  replacing  the  layer  by  a  stack 
of  homogeneous  sheets.  This  model  lends  Itself,  in  a  natural  way,  to  the  application 


73 


of  step-by-step  numerical  Integration  of  the  field  equations  expressed  In  the  form 
of  difference  equations.  By  assuming  a  value  for  the  amplitude  of  the  transmitted 
field,  and  taking  backward  differences  through  the  slab,  the  coupled  system  of 
equations  which  describe  the  conductivity,  electron  density,  energy  balance  of  the 
electron  gas,  and  the  electromagnetic  field  distribution  may  all  be  solved  simul¬ 
taneously.  This  method  is  general  enough  to  incorporate  the  inhomogeneities  of 
electron  density  and  collision  frequency  into  the  analysis.  The  method  is  not 
restricted  to  cases  where  the  permittivity  gradient  is  small  compared  to  a  wave¬ 
length,  such  as  the  convenient  WKB  asymptotic  solution.  It  is  found  that  for  the 
case  of  a  plane  wave  of  arbitrary  polarization,  incident  at  an  arbitrary  angle  upon 
the  slab,  the  wave  may  not  be  entirely  uncoupled  into  two  waves:  one  polarized 
perpendicular  and  a  second  parallel  to  the  plane  of  incidence.  This  is  because  the 
dielectric  constant  for  a  point  in  the  plasma  slab  is  a  function  of  the  amplitude  of 
the  total  local  field. 

A  homogeneous  plasma  slab  does  not  constitute  a  physically  realistic  model  for 
problems  concerned  with  the  nonlinear  interaction  of  microwave  radiation  with  the 
plasma  sheath.  This  is  because  the  microwave  power  absorbed  by  the  plasma  de- 
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pends  critically  on  the  ratio  (wp/u>).  Klein  et  al.  have  shown  that  a  maximum 

transfer  of  energy  between  an  electromagnetic  field  and  plasma  occurs  at  a  point 

2  2 

just  behind  the  plasma  resonant  density  (Np  *  u  m eQ/e  ).  This  implies,  of  course, 
that  the  interaction  of  the  microwave  field  with  an  inhomogeneous  plasma  results  in 

a  nonuniform  deposition  of  energy. 
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Bloembergen  and  Pershan  have  already  investigated  a  somewhat  similar  prob¬ 
lem  concerning  the  solution  to  Maxwell's  equations  at  a  plane  interface  between  a 
linear  and  nonlinear  medium.  However  their  analysis  proceeds  from  a  relationship 
between  the  nonlinear  atomic  properties  of  a  medium  and  a  time  dependent  suscepti¬ 
bility  tensor.  In  the  case  of  a  plasma  under  the  influence  of  high-power  electro¬ 
magnetic  radiation,  this  corresponds  to  the  situation  when  the  relaxation  time 
(r *  1  iGv)  for  energy  transfer  between  the  electron  gas  and  neutral  gas  is  much 
less  than  the  period  of  the  impressed  field  (ten  «  T  or  ui  «  Gv).  Then  the  electron 
temperature  will  not  reach  a  steady-state  value  independent  of  time.  Rather,  the 
electron  temperature  will  have  a  sinusoidal  time  dependence  at  the  second  harmonic 
of  the  frequency  of  the  impressed  field  (see  Ginzburg  and  Gurevich8).  The  electron 
density  will  also  be  time  dependent.  The  characteristic  time  for  electron  density 

variations  will  be  a  complicated  function  of  the  rates  of  electron  production  and 
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electron  recombination.  Molmud  has  considered  inducing  time-dependent  pertur¬ 
bations  of  the  electron  density  and  collision  frequency  in  the  D-region  of  the  iono¬ 
sphere  by  utilizing  high-powered  ground-based  radio  transmitters.  The  perturbations 
may  be  observed  by  the  changes  produced  in  the  absorption  of  a  second  or  wanted 
wave.  Here,  photoionization,  which  is  the  dominant  electron  production  mechanism, 
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is  assumed  to  proceed  at  a  rate  which  is  independent  of  the  electron  temperature. 
The  chief  electron  loss  mechanism  is  the  temperature  sensitive  process  of  attach¬ 
ment  to  Og. 

The  time-dependent  perturbation  treatment  of  nonlinear  media  of  Bloembergen 
67 

and  Pershan  predicts  harmonic  wave  generation  at  the  boundary  between  a  linear 

and  nonlinear  medium.  A  reflected  second  harmonic  wave  will  propagate  in  the 

linear  medium  in  the  same  direction  as  the  reflected  fundamental  wave,  whereas 

the  transmitted  harmonic  will  in  general  propagate  in  a  direction  different  than  the 

transmitted  fundamental.  The  matching  of  the  field  amplitudes  at  the  boundary  for 

a  wave  incident  at  an  arbitrary  angle  leads  to  generalizations  of  the  usual  Fresnel 
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formulas.  The  treatment  of  Bloembergen  and  Pershan  is  applicable  to  plane  wave 
transmission  and  reflection  from  magnetoactive  plasma-vacuum  interfaces  for 
incident  waves  whose  frequency  is  such  that  u  «  Gv.  In  the  present  analysis,  only 
waves  whose  frequency  satisfies  the  condition  u  »  Gv  will  be  considered,  so  that 
harmonic  generation  will  not  occur. 


4. 1  The  Field  Equations  (MRS) 

The  propagation  of  electromagnetic  waves  in  a  nonlinear  plasma  are  governed 
by  the  usual  equations  of  phenomenological  electrodynamics,  that  is.  Maxwell's 
equations: 


V  X  H  =  c  +  J  =  j  e  wE  +  J 

O  at  *  O 


(91) 


where  the  diamagnetic  effects  of  the  plasma  have  been  ignored.  The  E  and  H  fields 
are  assumed  to  be  proportional  to  e+^ut.  If  there  is  no  external  magnetic  field 
present  in  the  plasma,  the  macroscopic  current  density  may  be  related  to  the  electric 
field  through  the  scalar  conductivity; 

7  =  <xE  (92) 


For  the  nonlinear  medium  under  consideration,  the  conductivity  depends  upon  the 
square  of  the  electric  field  amplitude.  Equation  (65)  is  a  general  expression  for  the 
scalar  conductivity  in  terms  of  the  basic  plasma  parameters  such  as  electron- elec¬ 
tron  collision  frequencies,  the  Maxwellian  electron  velocity  distribution  function, 
and  the  electron  density.  It  is  more  convenient  to  work  with  the  complex  dielectric 
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constant  of  th«  plasma,  K,  defined  In  terms  of  the  conductivity  by  Eq.  (66).  Max¬ 
well's  equations  may  be  written  in  terms  of  the  complex  dielectric  constant: 


^xE  *  -  ju  WB 

O 

^  X  if  »  j€  wKE 

4  O 


(93) 


2  2 

For  the  range  of  parameters  of  interest  (u>  »  v  ,  u  »  v  ),  the  approximate  ex- 

C8 

pression  (67)  for  the  dielectric  constant  is  adequate. 


4. 2  Wave  Propagation  in  Nonlinear  Inhomogeneous  Media 


4.  2.  1  A  LINEAR,  HOMOGENEOUS  SEMI-INFINITE  PLASMA 
70 

Graf  and  Bachynsld  have  treated  the  problem  of  plane  wave  reflection  and 
transmission  from  a  plane  interface  between  vacuum  and  a  homogeneous  plasma. 

The  presence  of  the  electromagnetic  field  is  assumed  not  to  change  the  properties 
of  the  plasma  medium.  The  plasma  is  characterized  by  a  complex  dielectric  co¬ 
efficient: 

K  *  Kr  -  JKr 

The  vector  Helmholz  equation  which  governs  wave  propagation  in  the  plasma  medium 
is  obtained  by  taking  the  curl  of  the  first  Eq.  (93)  and  substituting  curl  3  from  the 
second  Eq.  (93): 


curl  curl  E  ■  grad  div  E  -  V2E 

=■  -  jwu  curl  3  *  k2  kE 
o  o 

By  utilizing  the  condition: 

f  ■  5  »  c  f  •  kE  *  e  K  f  ■  E  +  e  2  •  V"R  «  0 

0  0  0 

together  with  the  fact  that  the  plasma  is  uniform:  V  3  ■  0,  Eq.  (94)  becomes: 

V2E  +  k2  kE  «  0 

o 

where  kQ  *  w/c. 


(94) 


(95) 


(96) 
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Figure  1. 


Figure  1  illustrates  the  geometry  of  the  problem.  Since  the  field  vectors  will  not 
vary  as  a  function  of  x  (parallel  to  the  plane  of  incidence),  and  the  phase  velocity  in 
the  y  direction  will  be  the  same  in  both  media,  solutions  of  the  wave  Eq.  (96)  for 
each  component  may  be  written: 

Ej  *  Aj  exp  [  PZ  +  j  QZ  -  j  kQ  sin  9  y] 

where  i  *  x,  y,  *. 
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Graf  and  Bachynakl  have  obtained  expressions  for  the  real  and  imaginary  part  of 
the  propagation  constant  in  terms  of  the  dielectric  constant  and  the  angle  of  incidence: 


P  •  k0  {<Kr  -  sin2  9)2  +K2}  ^  \  (Kr  -  sin2  9)J  1/2 

Q  «  kQ  [^{(Kr  -  sin2  9)2  +  K2}  1/2  +  ^  (Kr  -  sin2  9)]  1/2 , 


(97) 


The  surfaces  of  constant  phase  do  not  coincide  with  the  surfaces  of  constant  ampli¬ 
tude  in  the  plasma  medium.  The  surfaces  of  constant  amplitude  are  planes  parallel 
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to  the  interface,  whereas  the  planes  of  constant  phase  are  at  an  angle  (tan'  1kQsin0/Q) 
with  respect  to  the  interface. 

If  Maxwell's  Eq.  (93)  is  written  out  in  component  form  in  Cartesian  coordinates, 

it  is  easily  seen  that  E  .  H„,  and  H  are  independent  of  E  ,  E  ,  and  H  .  The  first 
x  y  z  yzx 

set  of  field  components  corresponds  to  a  wave  with  the  electric  vector  perpendicular 
to  the  plane  of  incidence  and  is  termed  the  transverse  electric  mode  (TE).  The 
second  set  of  components  corresponds  to  a  wave  with  the  electric  vector  in  the  plane 
of  incidence  and  is  termed  the  transverse  magnetic  mode  fTM).  An  incident  wave 
polarized  at  an  arbitrary  angle  may  be  resolved  into  these  two  modes.  For  a  linear 
homogeneous  medium,  the  propagation  characteristics  of  one  of  these  modes  in  un¬ 
influenced  by  the  presence  of  the  other  mode.  However  these  two  modes  are  not 
uncoupled  for  a  nonlinear  medium,  and  are  uncoupled  only  in  special  cases  for 
inhomogeneous  media. 

70 

Graf  and  Bachynski  have  derived  generalized  Fresnel  equations  for  each  mode 
of  polarization  by  matching  the  tangential  components  of  E  and  it  and  the  normal 
component  of  5  at  the  interface.  Graf  and  Bachynski*70  have  derived  expressions 
for  the  Poynting  vector  in  the  plasma  for  both  modes  of  polarization.  It  is  demon¬ 
strated  that  for  the  TE  mode  the  instantaneous  Poynting  vector  changes  magnitude 
and  direction  during  one  cycle,  so  that  the  wave  is  not,  in  general,  plane.  The 
average  Poynting  vector  is  perpendicular  to  the  planes  of  constant  phase.  For  the 
TM  mode,  it  is  shown  that,  in  general,  neither  the  instantaneous  nor  the  average 
Poynting  vector  is  perpendicular  to  the  planes  of  constant  phase.  Only  for  the  case 
of  a  lossless  plasma  or  normal  Incidence  are  the  waves  in  the  plasma  plane. 

For  the  case  of  a  nonlinear  plasma  slab  with  a  permittivity  gradient  perpen¬ 
dicular  to  the  plane  interfaces,  wave  propagation  may  be  considered  by  imagining 
the  medium  to  be  composed  of  a  stack  of  linear  homogeneous  plasma  sheets.  The 
solutions  of-the  field  equation  for  a  particular  sheet  will  then  be  representable  in 
the  form: 


Ej  =  Aj  exp[  PZ  +  jQZ  -  jkQ  sin  9  y) 

The  essential  point  is  that  the  real  and  imaginary  parts  of  the  propagation  constant 
( P  +  JQ)  will  depend  upon  the  magnitude  of  the  dielectric  constant  (K)  at  the  point  Z 
in  the  plasma.  The  dielectric  constant  will  depend  both  on  the  nature  of  the  inhomo- 
genieties  and  on  the  square  of  the  local  field  amplitude. 


4.  2.  2  EXACT  SOLUTIONS  OF  THE  WAVE  EQUATION  FOR  A  ONE¬ 
DIMENSIONAL  INHOMOGENEOUS  LINEAR  PLASMA  MEDIUM 


A  rather  small  number  of  exact  solutions  of  the  wave  equation  have  been  obtained 
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for  normal  incidence  on  an  interface  when  the  dielectric  constant  varies  only  in  a 
direction  perpendicular  to  the  plane  Interface.  With  reference  to  Figure  1,  exact 
solutions  exist  only  when  the  dielectric  constant  K  varies  with  Z  in  certain  speci¬ 
fied  fashions.  The  solutions  are  usually  expressible  in  terms  of  tabulated  functions, 
and  are  particularly  useful  when  the  geometric  optics  approximation  VK/K  «  2ir/\ 

is  not  valid.  The  geometrical  optics  approximation  is  inappropriate  when  K  »  0. 

17 

Ginzburg  has  discussed  the  solutions  of  the  wave  equation  for  normal  incidence 
when  the  dielectric  coefficient  has  a  variation  in  one  direction.  The  form  of  the 
dielectric  constant  and  the  solutions  to  the  wave  equation  are  as  follows: 

(1)  A  linear  layer  without  absorption,  K  =  1  -Z/Z ,,  solutions  in  terms  of 
Bessel  functions  of  one-third  order,  or  Airy  functions. 

(2)  An  absorbing  linear  layer,  K  =  1  -Z/Z ^  -  j(a  +  0  Z/Zj),  solutions  in  terms 

of  Bessel  functions  of  one-third  order  with  complex  arguments. 

2  2  2  2 

(3)  A  parabolic  layer  without  absorption  K  *  1  -  u  (1  -  Z  /Z^),  where 

is  the  plasma  frequency  corresponding  to  the  maximum  electron  density  and  Zm 

is  the  half-width  of  the  layer,  solutions  in  terms  of  Weber  functions,  or  parabolic 
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cylindrical  functions  (see  Whittaker  and  Watson  ). 

_9 

(4)  A  layer  characterized  by  a  dielectric  constant  of  the  form  K  =  a(b  +  Z)  , 
solutions  in  terms  of  polynomials  in  Z. 

(5)  A  layer  characterized  by  a  dielectric  constant  of  the  form  K  =  a  +  b  e  • 

y7  —  J  y7  y7  . 9 

(1  +  e  )  +ce  (1  +  e  )  solutions  in  terms  of  hypergeometric  functions  (see 
Ref.  74). 

7 1 

Pappert  and  Plato  have  treated  the  cases  of  a  linear  profile,  a  parabolic 

profile,  and  a  cosine  profile.  A  cosine  profile  in  dielectric  constant  leads  to  solu- 
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tions  of  the  wave  equation  in  the  form  of  Mathieu  functions.  Taylor  has  demon¬ 
strated  that  the  solutions  to  the  wave  equation  may  be  represented  in  terms  of 

cylindrical  Bessel  functions  for  an  exponential  profile  and  for  a  profile  varying  as 
-2  73 

Z  Buchsbaum  has  obtained  solutions  of  the  wave  equation  in  the  form  of  Hankel 
functions  for  a  dielectric  constant  of  the  form 

K*K  (l+aZ)n  . 
o 

It  should  be  mentioned  that  there  exist  cases  where  an  exact  solution  to  the 
wave  equation  can  be  obtained  for  normal  incidence  and  that  these  cases  are  capable 
of  being  extended  to  arbitrary  angles  of  incidence  for  the  TE  mode  of  polarization. 
This  is  always  true  when  an  exact  solution  to  the  wave  equation  has  been  obtained 
co-responding  to  a  dielectric  constant  K(z),  such  that  it  contains  an  arbitrary  con¬ 
stant  that  does  not  depend  upon  Z.  The  proof  of  this  statement  rests  upon  the  follow¬ 
ing  three  facts: 
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(1)  For  the  TE  mode  of  polarization,  Vl5(Z)  =0,  so  that,  by  Eq.  (95), 

f- ■  S  *  0. 

(2)  The  component  of  the  phase  velocity  of  the  incident  wave  parallel  to  the 
interface  (k0  sin  9)  is  the  same  as  the  component  of  the  phase  velocity  of  the  trans¬ 
mitted  wave  parallel  to  the  interface. 

(3)  The  formal  structure  of  the  wave  Eq.  (96) 

2 

•4-S-  +  k2  (K  -  sln2e)  E  *  0 
dZ*  ° 


is  not  changed  for  the  TE  mode  as  the  angle  of  incidence  (9)  is  varied. 


4.  2.  3  APPROXIMATE  METHODS  FOR  OBTAINING  SOLUTIONS  TO  THE 
WAVE  EQUATION  FOR  A  ONE -DIMENSIONAL  INHOMOGENEOUS 
UNEAR  PLASMA  MEDIUM 

Most  treatments  of  electromagnetic  wave  propagation  in  inhomogeneous  media 
assume  that  the  variations  of  dielectric  constant  occur  along  only  one  coordinate 
direction.  Studies  of  electromagnetic  wave  transmission  and  reflection  character¬ 
istics  from  plane  media  whose  dielectric  constant  varies  in  a  direction  perpendicular 
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to  the  plane  interfaces  have  been  conducted  by  Schelkunoff,  Penico,  Stickler, 

Richmond, 14  Albini  and  Jahn,77  Nicoll  and  Basu,78  and  Brekhovsldkh. 79  Klein 
69  80 

et  al.  and  Safran  and  Meltz  have  applied  the  WKB  asymptotic  approximation  and 
a  numerical  integration  of  the  Riccati  equation  for  the  complex  Impedance  of  a  di¬ 
electric  slab  to  the  problem  of  computing  the  transmission  characteristics  of  the 
inhomogeneous  re-entry  plasma  sheath.  In  problems  involving  propagation  between 
two  media  with  differing  dielectric  constant,  the  surface  at  which  the  discontinuity 
occurs  usually  may  be  represented  as  a  plane,  cylindrical,  or  spherical  boundary. 

A  general  problem  would  involve  obtaining  solutions  to  the  wave  Eq.  (94)  to  compute 

physical  observables,  such  as  phase  shift,  absorption  coefficients,  reflection  and 
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transmission  coefficients.  Schelkunoff  has  pointed  out  that  for  the  general  case 
of  waves  propagating  in  an  inhomogeneous  medium,  such  that  the  direction  of  propa¬ 
gation  coincides  with  the  direction  of  the  gradient  of  dielectric  coefficient,  it  is  not 
always  possible  to  interpret  the  solutions  of  the  wave  Eq.  (94)  as  outgoing  or  incom¬ 
ing  disturbances .  However,  for  the  present  analysis,  in  which  a  plane-layered 
plasma  slab  is  bounded  on  both  sides  by  vacuum,  this  difficulty  will  not  arise  ex¬ 
cept  perhaps  when  the  incident  wave  is  polarized  in  the  plane  of  incidence  and  the 
dielectric  constant  K  *>  0.  This  would  correspond  to  the  physical  situation  where 
plasma  waves  could  be  excited  in  the  ionized  medium. 
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The  general  form  of  the  Maxwell  Eqa.  (93)  indicates  that  the  electric  field  com¬ 
ponents  are  coupled.  These  field  components  become  uncoupled  only  for  very  spe- 

12 

cial  coordinate  geometries.  Penico  has  investigated  the  uncoupling  of  the  field 

components  for  several  orthogonal  curvilinear  coordinate  systems.  The  choice  of 

coordinate  system  is  usually  dictated  by  the  nature  of  the  boundary  surface  when 

more  than  one  dielectric  medium  is  involved.  This  facilitates  the  matching  of  the 

12 

field  components  at  the  boundary.  Penico  has  emphasized  the  fact  that  the  usual 

techniques  for  obtaining  solutions  to  the  wave  Eq.  (94)  require  that  two  of  the  three 

unknown  field  components  be  eliminated  from  the  equation,  so  that  a  scalar  equation 

o 

for  a  single  field  component  may  be  derived.  Also,  the  operator  V  which  appears 
in  Eq.  (94)  may  be  interpreted  as  the  Laplacian  only  for  Cartesian  coordinates  (see 
Stratton,76  p.  48).  Penico12  claims  that  the  reduction  of  the  wave  Eq.  (94)  to  a 
scalar  equation  for  a  single  field  component  is  possible  in  general  only  when  the 
dielectric  constant  K  is  a  function  of  one  coordinate. 

Consider  the  oblique  incidence  of  a  monochromatic  plane  wave  on  a  plane-layered 
medium  which  is  inhomogeneous  in  a  direction  normal  to  the  interface  (Z-axis). 
Outside  the  layer,  for  Z  <  O  (see  Figure  1),  the  incident  wave  may  be  written: 


E  *  E0  exp[  -jk^  (sin  9  y  +  cos  9  Z))  .  (98) 

Since  the  incident  wave  is  independent  of  the  x-coordinate,  both  reflected  and  re¬ 
fracted  wave  will  be  independent  of  the  x-coordinate  for  a  medium  in  which  the 
dielectric  constant  is  independent  of  x.  Using  Gauss'  law,  Eq.  (95)  may  be  written 
in  Cartesian  coordinates  in  the  form: 


The  field  components  not  only  satisfy  the  wave  Eq.  (99)  they  must  also  satisfy  Max¬ 
well's  Eq.  (93),  which  assume  the  following  form  in  Cartesian  coordinates: 
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By  examining  Eq.  (100),  it  may  be  clearly  seen  that  the  vector  components  E  .  E. 

y  z 

H  form  one  independent  set;  whereas  the  components  E  ,  H  ,  and  H  form  a  second 

*  X  jr  Z 

independent  set.  The  first  set  corresponds  to  the  TM  mode  of  polarization;  the 
second,  the  TE  mode  of  polarization.  An  incident  transverse  electric  mode  will 
excite  only  a  TE  mode  in  the  plasma.  Similarly  an  incident  TM  wave  will  excite 
only  a  TM  mode  in  the  plasma  slab.  This  is  true  for  a  nonlinear  one -dimensional 
inhomogeneous  medium.  Moreover  for  a  linear  inhomogeneous  medium,  the  propa¬ 
gation  of  the  TE  wave  is  completely  Independent  of  the  presence  of  the  TM  wave,  so 
that  the  modes  are  entirely  uncoupled.  However,  for  a  nonlinear  medium,  whether 
it  is  inhomogeneous  or  not,  the  propagation  characteristics  of  the  TE  wave  depend 
upon  the  TM  mode,  for  the  dielectric  coefficient  is  a  function  of  the  square  of  the 
total  field  amplitude. 

The  phase  velocity  of  the  electromagnetic  wave  in  the  y-direction  (parallel  to 
the  interface)  is  the  same  in  vacuum  as  in  the  plasma  slab,  even  when  the  medium 
is  nonlinear  and  inhomogeneous.  This  is  merely  a  statement  of  Snell's  law  of 
refraction,  and  may  easily  be  proved  by  considering  the  nonlinear  inhomogeneous 
medium  as  the  limiting  case  of  a  medium  composed  of  a  stack  of  linear  homogeneous 
layers: 


k  sin  9  3  k  sin  6_ 
o  r 


where 
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8  *  angle  of  Incidence  and  0r  is  the  angle  of  refraction  at  the  first  interface.  One 
important  consequence  of  Snell's  Law  is  that: 

^-Bi--jko.in«Ei 

where  i  ■  x,y,  or  z.  Another  consequence  is  that  the  angle  of  incidence  equals  the 
angle  at  which  the  transmitted  wave  emerges  from  the  far  face  of  the  inhomogeneous 
nonlinear  slab. 

Richmond14  has  set  up  the  field  equations  in  the  form  of  difference  equations  so 

that  the  field  distribution  in  one -dimensional  Inhomogeneous  slabs  may  be  computed 

using  a  step-by-step  numerical  integration  technique.  By  assuming  the  value  for 

the  amplitude  of  the  transmitted  wave  on  the  source-free  Bide  of  the  slab,  the  field 

distribution  in  the  medium  may  be  calculated  by  taking  backward  differences.  For 

14 

the  TE  mode  of  polarization,  Richmond  has  utilized  the  wave  Eq.  (99)  for  the  x- 
component  of  the  electric  field.  To  obtain  the  solution  to  the  difference  equation, 
the  field  must  be  known  at  two  points.  The  amplitude  of  the  transmitted  wave  at  the 
source-free  aids  of  the  slab,  corresponding  to  the  point  Z  »  d,  is  assumed.  The 
field  at  the  location  of  the  first  backward  increment,  corresponding  to  Z  -  d-h,  is 
obtained  by  making  a  Taylor  series  expansion  of  the  field  about  the  point  z  -  d  in 
powers  of  the  increment  size  h.  The  difference  equation  form  of  the  wave  Eq.  (99) 
then  prescribes  the  field  amplitude  at  the  point  z  »  d-2h  in  terms  of  the  field  ampli¬ 
tudes  at  the  points  z  «  d  and  z  *  d-h.  Each  backward  step  corresponds  to  the  numeri¬ 
cal  integration  of  the  wave  equation  across  a  layer  whose  dielectric  constant  K(z) 
is  uniform  over  the  step-size  h.  The  boundary  conditions  for  the  field  components 
must  be  satisfied  at  each  transition  across  a  discontinuity  in  dielectric  constant. 

The  matching  of  the  field  components  at  each  boundary  of  the  plane -layered  medium 

presents  no  difficulty  for  the  TE  mode,  since  the  tangential  components  E  and  H 

*  y 

and  the  normal  component  H  are  all  continuous  across  a  discontinuity  in  dielectric 

constant  K.  The  treatment  for  the  TM  mode  is  somewhat  different.  Here  there  is 

a  component  of  the  electric  field  E„,  which  is  normal  to  the  plane  interfaces.  Since 

z 

the  normal  component  of  the  electric  displacement  must  be  continuous  across  a 
discontinuity  in  dielectric  constant,  the  boundary  condition  for  the  D  vector  would 
require  that  the  condition 

K(Z+)  Ez  *  K(Z.)  Ez 


(101) 
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be  satisfied  at  each  step  in  the  numerical  integration,  where  the  plus  subscript  re¬ 
fers  to  the  layer  just  to  the  right  of  the  discontinuity  and  the  minus  subscript  refers 
to  the  layer  at  the  left.  In  the  case  of  a  linear  inhomogeneous  medium,  this  added 
difficulty  may  be  circumvented  by  considering  only  the  field  component  tangent  to 
the  plane  interfaces,  E^.  Once  the  field  distribution  for  the  tangential  component 
Ey  is  determined,  there  is  information  sufficient  to  determine  the  total  reflection 
and  transmission  coefficients  of  the  slab.  This  is  because  the  angle  of  incidence 
must  equal  the  assumed  angle  of  transmission.  However,  for  a  nonlinear  medium, 
the  case  of  the  TM  mode  of  polarization  presents  special  difficulties.  It  is  impos¬ 
sible  to  find  the  field  distribution  for  the  component  alone,  because  all  the  elec¬ 
tric  field  components  are  coupled  through  the  dielectric  coefficient  K,  where  K  is  a 
function  of  all  the  components.  For  a  nonlinear  plane-layered  medium,  the  step-by- 
step  numerical  integration  must  be  carried  out  for  all  the  electric  field  components 
simultaneously,  with  the  condition  of  Eq.  (101)  imposed  at  each  dielectric  coefficient 
discontinuity  for  the  normal  electric  field  component  Ez_  This  involves  an  iteration 
at  each  boundary,  since  the  dielectric  coefficient  to  the  left  of  each  discontinuity, 
K(z_),  is  unknown.  First,  the  value  that  K(z_)  would  have  if  the  medium  were  linear 
is  assumed,  and  this  determines  Ez_.  The  fields  thus  obtained  are  used  to  compute 
a  new  value  for  K(z  ),  and  then  a  new  value  for  E„  is  obtained  from  Eq.  (101). 
Convergence  of  the  Iterative  procedure  should  be  rapid  since  the  nonlinear  departures 
from  the  linear  case  will  be  relatively  small. 


4.  2.  4  TECHNIQUE  FOR  COMPUTING  THE  FIELD  DISTRIBUTION, 

TRANSMISSION  AND  REFLECTION  COEFFICIENTS  OF  A  NONUNEAR 
ONE -DIMENSIONAL  INHOMOGENEOUS  PLASMA  SLAB 


The  procedure  used  to  compute  the  field  distribution  in  an  inhomogeneous  non- 
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linear  plaSma  slab  is  similar  to  the  method  used  by  Richmond  for  the  case  of 
linear  inhomogeneous  media.  A  major  modification  involves  the  direct  use  of  Max¬ 
well's  equations,  Eq.  (100),  to  find  the  field  distribution,  instead  of  working  with  the 
wave  Eq.  (99).  This  facilitates  the  numerical  integration,  for  a  Taylor  series 
expansion  of  the  field  amplitude  does  not  have  to  be  made  in  order  that  the  numerical 
integration  be  started. 

Equation  (100)  may  be  written  in  dimensionless  form: 


9% 

-§f  +  j  sin  9  uz  =  -J  vx 


J  Vv 


(102) 
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where 


sin  0  v  ‘  ■  Ku 
X  z 

ui  *  ^  Ei 

vi  *  ^  Hi 

T  » k0z 

Here  the  dielectric  constant  K  is  a  function  of  the  electron  density  Ng  and  effective 
collision  frequency  ><EFF.  For  the  present  analysis  the  approximate  expression 
(67)  for  K  is  adequate.  At  each  step  in  the  numerical  integration,  the  energy  balance 
equation,  Eq.  (75),  must  be  solved  for  the  electron  temperature  Te>  This  permits 
the  determination  of  the  electron  density  from  Eq.  (72)  and  the  effective  collision 
frequency  from  Eq.  (2).  The  field  amplitude  squared,  which  appears  in  the  ohmic 
heating  term  of  the  energy  balance  Eq.  (75),  is  the  square  of  the  total  field: 


E2  +  E2  +  E2 
x  y  z 


The  solution  is  started  by  assuming  values  for  the  amplitude  of  the  transmitted 

2  2  + 1/2  u“  1 

field  u,  the  angle  of  polarization  0  [tan  0  *  (u  +  u  )  x  ],  and  the  angle  of  trans- 

-i-i  y  z 

mission  0(0  *  tan  u  ai  )  at  the  source -free  side  of  the  slab  z  ■  d. ,  where  A  *  slab 
z  y  T 

thickness  and  plus  subscript  refers  to  a  location  on  the  side  of  a  dielectric  constant 
discontinuity  furthest  from  the  source.  At  a  point  just  inside  the  plasma  slab, 
corresponding  to  z  «  d_,  the  tangential  electric  and  magnetic  fields  u  ,  u  ,  v  ,  v  , 

-  x  y  x  y 

and  the  normal  magnetic  field  vz  have  the  same  magnitude  as  the  assumed  values 
at  the  point  z  *  d+.  If  there  is  a  normal  component  of  electric  field  uz>  its  value  at 
a  point  just  Inside  the  slab  (u^J  must  be  obtained  by  making  successive  approxima¬ 
tions  to  the  solution  of  Eq.  (101), 
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ud_  K(d->  *  ud+  K<V  '  ud+  * 

where  the  first  approximation  to  K(d  J  may  be  taken  as  the  value  of  the  dielectric 
constant  for  a  linear  medium.  If  there  is  no  component  of  the  electric  field  normal 
to  the  interfaces  (TE  mode),  there  is  no  change  in  the  value  of  the  amplitude  of  any 
of  the  field  components  as  a  discontinuity  in  dielectric  constant  is  traversed.  How¬ 
ever,  whether  the  mode  is  TE,  TM,  or  mixed,  the  calculation  of  the  dielectric  con¬ 
stant  at  the  point  *  *  d_  involves  the  solution  of  the  energy  balance  Eq.  (75)  for  the 
electron  temperature,  the  determination  of  the  electron  density  from  Eq.  (72),  and 
the  determination  of  the  effective  collision  frequency  from  Eq.  (2).  Once  a  self- 
consistent  set  of  values  have  been  obtained  for  the  field  components  and  dielectric 
constant  at  the  point  z  •  d-.  Maxwell's  equations  expressed  in  the  form  of  difference 
equations  give  a  prescription  for  finding  the  fields  at  the  point  z  *  (d-h)+,  where  h 
is  the  increment  size  used  in  the  difference  equation.  The  distance  h  corresponds 
to  a  plane  layer  thickness  over  which  the  electron  temperature,  density,  and  col¬ 
lision  frequency  are  assumed  to  remain  constant.  The  dimensionless  form  of  Max¬ 
well's  equations,  Eq.  (102),  expressed  as  difference  equations  are: 

uy(r-h)  -  uy(T)  +  Jh  (sin  0  -  ut(r) 

ux  (r-h)  «  ux(t)  -  Jh  vy(r) 

v  (r-h)  *  sin  8  u_(T-h)  (103) 

Z  X 

v  (T-h)  *  v  (t)  +  jh  [sin2e  -  KJ  u  (r) 

J  J  * 

vx<T-h)  *  vx(r)  +  JK  uy(r) 
u  (T-h)  •  -  sin  8  K'1  v  (t -h) . 

Z  X 

Equations  (103)  yield  the  values  pf  the  field  amplitudes  at  the  point  z  *  (d-h)+, 
once  the  fields  at  the  point  z  *  d_  are  known.  All  the  field  components  st  the  point 
z  •  (d-h)+  are  equal  to  the  fields  at  the  point  s  •  <d-h)_,  except  the  normal  component 
of  the  electric  Held,  uz,  which  is  again  determined  by  obtaining  a  self-consistent  set 
of  solutions  to  Eq.  (101): 


K(r-h).  ux(r-h)_  -  K(T-h)+  uf(T-h)+ 
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Eq.  (67)  for  K(r-h)_,  the  energy  balance  Eq.  (75),  the  electron  density  Eq.  (72), 
and  the  Eq,  (2)  for  effective  collision  frequency.  If  Maxwell's  equations  ( Eq.  (103)) 
are  combined  with  the  boundary  condition  for  the  normal  electric-field  component, 
the  results  may  be  represented  in  the  form: 


Uy(T'h) 
u  (r-h) 

Z 


uy(r)  +  jh[sin9-  j{|^  )Ul(r). 
K(r-h)+ 

*  WFW  u*<T'h)+ 


u8(r-h)+  -  -  vx(T-h> 

(104) 

vx<r-h)  *  vx(t)  +  j  K  uy(f) 

ux(f-h)  *  ux(r)  -  J  h  vy(r) 
v  <r-h)  *  sin  9  u  (r-h) 

vy(r-h)  -  vy(T)  +  Jh[  sin29  -  K]  ux(t) 


The  difference  equation  for  the  normal  component  of  the  electric  field  is  singular 

when  K  •  0.  This  situation  corresponds  to  the  excitation  of  plasma  waves  in  the 

medium  by  the  component  of  the  electric  field  normal  to  the  interface  fTM  mode). 

If  collisions  are  present  the  energy  associated  with  these  plasma  waves  will  result 
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in  heating  of  the  plasma  (see  Gins  burg  ). 

The  step-by-step  numerical  integration  is  performed  until  the  point  z  *  0_  is 
reached,  corresponding  to  the  face  of  the  slab  upon  which  the  plane  wave  is  incident. 
Once  the  entire  field  distribution  has  been  obtained,  the  reflection  and  transmission 
coefficients  for  each  mode  of  polarization  may  be  found  immediately  from  the  fol¬ 
lowing  relations: 


Transverse  Electric  Mode: 
sin  9  ( (Sim, I  +  ISj^l)  ■  -  Hx(o) 

cos  9  (|3jncl  *  l*S»l>  '  V°> 


(105) 
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Transverse  Magnetic  Mode: 
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where 

are  the  magnitudes  of  the  Incident  magnetic  and  electric  vectors,  |h_  |  and  |lLJ 
are  the  magnitudes  of  the  reflected  magnetic  and  electric  vectors,  |hL|  and  |  E_  | 
are  the  magnitudes  of  the  transmitted  magnetic  and  electric  vectors.  R^  and  Tj^ 
refer  to  the  reflection  and  transmission  coefficients  for  the  TE  inode  and  R 1 1  and 
T 1 1  refer  to  these  quantities  for  the  TM  mode.  The  angle  0  is  the  angle  which  was 
assumed  as  the  angle  of  transmission.  It  should  be  noted  that  the  angle  of  trans¬ 
mission  must  equal  the  angle  of  Incidence. 

In  this  analysis,  calculations  will  be  performed  only  for  normal  incidence.  In 
this  case,  Eq.  (105)  assumes  an  Indeterminate  form.  For  normal  incidence  there 
is  no  distinction  between  the  TE  mode  and  the  TM  mode.  The  following  relations 
should  be  used  for  computing  the  reflection  and  transmission  coefficients  of  a  non¬ 
linear  slab  for  normal  incidence: 


Einc '  n+siy  {3- 37  x  1<)5>  volt8/m 


Er  *  R  fl+RT(3'  37  X  105)  volta/m  . 

The  computations  may  be  performed  for  several  step  sizes  (h)  to  obtain  an  estimate 
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of  the  accuracy  of  the  solution.  Stickler  and  Richmond  have  discussed  the  de¬ 
pendence  of  the  accuracy  of  the  solution  on  the  Increment  size  for  the  case  of  trans¬ 
mission  through  plane-layered  dielectric  slabs. 

5.  CONCLUSIONS 

The  results  of  the  computations  for  the  reflection  and  transmission  coefficients 

of  a  nonlinear  plasma  slab  which  simulates  the  composition  of  high-temperature  air 

(5000TC)  will  be  made  available  In  a  future  report.  Order  of  magnitude  estimates 

o 

reveal  that  power  fluxes  of  the  order  100  to  1000  watts/cm  are  capable  of  produc¬ 
ing  factors  of  four  or  five  changes  In  electron  temperature.  This  will  result  in 
about  factors  of  two  change  in  collision  frequency  and  factors  of  two  or  more  changes 

in  electron  density.  From  the  graphs  of  attenuation  constant  vs  (u>  /w)  for  various 

31  P 

ratios  of  v/u  presented  by  Bachynski  et  al.,  it  may  be  inferred  that  these  changes 
will  Induce  changes  in  absorption  coefficient  of  factors  of  two  or  more  for  certain 
ranges  of  the  parameters  (u^/w)  and  ( v/u ).  These  effects  are  certainly  important 
from  a  re-entry  communications  system  point  of  view. 

Unfortunately  the  order  of  magnitude  changes  which  are  expected  to  result  from 
this  type  of  nonlinear  interaction  are  of  the  same  order  as  the  uncertainties  in  the 
rate  coefficients  for  electron-ion  recombination.  The  present  state  of  the  art  is 
such  that  even  for  calculations  of  electron  density  profiles  about  blunt-nosed  re¬ 
entry  vehicles  under  equilibrium  flow  regimes,  only  order  of  magnitude  estimates 
may  be  obtained.  Thus,  the  effects  of  the  nonlinear  interaction  of  microwave  radia¬ 
tion  with  ionized  flow  fields,  as  presented  in  this  report,  will  result  in  changes  of 
electron  density  which  are  of  the  same  order  of  magnitude  as  the  uncertainties  in  the 
electron  density.  However,  as  more  accurate  information  on  the  rate  coefficients 
for  the  many  processes  that  occur  in  high-temperature  air  become  available,  the 
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theory  presented  will  enable  accurate  predictions  to  be  made  on  changes  in  power 
transmitted  through  ionized  flow  fields  at  relatively  high  power  levels. 

The  present  analysis  has  been  based  upon  the  assumption  that  the  plasma  medium 

reaches  a  steady  state  under  the  Influence  of  a  plane  monochromatic  electromagnetic 

wave.  The  boundaries  are  assumed  to  be  so  steep  compared  to  a  wavelength  that 

they  may  be  represented  as  plane  Interfaces.  Figures  10a,  b,  and  c  of  Lin  and 
4 

Teare's  report  indicate  that  the  electron  density  has  essentially  reached  its  maxi¬ 
mum  value  one  hundred  upstream  mean  free  paths  behind  a  normal  one -dimensional 
shock  front  for  Mach  numbers  between  about  15  and  25.  At  an  altitude  of  200,  000 
ft,  one  upstream  mean  free  path  is  about  0.  05  cm  so  that  within  a  distance  of  5  cm 
the  electron  density  builds  up  to  about  5  X  10  particles /cm  .  If  this  region  over 
which  steep  electron  density  gradients  occur  is  less  than  about  a  quarter  of  a  wave¬ 
length,  they  may  be  represented  as  plane  interfaces  to  a  good  approximation.  Of 
course,  at  altitudes  below  200,  000  ft  the  gradients  are  even  steeper;  while  at  about 
250,  000  ft  the  shock  front  becomes  diffuse.  At  100,  000  ft,  one  hundred  upstream 
mean  free  paths  is  about  0.  1  cm.  The  one  hundred  mean  free  paths  criterion  for 
electron  density  buildup  applies  only  to  a  normal  one -dimensional  shock.  The  model 
considered  in  this  report  of  a  planar  slab  is  more  appropriate  for  regions  that  lie 
at  a  distance  of  more  than  one  nose  radii  behind  the  tip  of  a  hemisphere -cylinder 
body.  Naturally  the  plasma  sheath  in  this  region  is  cylindrical  in  shape,  but  a  planar 
geometry  will  give  an  accurate  description  of  the  electromagnetic  reflection  and 
transmission  characteristics  if  the  radius  of  the  cylinder  is  large  compared  to  a 
wavelength. 

Since  the  boundaries  defining  the  plasma  medium  at  altitudes  less  than  200, 000 
ft  are  very  steep  compared  to  wavelengths  for  X-band  radiation  (3  cm)  and  frequencies 
lower  than  10  kMcps,  the  boundaries  will  be  represented  as  plane  interfaces.  For 
an  abrupt  discontinuity  in  dielectric  constant  of  the  medium,  there  is  no  energy 
deposition  at  the  discontinuity,  only  reflection  or  transmission  of  electromagnetic 
radiation.  Hence  the  assumption  will  be  made  that  the  boundaries  of  the  plasma 

medium  are  not  disturbed  by  electromagnetic  radiation  at  X-band  frequencies  for 
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power  levels  up  to  1000  watts /cm  .  This  appears  to  contradict  the  result  of  the 
16 

investigations  of  King,  who  found  that  the  plasma-air  boundary  representing  a  shock 
front  will  move  toward  a  source  of  high -intensity  microwave  radiation.  This  study 
is  based  upon  the  WKB  solution  of  the  wave  equation  for  the  field  distribution  in  the 
plasma  medium,  and  is  valid  only  when  the  electron  density  gradients  are  small  com¬ 
pared  to  a  wavelength  (when  u>  >  Wp  and  d  in  Ne/dZ  «  1/X).  The  application  of  the 
WKB  approximation  to  problems  of  propagation  through  the  re-entry  plasma  sheath 
appears  to  be  restricted  to  altitudes  above  200,  000  ft  and  to  electromagnetic  radia¬ 
tion  wavelengths  of  about  1  cm  or  less.  The  analysis  as  presented  in  this  report 
pertains  to  electromagnetic  propagation  from  the  vehicle  out  through  the  re-entry 
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•heath,  whereas  King  is  concerned  primarily  with  propagation  in  the  opposite  sense. 
Thus,  insofar  as  propagation  out  through  the  re-entry  sheath  is  concerned,  the  inter¬ 
face  of  the  plasma  medium  defined  by  the  boundary  layer  will  not  physically  move 
because  of  the  proximity  of  the  vehicle  surface.  Since  the  power  level  of  the  electro¬ 
magnetic  wave  will  be  considerably  reduced  by  the  time  the  wave  has  traversed  the 
plasma  layer,  one  may  assume  that  the  boundary  of  the  plasma  medium  at  the  shock 
front  will  not  move.  One  of  the  basic  assumptions  in  this  report  is  that  the  electron 
density  and  temperature/ profiles  will  redistribute  themselves  behind  the  shock  front 
under  the  influence  of  high-power  electromagnetic  radiation  in  such  a  fashion  that 
no  mass  motion  of  electrons  will  result,  nor  will  the  boundaries  defining  the  plasma 

medium  be  seriously  perturbed. 

16 

King  has  found  that  an  ionization  front  will  move  toward  the  source  of  high- 
power  electromagnetic  radiation.  This  result  was  obtained  under  the  assumption 
that  the  WKB  approximation  adequately  represents  the  solution  to  the  wave  equation 
in  the  plasma  medium.  Thus,  King  has  obtained  the  result  of  a  moving  plasma 
boundary  precisely  for  the  case  where  the  electron  density  gradients  at  the  boundary 
vary  slowly  compared  to  a  wavelength.  The  model  King  has  chosen  to  represent  the 
interaction  of  high-power  radio-frequency  radiation  with  a  plasma  is  somewhat 
idealized,  in  that  a  full  description  of  the  many  microscopic  processes  is  not  taken 
into  account.  For  instance  only  three  parameters  are  used  to  characterize  the 
plasma  medium:  the  ionization  cross  section  (Oj),  the  total  cross  Bection  (a),  and 
the  ionization  potential  (Wj).  Not  only  is  such  a  model  restricted  to  plasmas  con¬ 
sisting  of  but  a  single  neutral  specie,  but  also  King  uses  the  same  parameters 
(a,  Cj,  and  Wj)  to  represent  the  medium  in  front  of  and  behind  the  ionization  profile. 
This  is  certainly  not  realistic  with  regard  to  air  shocks,  where  the  constituent  NO 
(which  has  a  low  ionization  potential)  exists  behind  the  shock  front  but  not  in  the  un¬ 
perturbed  air  ahead  of  the  shock  front.  Then,  too,  the  assumption  that  the  quantity 

&L  JV21-  B*1 

LW0  m(<A„2)  J 

represents  the  number  of  ion  pairs  created  per  unit  time  is  founded  on  the  supposition 
that  the  collision  frequency  for  momentum  transfer  v,  which  appears  in  the  formula 
for  the  conductivity,  is  identical  to  the  'effective'  collision  frequency  given  by  Eq. 

(2).  But  this  is  approximately  true  only  if  the  velocity -dependent  cross  section  for 
momentum  transfer  e(v)  is  not  too  strongly  dependent  on  velocity  and  if  the  isotropic 
part  of  the  electron  distribution  function  fQ  is  Maxwellian.  King  has  avoided  a  dis¬ 
cussion  of  the  conditions  under  which  the  isotropic  part  of  the  electron  distribution 
function  will  remain  Maxwellian  and  hence  there  is  no  direct  employment  of  the 
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electron  temperature  as  a  parameter.  This  was  possible  by  coupling  the  solution 
to  the  wave  equation  directly  to  the  particle  conservation  (continuity)  equation,  by¬ 
passing  any  consideration  of  an  energy  balance  equation.  However  the  conditions 
under  which  the  isotropic  part  of  the  electron  distribution  function  will  remain 
Maxwellian  under  the  influence  of  an  electromagnetic  field  (i<  »Gi<)  are  precisely 

the  same  conditions  that  lead  to  a  large  thermal  conductivity  coefficient  of  the  electron 
gas  (Keeh  Electron-electron  collisions  are  effective  in  transporting  heat  from  one 
part  of  the  electron  gas  to  another,  and  a  high  interelectron  collision  frequency 
implies  a  large  thermal  conductivity  coefficient. 

These  preliminary  facts  lead  to  the  important  observation  that  the  electromag¬ 
netic  energy  deposition  in  a  slowly  varying  electron  density  gradient  is  not  always 
69 

as  selective  as  Klein  et  al  claim,  because  the  pattern  of  electromagnetic  energy 
deposition  in  a  plasma  may  be  determined,  in  certain  important  instances,  by  the 
thermal  conductivity  of  the  electron  gas  in  addition  to  its  dependence  on  the  variation 
of  dielectric  constant  of  the  plasma  (due  to  electron  density  gradients).  Klein  et  al68 
have  found  that,  for  a  slowly  varying  electron  density  gradient,  a  maximum  transfer 
of  electromagnetic  energy  to  the  plasma  will  occur  just  anterior  to  the  resonant 
plasma  depth.  The  resonant  plasma  depth  is  the  depth  of  the  plasma  for  which  the 
electron  density  becomes  critical,  Nc  *  (u^^lmt^el  A  high  interelectron  col¬ 
lision  frequency  (i>  »  Gv,  a  condition  that  is  met  for  the  stagnation  region  of  the 

re-entry  plasma  sheath  and  many  radio-frequency  plasma  discharges  in  the  labora¬ 
tory)  will  produce  an  equilibration  within  the  entire  body  of  the  electron  gas  of  the 
electromagnetic  energy  deposited  at  any  point  on  the  electron  density  profile.  If 
vge  »  Gv,  the  time  for  energy  equilibration  within  the  electron  gas  Tge  -  l/vee 
is  much  faster  than  the  time  for  energy  transfer  between  electron  gas  and  neutral 
gas.  Electromagnetic  energy  which  is  most  efficiently  deposited  at  the  critical 
density  on  the  electron  density  profile  will  very  quickly  be  redistributed  within  the 
entire  electron  gas  behind  the  ionization  front.  The  effect  of  the  high  thermal 
conductivity  of  the  electron  gas,  together  with  the  rapid  electron-electron  thermal!  - 
zationtime,  has  not  been  included  in  King's  model.  It  is  expected  that  the  electron 
thermal  conductivity  effects  will  tend  to  reduce  the  local  rate  of  electron  growth 
at  the  resonant  plasma  depth,  so  that  the  picture  of  an  ionization  front  moving  toward 
the  source  of  electromagnetic  energy  will  be  modified.  It  is  suggested  that  the 
whole  phenomenon  of  moving  ionization  fronts,  including  the  effects  of  electron 
thermal  conductivity,  should  be  more  intensively  investigated  both  from  the  experi¬ 
mental  and  theoretical  point  of  view. 

It  should  be  stressed  at  this  point  that  the  important  effects  of  electron  gas 
thermal  conductivity  may  be  incorporated  into  the  scheme  presented  in  this  report. 

In  the  first  approximation  the  effects  of  the  thermal  conductivity  of  the  electron  gas 
are  neglected.  This  is  justified  if  the  electron  temperature  gradients  established 
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by  the  presence  of  the  electromagnetic  field  are  not  greater  than  one  hundred  times 
the  gradients  that  prevailed  in  the  absence  of  the  electromagnetic  field.  It  may  then 
be  easily  shown  that  the  terms  in  the  energy  balance  which  represent  heat  flow  from 
electron  gas  to  neutral  and  ion  gas  are  one  or  two  orders  of  magnitude  greater  than 
the  terms  which  represent  heat  transport  from  one  part  of  the  electron  gas  to  an¬ 
other  (when  i>  »  Gv,  the  term  ?  K  VT  is  dominant).  Hence,  in  the  first 
66  ce  © 

approximation  the  plasma  medium  may  be  divided  into  a  stack  of  homogeneous  slabs 
with  no  thermal  coupling  between  them.  Then,  the  steady-state  electron  density  and 
temperature  distribution  due  to  a  high-power  electromagnetic  wave  may  be  found  by 
using  the  energy  balance  equation,  the  equation  for  the  effective  dielectric  constant 
of  the  plasma  medium,  and  Maxwell's  equations  written  in  the  form  of  difference 
equations.  The  energy  balance  equation  is  solved  for  the  electron  temperature  at 
each  step  in  the  numerical  integration,  where  only  terms  that  represent  heat  flow 
from  electron  gas  to  neutral  and  ion  gas  are  Included.  Once  the  new  electron 
density  and  temperature  distribution  have  been  found  due  to  the  perturbing  electro¬ 
magnetic  field,  a  correction  to  the  temperature  distribution  due  to  thermal  conduc¬ 
tivity  may  be  taken  into  account  by  an  iterative  procedure.  Consider  a  thin  sheet 
of  plasma  within  the  slab,  located  at  the  position  Z  and  with  thickness  dZ.  The 
heat  flowing  into  the  face  at  Z  is  Jz  Kee(dT  J dZ)  dA,  whereas  the  heat  flowing 
into  the  face  at  Z  +  dZ  is  Jz+dz  Kee  (dTe/dZ)  dA.  By  using  the  divergence  theorem, 
the  net  heat  flowing  into  the  thin  sheet  may  be  written: 

j  dT  dK  /  dT  \2  d2^ 

HZ  Kee  ~HZ~  ("d Z  )  +  Kee 

Since  only  the  heat  flow  produced  by  the  perturbing  electromagnetic  field  should  be 
considered  in  the  energy  balance  equation,  the  quantity 


represents  the  heat  flow  into  a  thin  sheet  of  plasma  due  to  electron  thermal  conduc¬ 
tivity. 

After  finding  the  temperature  distribution  in  the  slab  by  neglecting  the  electron 
thermal  conductivity,  a  second  step-by-step  numerical  integration  of  the  field  equa¬ 
tions  may  be  performed,  where  now  the  term 
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is  added  to  the  energy  balance  equation  at  each  step  in  the  numerical  integration. 

This  procedure  will  give  a  new  electron  temperature  distribution  (and  a  new  electro¬ 
magnetic  field  distribution)  that  may  now  be  used  to  perform  a  third  step-by-step 
numerical  integration.  The  entire  procedure  may  be  repeated  until  a  completely 
self-consistent  set  of  solutions  to  Maxwell's  equations,  the  energy  balance  equation, 
and  the  equation  for  the  effective  dielectric  constant  of  the  medium  have  been  ob¬ 
tained. 

The  assumption  that  the  neutral  and  electron  gas  remain  stationary  during  the 

interval  of  high-power  interaction  is  reasonable,  so  that  there  will  be  little  cooling 

of  the  plasma  as  a  result  of  convective  flow.  The  characteristic  time  for  energy 
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equilibration  within  the  electron  gas  is  r  _  =  l/»»  »  10  sec  and  the  character- 

66  66  ^  g 

istic  time  for  energy  transfer  between  electron  gas  and  neutral  gas  is  rEN  “  10  sec 
at  T  s  5000TC  and  r_„  *  10*7  sec  at  T  =  20,  000TC  for  the  stagnation  region  of  an 
air  shock  at  200,  000  feet.  This  relaxation  time  for  energy  transfer  between  electron 

n  _Q 

gas  and  neutral  gas  will  become  as  fast  as  10”  sec  at  Tg  *  5000®K  and  10  sec  at 
Te  -  20,  000°K  for  shocks  at  100,  000  feet.  Assuming  a  characteristic  length  of  10  cm 
corresponding  to  the  width  of  an  electromagnetic  beam,  and  a  maximum  flow  velocity 
of  106cm/sec,  the  time  for  a  portion  of  the  plasma  sheath  to  flow  out  of  the  influ - 

_5 

ence  of  the  perturbing  beam  is  approximately  10  seconds.  Thus  the  very  important 
conclusion  has  been  reached  that  the  plasma  sheath  may  be  considered  as  a  stationary 
gas  when  considering  the  perturbing  influence  of  high-power  electromagnetic  radia¬ 
tion  on  the  ionized  flow  field,  especially  at  altitudes  below  200,  000  ft  and  at  power 

levels  that  may  raise  the  electron  temperature  to  20,  000TC  (corresponding  to  an 

-2 

electron  fractional  energy  loss  G  *  10  ). 

The  entire  analysis  has  been  carried  as  far  as  the  present  state  of  the  art 
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permits.  At  power  levels  higher  than  about  1000  watts/cm  ,  electron  temperatures 
will  be  brought  into  a  range  (above  30,  000*K)  where  the  cross  sections  for  elastic 
and  inelastic  electron-neutral  collisions  and  the  rate  coefficients  for  electron-ion 
recombination  are  unknown. 

A  basic  parameter  in  the  analysis,  the  relative  fractional  energy  loss  of  an 
electron  per  collision  (G)  is  known  for  most  constituents  only  up  to  electron  energies 
of  about  3  ev  (30,  000°K).  Of  course,  it  has  been  implicitly  assumed  during  this 
investigation  that  the  values  of  the  relative  fractional  energy  loss  of  an  electron  per 
collision  with  neutrals  which  have  been  experimentally  determined  for  cold  gases 
(T  »  300TC)  will  still  have  the  same  dependence  on  electron  temperature  for  high 
temperature  air  (T  *  5000*K).  In  other  words  it  has  been  assumed  that  G  is  only 

a  function  of  T  and  not  T.  This  is  similar  to  the  implicit  assumption  of  Shkarofsky 
32  ® 

et  al  that  the  collision  cross  sections  for  electron  scattering  of  the  neutral  con¬ 
stituents  is  a  function  only  of  T#  and  not  T.  The  gas  temperature  is  assumed  to 
have  an  effect  only  in  that  the  degree  of  dissociation  of  air  changes  with  gas 
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temperature,  so  that  the  effective  scattering  cross  section  or  effective  G  factor  of 

the  gas  mixture  must  be  determined  by  weighting  the  cross  sections  or  C  factors 

by  the  relative  concentrations  of  the  appropriate  neutral  constituents.  However  the 

th 

cross-section  cr  or  the  parameter  for  the  particular  j  1  neutral  constituent  is 
assumed  to  vary  only  with  Tg  and  not  T . 

That  this  assumption  may  not  be  entirely  correct  can  be  seen  from  the  following 
considerations.  A  homonuclear  diatomic  molecule  such  as  Ng  or  02  has  no  perman¬ 
ent  dipole  moment,  so  that  the  vibrational  levels  corresponding  to  each  electronic 
state  of  the  molecule  are  metastable  and  may  have  a  lifetime  aB  long  as  10  seconds. 

-  _Q 

The  rotational  states  have  a  lifetime  of  about  10  seconds.  When  the  translational 

temperature  behind  an  equilibrium  shock  in  air  is  about  T  «  5000®K  (0.  5  ev),  most 

of  the  0„  will  be  dissociated  but  there  will  be  an  appreciable  amount  of  N„.  Molecule- 
1  +*l  - 1 

molecule  and  atom-molecule  collisions  will  occur  at  the  rate  of  about  10  sec 

It  requires  only  one  or  two  molecule-molecule  collisions  to  achieve  translational- 

rotational  equipartition  and  about  100  molecule -molecule  collisions  to  achieve  trans- 

and  for  T  =  500D»K, 
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lational -vibrational  equipartition.  For  T  *  5000*K,  p/p„  -  10 
g  •  c  o 

Gv  w  10  sec*1.  At  about  one  hundred  upstream  mean  free  paths  behind  a  normal 
one -dimensional  shock  at  200,  000  ft,  corresponding  to  a  velocity  of  about  20,  000 
ft/ sec,  an  appreciable  fraction  of  the  molecules  (mainly  NO  and  N,.)  will  be  in  ex- 

no  * 


cited  vibrational  states.  Harries 
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found  that  the  lowest  vibrational  state  of  N2  is 


0.  29  ev.  Since  NO  is  a  nonsymmetric  diatomic  molecule,  it  has  a  permanent  dipole 

_  O 

moment;  thus  the  lifetime  of  a  vibrational  level  is  expected  to  be  about  10  sec. 


The  time  between  electron-neutral  inelastic  collisions  is  TgN  »  1/Gv  »  10 


-6 


sec, 

so  that  electron  collisions  with  the  molecule  NO  can  result  in  collisions  which  ex¬ 
cite  at  least  the  lowest  vibrational  level  of  NO.  However  the  lifetimes  of  the  vibra- 
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tlonal  levels  of  the  homonuclear  molecules  N2  and  Og  are  about  10  aec  and  10 
sec  respectively,  which  are  longer  than  1/Gv  »  10"6  seconds.  Hence  electron 
energy  loss  under  impact  with  the  molecules  N2  and  02  can  occur  only  by  excitation 
of  the  rotational  levels  of  the  associated  metastable  vibrational  levels  of  these 
molecules,  and  not  by  excitation  of  the  vibrational  levels  themselves  because  the 
second  vibrational  level  would  not  be  available  for  excitation. 

For  cold  gases  (T  »  300*K)  the  electron  energy  loss  per  collision  with  a  homo¬ 
nuclear  diatomic  gas  is  independent  of  pressure  and  depends  only  on  electron  tempera¬ 
ture.  Since  crosB -modulation  measurements  of  Gr  are  usually  made  in  gases  that 
are  less  than  0.  1  per  cent  ionized,  relatively  few  molecules  are  involved  in  inelastic 
encounters  and  there  is  always  a  plentiful  supply  of  molecules  which  will  not  be  in 
metastable  levels.  This  implies  that,  in  the  case  of  cold  (T  K  300°K)  homonuclear 
diatomic  gases,  cross-modulation  measurements  of  Gr  for  Tg  »  5000°K  would  not 
depend  upon  whether  rEN  =  1/Gv  is  less  or  greater  than  tvib> 
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This  analysis  implies  that  whereas  for  a  cold  gas  (T  »  300TC,  Tg  »  5000*K)  the 
measurements  of  Gv  are  dependent  only  on  Te>  for  a  high  temperature  (T  >  5000TC, 

Tg  3  5000*K)  homonuclear  gas  (N2)  the  fractional  energy  loss  of  an  electron  depends 
both  on  electron  temperature  and  whether  rEN  3  1/Gv  is  much  greater  or  much 
less  than  Tyjg.  For  N2  and  02  the  following  relations  should  apply: 

When  ten  3  1/Gv  «  TyIB  (High  Pressure) 

GfT  3  5000*K,  T  «  5000TO  <  G(T  3  300TC,  T  3  5000*K) 
e  e 

When  ten  3  1/Gv  »  Tyjg  (Low  Pressure) 

GO1  «  5000*K,  Te  3  5000TO  3  GfT  3  300*K.  T#  3  5000TC)  . 

For  a  diatomic  molecule  with  a  permanent  dipole  moment  such  as  NO, 

G(T  3  5000*14,  T  3  5000*K)  3  GfT  3  300*K,  T  3  5000*K) 
e  e 

regardless  of  the  value  of  Gv  (rEN  3  1/Gv  is  almost  always  longer  than  TyiB(NO)  * 

10" 8  sec  except  at  very  high  pressures). 

7 

Massey  and  Burhop  contend  that  the  G  values  of  02  are  somewhat  higher 
than  the  G  values  of  N_  for  cold  gases  at  T  »  0.  5  ev  because  of  low-lying  electronic 
states  of  the  oxygen  molecule  which  may  be  excited  under  electron  impact.  The 
low-lying  electronic  states  of  02  (aA^  at  0.  8  ev  and  b  E*  at  1.  3  ev)  are  metastable 
with  enormous  lifetimes  of  the  order  of  10  seconds.  At  T  3  3000*K  (O,  is  fairly 
well  dissociated  at  T  3  4000*K  and  p/pQ  *  10  ),  about  10  per  cent  of  the  02  molecules 

will  be  in  the  excited  metastable  low-lying  electronic  states.  Only  the  rotational 
levels  associated  with  the  lowest  vibrational  state  of  the  metastable  electronic  state 
of  02  will  be  available  for  excitation  under  electron  impact.  Hence  it  is  expected 
that  at  high  gas  temperatures  (T  3  3000*10,  high  pressures  tEjj  3  1/Gv 
and  for  **  0.  5  ev,  the  G  values  for  02  will  decrease  from  the  cold  gas  G  values 
corresponding  to  the  same  electron  temperature.  The  interesting  point  is  that  the 
G  values  for  02  should  exhibit  a  more  pronounced  decrease  than  the  G  values  for 
as  the  gas  temperature  is  raised.  The  G  values  for  NO  should  not  show  much  of  a 
change  as  the  gas  temperature  is  raised.  A  final  point  should  be  added  about  the 
effects  of  electron  attachment  on  the  measured  G  values  of  electronegative  gases 
(NO  and  02).  Although  electron  attachment  to  neutrals  constitutes  a  loss  of  energy 
so  far  as  the  electron  gas  is  concerned,  the  rates  of  electron  attachment  (three  body) 
and  detachment  are  so  slow  compared  with  the  electron  energy  loss  due  to  rotational, 
i  vibrational,  and  electronic  excitation  that  electron  attachment  does  not  contribute 
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to  the  relaxation  time  rEN  *  HGv  as  measured  in  a  cross -modulation  experiment. 
These  interesting  suggestlons*lead  to  the  possibility  of  a  whole  new  host  of  experi¬ 
mental  investigations  of  the  G  values  of  diatomic  gases  at  high  gas  temperatures 
fT  >  3000*K). 

It  is  hoped  that  the  type  of  approach  presented  in  this  report  will  bring  about  a 
better  understanding  of  the  difficult  problem  of  breakdown  of  missile  antennas  at 
high-power  levels.  The  model  of  a  plane-wave  incident  at  arbitrary  angle  upon  a 
nonlinear  one -dimensional  inhomogeneous  multicomponent  plasma  slab  may  also  be 
applicable  to  the  interpretation  of  ionospheric  radio-scattering  data  performed  at 
high-power  levels.  Molmud68  has  already  suggested  using  high-powered  ground- 
based  radio  transmitters,  such  as  the  Puerto  Rican  Arecibo  facility,  to  perturb 
the  D  region  of  the  ionosphere.  Induced  changes  in  electron  density  and  collision 
frequency  may  be  observed  by  the  changes  in  reflection  coefficient  they  produce  on 
the  high-power  transmitting  wave  or  on  a  second  wanted  wave.  From  an  examina¬ 
tion  of  the  shape  of  pulse  returns,  it  may  be  possible  to  gain  information  on  the 
dominant  rate  processes  for  electron  attachment  and  recombination. 


*  ■ 

The  author  is  indebted  to  Mr.  Michael  Mentzoni  for  several  stimulating  discussions. 
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